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SUMMARY 


A  new  analytical  approach  to  the  problem  of  the  accurate  calibration  of  metric 
cameras  is  developed  and  specific  applications  are  reported.  The  method  permits  an 
Indefinitely  large  number  of  frames  from  a  given  camera  to  be  reduced  simultaneously, 
yet  efficiently,  to  produce  common  parameters  of  the  inner  cone  for  all  frames,  as  well 
as  independent  elements  of  exterior  orientation  for  each  frame.  Because  control  points  may 
be  exercised  repeatedly  in  a  common  reduction,  very  large  sets  of  well  distributed  residuals 
can  be  generated  from  a  relatively  small  set  of  control  points  (in  principle,  a  complete  and 
accurate  calibration  could  be  performed  from  as  few  as  three  control  points).  From  such 
large  sets  of  residuals  empirical  functions  can  be  derived  to  account  for  persistent,  unmod- 
elbd  systematic  error.  In  addition,  it  becomes  feasible  to  establish  the  variation  in  the 
accuracies  of  the  radial  and  tangential  components  of  plate  coordinates  throughout  the 
format  and  thus  to  establish  appropriate  empirical  weighting  functions  for  subsequent 
applications.  The  method  is  of  universal  applicability  and  encompasses  (a)  calibration' 
from  aerial  photographs,  (b)  calibrations  from  stellar  photographs,  and  (c)  calibrations 
from  multicollimator  photographs.  Applied  to  stellar  calibrations,  the  method  leads  to 
improved  accuracies  and  convenience  by  completely  doing  away  with  conventional  re¬ 
quirements  for  precise  timing  of  the  shutter  and  for  stability  of  the  camera  throughout 
successive  exposures.  Applied  to  multicollimator  calibrations,  the  method  has  far  reaching 
implications  concerning  both  the  use  of  existing  multicollimators  and  the  design  of  future 
multicolfimators. 


FOREWORD 

This  is  Part  1  of  o  final  report  prepared  by  Dyane  C.  Brown  of  DBA  Systems,  Inc., 
under  Contract  DA-44-009-AMC-1457(X)/  project  number  4A001 45001 A52C,  task  01, 
subtask  006.  ETL  project  engineer  was  Lawrence  A.  Gambino,  Chief,  Advanced 
Technology  Division,  Computer  Sciences  Laboratory. 
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ADVANCED  METHODS  FOR  THE 
CALIBRATION  OF  METRIC  cameras 

By 

Duane  C.  Brcwn 


1.  INTRODUCTION 

It  is  generally  appreciated  that  calibration  is  most  meaningful  and  effective  when 
it  is  performed  on  a  system  in  its  operating  environment.  Despite  the  acknowledged  soundness 
of  this  principle,  the  calibration  of  aerial  mapping  cameras  has  almost  universally  been  rele¬ 
gated  to  the  laboratory.  Laboratory  methods  were  long  ago  found  to  be  inadequate  for  the 
stringent  calibration  required  for  ballistic  cameras  and  have  been  supplanted  by  more  power¬ 
ful,  and  more  comprehensive  methods  of  stellar  calibration.  The  outstanding  success  of  Stella* 
calibration  with  ballistic  cameras,  some  of  which  employ  the  same  type  of  lenses  as  certain 
aerial  cameras,  has  naturally  led  to  suggestions  that  the  stellar  technique  also  be  applied  to 
aerial  mapping  cameras.  This  has  been  doitu  experimentally  to  a  limited  extent.  However, 
despite  its  potential  advantages,  the  stellar  method  does  indeed  share  with  laboratory  methods 
Hie  shortcoming  of  failing  to  simulate  the  typical  operational  utilization  and  environment  of 
mc  pping  cameras.  It  follows  that  for  the  stellar  method  to  become  fully  acceptable  for  cali¬ 
bration  of  mapping  cameras,  it  must  Le  shown  to  yield  results  equivalent  to  those  obtainable 
from  definitive  operational  calibrations. 

In  mid  1966,  D.  Brown  Associates,  Inc.  (DBA)  received  a  contract  from  the  Army’s 
Engineering  Test  Laboratories  (ETL)  to  implement  on  unsolicited  DBA  proposal  to  develop  an 
advanced  procedure  for  the  definitive  metric  calibration  of  photogrammetric  cameras  in  their 
operating  environment.  Although  also  applicable  to  stellar  calibrations,  the  proposed  technique, 


called  SMAC  (Simultaneous  Multiframe  Analytical  Calibration),  was  especially  directed 
at  the  probbm  of  calibrating  aerial  cameras  frcm  photography  taken  ov  ,f  a  targeted 
calibration  range.  Here  it  was  recognized  that  in  many  situations,  any  one  frame  would 
be  unlikely  to  contain  enough  target  images  to  permit  a  definitive  calibration.  Even  with 
highly  dense  control,  it  was  appreciated  that  such  factors  as  residual  uncompensated  film 
deformation  and  refractive  effects  of  air  turbulence  or  shock  waves  could  conceivably 
compromise  the  results  from  a  single  frame.  Accordingly,  it  was  considered  to  be  desirable 
to  be  able  to  merge  observations  from  a  large  number  of  frames  in  order  to  arrive  at  a 
comprehensive  result.  Inasmuch  as  each  frame  would  have  a  unique  set  of  elements  of 
exterior  orientation  (corresponding  to  the  instantaneous  position  and  angular  orientation  of 
the  cameras),  such  an  approach,  if  rigorously  applied,  would  entail  the  introduction  of  six 
fresh  unknowns  for  each  frame  carried  in  the  reduction.  Seemingly,  then,  the  reduction 
for  a  moderately  large  number  of  frames  could  well  become  computationally  prohibitive. 

The  SMAC  reduction  overcomes  such  difficulties  by  exploiting  certain  patterned  characteristics 
of  the  general  system  of  norma!  equations  to  produce  a  practical  and  fully  rigorous  solution  to 
the  problem  no  matter  how  many  frames  are  carried  in  the  reduction.  By  virtue  of  this,  the 
computational  requirements  of  a  SMAC  reduction  increase  only  linearly  with  the  number  of 
frames  exercised  rather  than  as  the  cube  of  the  number  as  in  a  'brute  force'  reduction.  The 
number  of  frames  that  can  be  accommodated  in  a  simultaneous  adjustment  is  therefore  without 
set  limit. 

In  this  report  we  shall  outline  the  development  of  the  SMAC  reduction  and  shall  present 
the  results  of  its  application  to  the  calibration  of  the  Air  Force  KC-6A  aerial  mapping  camera. 
We  shall  also  discuss  results  of  other  specific  applications  and  shall  indicate  how  the  method  can 
be  applied  to  laboratory  calibrations.  As  we  shall  see,  exceptionally  accurate  and  meaningful 
calibrations  are  mode  possible  by  the  application  of  SMAC  to  sets  of  photographs  taken  over 
anoerial  calibration  range.  Such  calibrations  con  be  particularly  useful  as  standards  for 
evaluating  the  adequacy  of  laboratory  and  stellar  methods  of  calibration. 


2. 


THEORETICAL  DEVELOPMENT 


2.1  Symmetric  Radio!  Distortion 

The  SMAC  method  of  calibration  is  a  direct  and  natural  outgrowth  of  our  earlier 
work  (Brown,  1956,  1957,  1958,  1964,  1965;  Brown,  Davis,  Johnson,  1964).  In  Brown 
(1956)  the  stellar  method  of  calibration  of  symmetric  radial  distortion  was  placed  on  a 
rigorous  footing  through  the  extension  of  the  basic  plate  reduction  developed  by  Schmid 
(1953).  This  extension  exploited  the  result  from  optical  ray  tracing  that  the  radial  distortion 
6r  of  a  perfectly  centered  lens  can  be  expressed  as  an  odd  powered  series  of  the  form 

(1)  6r  -  Kji»+Kgi*+K,r’+... 

in  which 

r  =  (x3  +  y3)^  =  radial  distance 

x,y=  coordinates  of  image  referred  to  principal  point  as  origin. 

By  virtue  of  the  fact  that  the  x,y  components  of  distortion  can  be  expressed  as 


it  was  shown  that  coefficients  of  distortion  (K; ,  Ks,  Ka , . . .)  could  be  introduced  directly 
into  the  projective  equations  and  that  these  parameters  could  be  sharply  recovered  simultaneously 
with  the  elements  of  orientation  in  a  rigorous  least  squares  adjustment.  This  reduction  marked 
the  beginning  of  *he  strictly  analytical  approach  to  camera  calibration.  To  be  sure,  stellar 
techniques  had  previously  been  used  for  camera  calibration,  but  they  were  dichotomous  in 
character.  That  is,  the  elements  of  orientation  were  determined  from  an  initial  reduction 
Involving  a  carefully  balanced  selection  of  control,  but  otherwise  ignoring  the  existence  of 
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radial  distortion.  From  the  resulting  residual  vectors,  the  distortion  function  was  then 
established  by  graphical  means  based  on  a  plot  of  the  radial  components  of  the  residuals 
versus  radial  distance . 

The  method  developed  in  Brown  (1956)  immediately  became  standard  for  Air  Force 
ballistic  cameras.  Although  the  method  constituted  a  significant  advance,  it  was  found 
with  some  cameras  to  be  impossible  to  achieve  an  rms  error  of  plate  residuals  of  much  better 
than  four  microns  when  stars  throughout  the  entire  plate  format  were  carried  in  the  reduction. 

By  contrast,  when  stars  were  limited  to  areas  of  about  one  third  the  diameter  of  the  format, 
no  difficulty  was  encountered  in  achieving  rms  errors  within  the  2  to  3  micron  range  considered 
to  characterize  random  plate  measuring  errors.  This  clearly  meant  that  for  such  cameras  a 
significant  unmodelled  systematic  error  was  affecting  the  results.  Although  the  source  of  the 
systematic  error  was  not  known  at  the  time  (plate  measuring  comparators  were  suspected,  but 
were  later  vindicated),  adverse  effects  were  minimized  in  trajectory  reductions  by  the  expedient 
of  performing  independent  plate  reductions  for  two  or  three  slightly  overlapping  regions 
covering  the  trajectory  on  each  plate. 


2.2  Decentering  Distortion 

It  was  not  until  1962  that  the  underlying  difficulty  was  resolved  in  an  investigation 
reported  in  Brown  (1964)  and  later  extended  in  Brown  (1965).  The  unmodelled  systematic 
error  was  identified  as  decentering  di'tortion,  i.e.,  the  distortion  resulting  from  the  imperfect 
centration  of  lens  elements.  This  type  of  distortion  had  been  considered  earlier  iri  terms  of 
the  thin  prism  model  as  described  in  both  the  second  arid  third  editions  cf  the  Manual  of 
Photogrammetry.  However,  it  was  not  possible  to  reconcile  the  observed  systematic  components 
of  residuals  from  plate  reductions  with  this  model.  The  reason  was  uncovered  in  Brown  ( 1964) 
by  means  of  analytical  ray  tracing  through  a  thin  prism  situated  in  front  of  a  hypothetically 
perfect  lens.  This  revealed  the  exact  analytical  expression  for  the  thin  prism  model  and  showed 
that  the  generally  adopted  description  of  the  thin  prism  effect  was  incomplete.  The  full  state¬ 
ment  of  the  model  was  found  to  demand  an  asymmetric  radial  component  in  addition  to  the 
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previous!/  accepted  tangential  component.  As  thus  revised,  the  thin  prism  model  was  found 
to  be  consistent  with  observed  residuals  when  effects  of  projective  compensations  were  duly 
taken  into  account . 

In  Brown  (1965),  the  thin  prism  model  was  compared  with  a  mo  derived  by  Conrady 
(1917).  Although  the  Conrady  model  was  derived  by  means  of  direct  analytical  ray  tracing 
through  a  decentered  lens  (rather  than  from  arguments  of  analogical  equivalence  underlying 
the  thin  prism  model),  it  had  been  neglected  by  the  photogrammetric  community.  The  thin 
prism  and  Conrady  models  were  found  to  be  in  agreement  with  regard  to  the  magnitude  of  the 
tangential  component  of  decentering  distortion,  but  the  radial  components  of  the  two  models 
were  found  to  differ  by  a  factor  of  three.  It  was  shown  that  this  discrepancy  could  be  resolved 
by  a  projective  transformation  insofar  as  first  order  effects  were  concerned.  However,  the  pro- 
jective  equivalence  of  the  two  models  was  found  no  lunger  to  hold  when  higher  order  effects 
were  of  significance.  Accordingly,  ri.c  rhin  prism  model  was  supplanted  in  Brown’s  plate  re¬ 
duction  by  an  extended  form  of  the  more  widely  applicable  model  of  Conrady. 

In  terms  of  radial  and  tangential  components,  Conrady's  model  assumes  the  form 

Ar  *  3  P„  sin  (<p-  tf>0) 

(3) 

At  *  Pt  cos(;p-  <pc) 
in  which 

Pr  '  -  Ji  r2  +  J2r*  +  +  ...  -  profile  function  of  tangential  distortion 

<p  =  angle  between  positive  x  axis  and  radius  vector  to  point  x,y  =  arc  sin  (x/r) 

=  arc  cos  (y/r) 

< po  =  angle  between  positive  x  axis  and  axis  of  maximum  tangential  distortion. 

In  the  thin  prism  model,  the  factor  of  three  in  the  expression  for  Ar  »s  replaced  by  unity. 

In  both  models,  the  quantities  requiring  determination  in  the  process  of  calibration  e  the 
phase  angle  <p0  and  the  coefficients  of  the  profile  function  J],  J2,  •  •  •  -  In  Brown  (1%5) 
it  is  shown  that  Conrady’s  model  can  be  expressed  in  terms  of  x  ar.d  y  components  as 
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Ax  =  [  Pj  (r2  +  2x2)  +  2  P2  x  y]  [  1  +  P3  r2  +  . . .  ] 

(4) 

Ay  =  [2^xy  +  P2  (r2+ 2)^]  [1  +  P3r2  +  ...] 

in  which  the  new  coefficients  Pj,  P2,  P3  are  defined  by 

Pj  -  -J|sin0o 
P2  =  Jj  cos  0O 

<5)  -  V Jt 

P4  =  yj, 

This  formulation  has  the  advantage  of  being  a  linear  expression  in  the  leading  coefficients 
P|,  P2  when  the  higher  order  coefficients  P3#  P4  are  zero. 

When  decentering  exists,  the  principal  point  of  autocollsmation  does  not  coincide 
with  the  principal  point  of  photogrammetry.  With  mapping  lenses  this  separation  is  likely 
to  amount  to  only  a  few  tens  of  microns  and  is  of  no  metric  consequence.  With  objectives 
having  very  long  focal  lengths  (such  as  those  of  astronomical  telescopes)  the  separation  could 
amount  to  several  millimeters  and  thus  assume  a  degree  of  significance.  The  means  for  handling 
this  situation  is  discussed  in  Brown  (1965). 

When  the  model  for  decentering  distortion  was  incorporated  into  the  plate  reduction, 
it  became  possible  with  good  consistency  to  achieve  rms  errors  within  the  desired  two  to 
three  micron  range  without  resorting  to  piecewise  reductions. 

2.3  Single  Frame  Calibration 
2-3.1  introduction 

We  shall  define  a  frame  to  consist  of  an  exposure  or  a  series  of  exposures  for  which 
the  position  and  orientation  of  the  camera  may  be  considered  to  be  invariant.  In  this  context. 
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a  stellar  plate  containing  several  different  exposures  may  be  regarded  as  consisting  of 
either 


(°)  a  single  frame  if  the  orientation  is  considered  to  remain  unchanged  through¬ 
out  the  series  of  exposures; 

or 

(b)  several  frames  it  the  orientation  is  considered  to  change  (however  slightly) 
from  one  exposure  (or  group  of  exposures)  to  the  next. 

In  the  case  of  aerial  photographs,  each  exposure  would  ordinarily  constitute  an  independent 
frame. 

In  this  section  we  shall  formulate  the  plate  reduction  for  a  single  frame  under  the 
assumption  that  the  photographed  control  (whether  stellar  or  ground  control)  is  flawless. 

Such  a  reduction  is  a  special  case  of  the  advanced  plate  reduction  formulated  in  Brown  (1964). 
As  v/e  shall  presently  see,  the  single  frame  reducrion  generates  the  basic  building  blocks  needed 
in  the  Simultaneous  Multiframe  Analytico!  Calibration  SMAC. 


2.3*2  Observational  Equations 

The  projective  equations  resulting  from  on  undistorted  central  projection  may  be  written 
as  (Brown,  1957) 

AX  +  Bu  +  Cv 
X  "  XP  C  DA  +  EH  +  FiT 

A*X  +  9u  +  Cv 
y"y-  =  c  DTTFpTTr 


=  elements  of  interior  orientation, 

_  orientation  matrix,  elements  of  which  are  functions  of  three  independent 
angles  at  W,  K  referred  to  arbitrary  X,Y,Z  frame  in  object  space, 

=  X,Y,  Z  direction  cosines  of  roy  joining  corresponding  image  and  object  points. 


in  which 


*p,  yp*  c 

A  B  C” 
A*  Bf  C 
DE  F_ 

A,p,v 
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When  points  in  object  space  are  at  a  sensibly  finite  distance,  the  direction  cosines  may  be 
expressed  as 


X  =  (X-XC)/R 

(7)  p  =  (Y-YC)/R 
v  =  (Z-ZC)/R 

C  jC  c 

in  which  X,Y,Z  and  X  ,  Y  ,  Z  are  obiect  space  coordinates  of  the  point  and  the  center  of 
projection,  respectively,  and  R  is  the  distance  between  these  points,  i.e., 

(8)  R  =  [(X-XC)2+  (Y-YC)2+  (Z-ZC)23*. 

When  stellar  or  other  strictly  directional  control  is  employed,  the  direction  cosines  may  be 
expressed  as 

X  =  sin  Of*  cos  to* 

(9)  p  =  cos  a*  cos  to* 

V  -  sin  to* 

in  which  the  angles  a* ,  to*  are  measured  in  the  same  sense  as  the  pair  of  angles  a,  todefining 
the  direction  of  the  camera  axis. 


We  shall  assume  that  the  direction  cosines  X,  p,j/  have  been  corrected  for  atmospheric 
refraction  by  means  of  standard  formulas.  If  we  then  let  x®,  y®  denote  the  observed  plate 
coordinates,  the  left  hand  sides  of  the  projective  equations  (6)  can  be  replaced  by 


(10) 


x  -Xp  =  x  +  vx  +  x  (Kj  r2  +  l<2  r*  +  K3r*  +  . . .) 

+  CP!  (r2  +  2x2}  +2P2Ry][l  +Par2  +  P4  r*  +  ...  ] 

y-yp  =  y  +  vy  +  y  (Kt  r2  +  K2  r4  +  Kjr4  +  ...) 

+  [2F\x  y  +  P2  (r2  +  2y2)  ]  C  1  +  P3f2  +  P4  r4  +  ...  ] 


in  which  vx  and  vy  are  plate  measuring  residuals  and 


x  =  x°  -Xj, 

y  =  y°  -  y, 

r  =  +y3)^ 


observed  plate  coordinates  referred  to  principal  point. 
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When  these  relations  are  substituted  Into  (6),  the  projective  equations  explicitly  Involve  the 
coefficients  of  radial  and  decentering  distortion.  It  is  to  be  noted  that  in  formulating 
equation  (10)  v/e  implicitly  regard  the  optical  axis  and  principal  axis  as  being  coincident. 

This  formulation  can  introduce  no  error  of  metric  consequence  with  conventional  metric 
cameras  focussed  at  or  near  infinity,  for  here  the  optical  axis  and  principal  axis  must,  in 
fact,  nearly  coincide  ?f  acceptably  sharp  imagery  is  to  be  maintained  throughout  the  format. 

A  more  complicated  formulation  accounting  for  the  variation  of  distortion  with  object 
distance  would  be  needed  for  close-in  terrestrial  photogrammetry.  A  still  more  complicated 
formulation  would  be  needed  when  the  Scheimpflug  condition  is  exercised  in  terrestrial 
photogrammetry  to  maintain  perfect  focus  throughout  a  selected  plane  in  object  space  that 
is  not  perpendicular  to  the  optical  axis.  Although  the  present  treatment  will  be  limited  to 
cameras  focussed  at  or  near  infinity,  it  is  capable  of  being  extended  to  embrace  specialized 
applications  to  terrestrial  photogrammetry. 

We  shall  use  the  term  projective  parameters  to  denote  collectively  the  nine  elements 
of  orientation  (xp,  yp,  c,  or,  to,K,  XC,  Y®,  Z°)  and  the  coefficients  of  radial  and  decentering 
distortion  (K,,  K2,  . ..,  Pj ,  P2,  •••)•  We  shall  also  distinguish  between  'interior'  projective 
parameters  consisting  of  xp,  yp,  c,  Kj ,  K2,  •  •  • ,  P] ,  P2/  •  •  • ,  °nd  'exterior'  projective  param¬ 
eters  consisting  of  or,  to,  K,  XC,  Y0,  Z° .  Anticipating  later  results,  we  find  it  convenient  to 
introduce  the  symbol  uj  to  denote  the  i™  interior  projective  parameter  and  the  symbol  u’j  to 
denote  the  jt"  exterior  projective  parameter.  Specifically, 


i'¬ 

ll 

A 

Pi 

(12) 

«2=  yp 

Uj  =  K2 

u8  =  P2 

03=  c 

o$  =  K3 

p3 

and 

•• 

U|=  Of 

°x 

11 

(13) 

3 

ll 

• 

ii5  *  Y° 

•  • 

Uj=  K 

°N 

II 

:3* 

• 

In  02)  we  take  note  of  practical  experience  involving  the  calibration  of  a  wide  var.ety  of 
lenses  which  indicates  that  no  more  than  three  coefficients  are  likely  to  be  required  for  either 
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radial  or  decentering  distortion  (indeed,  with  but  few  exceptions,  two  coefficients  have  been 
found  to  suffice  for  decentering  distortion  and  one  or  tv/ o  coefficients  have  beer,  found 
to  suffice  for  radial  distortion). 

With  the  above  symbolism,  the  projective  equations  may  be  expressed  functionally 
as 


f,  (x°  +  vx,  y®  +  vy;  U| ,  u2, 
f2  (x®  +  vx,  y®  +  vy;  U| ,  u2. 


•  *  •  /  W?  U1  >  u2>  ’ '  ‘  •  ^6  )  =  0  # 
U9;  uj ,  u2,  v6)  -  0  . 


We  shall  now  proceed  as  if  ail  parameters  were  unknown  and  had  to  be  recovered 
from  plate  measurements  of  images  of  photographed  control  points.  Accordingly,  we  follow 
the  usual  piocedure  of  operating  on  observational  equations  that  have  been  linearized  about 
initial  approximations  by  means  of  Taylor's  series.  Thus,  in  (14)  we  replace  the  u's  and  u ‘s 
by  the  expressions 


U|  =  u®°  +  6  Ui  ,  i  =  1 , 2,  .  •  • ,  9  , 

(15) 

Uj  -  u®°  +  6uj  ,  J=1,2,  ...,6, 

in  which  the  u°®  and  u  00  denote  approximations  and  the  6  u  and  6u  are  associated  corrections 
When  equations  (15)  are  substituted  into  (14)  and  the  resulting  expressions  are  linearized,  we 
obtain  the  following  result 


(16) 


•*  +  •  * 
ajvx  +  °2Vy  +  b|6uj  +  b28u2  +  ...  +  bp80y 

+  Kj6U,  +  &2  8u’2  +  ...  +  b$S*u$  =  fj® 

a*]Vx  +  a2Vy  +  bj86j  +  b2  Su2  +  ...  +  b^8uy 

+  b\  8uj  +  b*2  6u2  +  ...  +  bj  SUj  =  f2°  . 


At  this  point,  we  introduce  the  subscript  i  to  distinguish  the  equatior>s  arising  from  the 
measurements  of  the  ifh  image.  The  pair  of  observational  equations  (16)  can  then  be 
repiesented  in  matrix  form  as 
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Explicit  expressions  for  the  elements  of  the  matrices  defined  in  (19)  are  provided 
in  Appendix  A.  There  it  is  shown  thut  the  coefficient  matrix  Aj  of  the  residual  vector 
becomes  equal  to  the  unit  matrix  when  negligible  first  order  terms  ore  dropped.  Accordingly, 
the  entire  set  of  linearized  observational  equations  arising  from  all  n  measured  control  points 
may  be  merged  into  the  single  matrix  equation 


(21)  v  +  B  6  +  B  *6  =  € 


in  which 


(22) 


v 

(2n,  1) 


V1 

v2 

,  B  = 

Bi 

b2 

•  • 

r—  —i 

Bl 

*1 

€2 

• 

m 

,  B  = 

• 

,  *  = 

• 

• 

vn 

(2n,9) 

• 

• 

Bn 

(2n,6) 

1 

• 

«  • 

Br* 

(2n, 1) 

• 

• 

e 

Kn 

We  shall  assume  that  the  covariance  matrix  A  of  the  observational  vector  is  of  The  form 

(23)  (2n/2n)  =  di°9  ( (2,2)  '  (2,2)  '  *  *  *  '  (2,V 


in  which 


(24) 


Ai 

(2,2) 


CTv 


*iyi 


xiyi 

V 


Thus,  A  is  c  block  diagonal  matrix  of  2  x  2  matrices  which  reed  not  necessarily  be  diagonal, 
a  fact  that  admits  consideration  of  plate  coordinates  that  have  been  derived  from  more  primary 
measurements  (e.g.,  measurements  made  on  a  polar  coordinate  comparator). 
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2.3.3.  Introduction  of  A  Priori  Constraints  on  Projective  Paramters 


Thus  far,  we  have  regarded  the  projective  parameters  as  being  completely  unknown. 

For  greater  generality  and  flexibility  we  shall  now  abandon  this  assumption  and  assume  that 
each  parameter  is  subject  fo  a  priori  constraints.  Accordingly,  we  introduce  the  'supplementary' 
sets  of  observational  equations 

Of  +  vj  =  uj  (5-1,2,...,  9) 

(25) 

Uj°  +  Vj  -  uj  (l  =  1 , 2, . . . ,  6) 

The  quantities  Uj  and  u'j  in  (25)  can  be  replaced  by  expressions  given  in  (15)  to  produce  the 
alternative  set  of  supplementary  observational  equations: 


$ 

-  5^  = 

fi  - 

0  =  1,2 . 9) 

(26! 

«  t 

V! 

(1 

CO 

ii 

- 

(1=  1,2,...,  6) 

in  which 

• 

•  •  #0* 

*  o 

(27) 

'l 

®  Ui  - 

ui  f 

*1 

M 

=  U,  - 

••  o 

uj  • 

These  may  be  represented  in  matrix  form  as 


• 

V  „ 

8 

• 

-  € 

» 

(9,1)  (9,1) 

(9,1) 

t 

(28)  .. 

V 

•  • 

8 

•  • 

(6,1)  (6,1) 

(6,1) 

t 

in  which 

Vl 

* 

*  • 

V1 

V 

*2 

^2 

*• 

v2 

<2 

s 

<• 

II 

♦ 

* 

• 

7  «  = 

• 

• 

•  « 

,  v  = 

• 

• 

•  ♦ 

7  C  = 

* 

* 

• 

e 

—  — 

* 

e 

•  • 

_v6_ 

of*  * 
1 _ 
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•  •  4 

We  shall  let  A  (a  9  x  9  matrix)  and  A  (a  6  x  6  matrix)  denote  the  specified  covariance 
matrices  of  the  observational  vectors  u°  and  u  °,  respectively.  Later,  we  shall  discuss 
various  applications  of  these  covariance  matrices. 


2.3.4  Formation  of  Normal  Equations 

The  observational  equations  defined  by  (21)  and  (28)  can  be  merged  into  the 
composite  equation 


(30) 


V 

• 

B 

••*1 

B 

« 

8 

c 

4 

V 

+ 

-1 

0 

8* 

= 

t 

•  • 

V 

.0 

-I 

L  -j 

v  • 

€ 

which,  with  obvious  notation,  can  be  expressed  compactly  as 


(31)  v  +  B6  =  *•  . 

The  composite  covariance  matrix  associated  with  v  is  of  the  form 


(32)  A 

(2n  +  15,  2n  + 1 5) 


A  0 

(2n,2n) 


0 

0 


A 

(9,9) 


0 

0 

•  * 

A 

(6,6). 


The  normal  equations  arising  from  the  least  squares  adjustment  of  the  observations 
governed  by  (31)  and  (32)  are  of  the  form 


IN  6  C 

'  '  (15,15)05,1)  (15,1)  ' 


in  which 

N  =  B  5  A  B  , 
(34)  r  -1 

c  =  B  A  *  . 
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By  virtue  of  the  structure  of  B,  5,  A  and  f  ,  the  normal  equations  can  be  expressed  in  the 
more  expanded  form 


(35) 


|V|  +  W  N 

• 

6 

•  w 

i 

«o 

i _ 

FiT  n  +  W. 

LsJ 

_c  -  'A  vj 

in  which 


Ki  =  BTWB  ,  c  =  BT  W  f  , 

(36)  N  =  BT  W  B  ,  c  *  BT  W  €  , 

•  •  T  •  • 

n  =  r  WB 


wherein 

(37)  W  =  A"1  ,  W  =  A"1 ,  W  =  K~' 


By  virtue  of  the  block  diagonality  of  A ,  the  expressions  for  N,  N,  Kl,  c,  c  can  be  expressed 
as  sums 


N  =  £  N,  , 


N  =  £  N  , 

i=i  1 


in  which 


c  = 


(38)  N  =  £  N,  ,  c  = 
i=l 


i=l  *  ' 

I,  - 


w, 

K 

,  c  =  §T 

'  i  wi 

W,  , 

N,  -  6,T 

wt 

\ 

,  e  -  BT 

9  i  i 

W,  , 

**T 

<t  • 

Bt 

f  =  B 

1  i 

w, 

9 
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where 


(40)  W,  =  A"1 


It  should  be  appreciated  that  the  above  development  of  the  normal  equations  is 
strictly  derivational.  The  oomputatu  nai  flow  itself  is  actually  quite  short  and  consists  of 
the  following  major  steps: 

(1)  By  means  of  formulas  given  in  Appendix  A,  the  basic  matrices  Bj  and  ere 
evalued  from  the  following  given  data; 

(a)  The  observed  plate  coordinates  Xj0,  y°  ; 

(b)  The  direction  cosines  ,  pt  ,V  (in  the  case  of  stellcr  control)  or  the  Xt , 

Yt ,  Zj  coordinates  of  the  Jf-h  control  point; 

(c)  ~he  approximations  u^00  of  the  projective  parameters. 

(2)  In  terms  of  the  computed  Bt,  Bt ,  together  with  the  given  A  ,  eqs.  (39)  are 

evaluated  to  generate  ct ,  c  . 

(3)  As  each  ,  N.  ,  Nj ,  cl ,  c  is  formed,  it  is  added  to  the  sum  of  its  predecessors 
thereby  ultimately  generating  the  N,  N,  53,  c,  c  indicated  in  equation  (38). 

(4)  From  the  given  a  priori  values  u°,  u  °  and  the  associated  covariance  matrices  A 

and  A,  the  terms  W,  W,  W,  c  and  W  *  are  evaluated  in  accordance  with  equa¬ 
tions  (27),  (29),  and  (37). 


(5)  From  the  results  of  steps  (3)  and  (4),  the  normal  equations  (35)  are  formed. 


2.3.5  Computation  of  Residuals 

Once  the  normal  equations  have  been  formed,  they  may  be  solved  for  the  vectors  6,  8 
which  provide  corrections  to  the  initial  approximations.  The  improved  approximations  would 
then  be  used  to  initiate  an  iterative  process  which  would  continue  until  the  revised  corrections 
become  negligibly  small.  Upon  convergence  (which  is  normally  obtained  in  three  iterations), 
the  plate  residuals  for  the  i*h  point  may  be  computed  from 
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(41) 


t 


Vj  =  Cj 

(2,1)  (2,1) 

in  which  the  discrepancy  vector  ft  is  based  on  the  final  set  of  projective  parameters.  Similarly, 
the  residuals  of  the  a  priori  values  of  the  projective  parameters  may  be  computed  from 

•  *  ~  t 

(42) 

V  =  i 

•  «  • 

in  which  the  elements  of  i  and  *  are  generated  from  equations  (27)  with  *he  final  set  of 
projective  parameters  being  used  for  Uj00  and  uj30.  The  quadratic  form  of  the  composite  residual 
vector  is 

(43)  q  =  J  YVy  +  ‘v^  Wv  £  (vj  Wj  Vj )  . 

The  statistical  degrees  of  freedom  associated  with  q  may  be  taken  as  2n-k,  where  k  denotes 
the  number  of  projective  parameters  for  which  no  effective  a  priori  constraints  are  exercised. 


2.3.6  Practical  Exercise  of  A  Priori  Constraints 


Not  only  do  the  covariance  matrices  A  and  A  provide  proper  weighting  of  actual  a  priori 
values  of  projective  parameters,  but  they  also  can  be  exploited  for  other  purposes.  With  regard 
to  any  given  projective  parameter,  three  situations  arise; 

(1)  The  parameter  is  actually  known  in  advance  to  a  worthwhile  degree  of  accuracy; 
here,  the  pertinent  elements  covariance  matrix  would  realistically  reflect  ones 
knowledge  of  the  parameter. 

(2)  The  parameter  is  not  known  in  advance  to  a  worthwhile  degree  of  accuracy;  here, 
the  pertinent  variance  would  be  deliberately  inflated  by  several  orders  of  magni¬ 
tude  over  what  is  considered  io  be  a  reasonable  measure  of  the  uncertainty  of  the 
adopted  initial  approximation.  This  allows  the  parameter  virtually  unrestricted 
freedom  to  adjust. 
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(3)  The  parameter  is  not  to  be  exercised  in  the  adjustment  (most  commonly# 

the  case  with  higher  order  coefficients  of  radial  and  decentering  distortion); 
here  the  pertinent  variance  would  be  set  equal  to  a  value  effectively  equal 
to  zero.  An  appropriate  value  would  be  one  several  orders  of  magnitude 
smaller  than  the  variance  to  be  expected  for  the  recovery  of  the  parameter 
if  it  were  free  to  adjust. 

By  thus  using  the  covariance  matrix  of  the  projective  parameters  to  exercise  general  control 
over  the  adjustment#  one  not  only  avoids  the  need  for  programming  a  variety  of  special  options 
but  one  also  gains  added  simplicity  by  avoiding  the  redimensioning  of  operational  matrices. 

A  specific  application  of  the  process  of  exercising  a  priori  constraints  to  effect 
parametric  control  involves  the  solution  for  decentering  distortion.  When  the  higher  order 
coefficient  P3  is  equal  to  zero#  the  expressions  in  the  projective  equations  accounting  for 
decentering  distortion  become  linear  in  the  leading  coefficients  Pt  and  Pa.  This  means  that 
one  can  adopt  values  of  zero  as  initial  approximations  for  Pj  and  P8.  However#  by  doing 
this#  one  causes  the  entire  set  of  coefficients  corresponding  to  P3  in  the  B  matrix  to  be  zero# 
thereby  rendering  the  normal  equations  indeterminate.  The  remedy  consists  of  using  a  priori 
constraints  to  suppress  P3  to  zero  in  the  initial  iterations.  Once  stable#  nonzero  approximations 
have  been  obtained  for  PT  #  P8#  the  constraints  on  P3  cap  be  relaxed  in  subsequent  interations 
to  permit  recovery  of  this  parameter  together  with  final  refinement  of  Pj  and  P3. 

Another  application  of  the  special  exercise  of  a  priori  constraints  is  of  particular 
value  in  the  calibration  of  cameras  having  narrow  angular  fields.  Here#  the  projective  effects 
of  a  small  translation  6xp,  flyp  of  the  plate  are  very  nearly  equivalent  to  the  projective  effects 
effects  of  a  small,  suitably  directed  rotation  of  the  camera.  This  coupling  of  translation  and 
rotation  can  sometimes  lead  to  a  sufficiently  ill-conditioned  system  of  normal  equations  to 
prevent  convergence  of  the  iterative  process.  The  problem  is  compounded  when  decentering 
distortion  is  being  recovered,  for  decc  ntering  coefficients  also  interact  to  a  moderate  degree 
with  xp,  yp  .  The  remedy  consists  of  subjecting  xp,  yp  in  the  initial  reduction  to  fairly  tight 
a  priori  constraints  (e.g.,  a  few  hundred  microns)  and  of  relaxing  such  constraints  by  stages 
in  subsequent  iterations  until  they  ultimately  become  inconsequential. 
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2.3.7  Error  Propagation 

The  covariance  matrix  of  the  adjusted  projective  parameters  is  provided  by  the 
inverse  of  N,  the  coefficient  matrix  of  the  normal  equations.  If  n11  denotes  the  element 
in  the  |th  row  and  column  of  N"*,  it  follows  from  the  ordering  of  the  parameters  that  the 
standard  deviations  of  the  adjusted  elements  of  interior  orientation  are  given  by 

®xp  =  (n")1 

(44)  ayp  =  (n»)J 

vc  =  ("”)* 

The  variances  of  the  radial  and  decentering  parameters  are  provided  by  elements  n44  through 
ri».  However,  these  are  of  little  interest  in  themselves.  Rather,  what  is  of  interest  is  the 
uncertainty  of  the  calibrated  distortion  functions  throughout  4he  format. 

In  the  case  of  symmetric  radial  distortion,  the  propagation  of  errors  into  the  distortion 
function  is  readily  accomplished.  Inasmuch  as  the  covariance  matrix  of  the  coefficients  of 
radial  distortion  is 


$ 

1  c 

n« 

n44' 

(45)  ak  = 

rs45 

n55 

n* 

n44 

n« 

n« 

and  the  distortion  function  is  defined  by 

(46)  6r  =  Kji3  +Kar6 +Kgr7  +  ..., 

it  follows  from  the  theory  of  error  propagation  that  the  standard  deviation  of  Sr  for  an  arbitrary 
radial  distance  is  given  by 

(47)  a6r  =  (uknku') 


-1?- 


5n  which 


(48)  uk  -  =  03 rS 

By  means  of  this  result  one  con  generate  confidence  intervals  (or  error  bounds)  associated  with 
the  calibrated  distortion  curve. 

In  the  case  of  decentering  distortion  the  error  propagation  is  somewhat  more  involved 
since  the  end  results  are  most  conveniently  expressed  in  terms  of  the  tangential  profile  function 
Pr  and  its  associated  phase  angle  0O.  By  virtue  of  (5)  we  have 

J,  =  (P,2  +  P22)^ 

(49)  jj  =-<p?  +  fir 

0O  “  -arc  tan  (Pj/P2) 

Inasmuch  as  the  covariance  matrix  of  the  adjusted  decentering  parameters  Pj ,  P2,  P3  is 
given  by 


n77  n78  n77 


(50)  Qp  =  n78  n88  r.87 

n77  n»  n77 


It  follows  that  the  covariance  matrix  of  the  derived  parameters  Jj ,  J2,  0o  is  given  by 


(5!)  Cc  -  (up  f<pUa) 


;n  which 


(52)  u,  °  °  Iff  +  Pj2)"*  -P|Pj  -Pi2  +  ft)  ■ 


[_cos  0O  s*n  0e  0 


inasmuch  as  the  ptofile  function  is  defined  by 


(53)  Pr  =  Ji  r2  +  Jjr4  +  .  .  .  , 

the  covariance  matrix  of  the  profile  function  and  the  phase  angle  is  given  by 


(54) 


or, 


Pr 


Pr  0« 


CT, 


Pr0e 


4)  fib  Mo  ' 


in  which 


(55)  U0 


-  K)  _  fr2  r4 

*0]7 [0  o 


0 

1 


By  virtue  of  this  result,  error  bounds  can  be  established  for  the  representation  of  decenterina 
distortion  in  terms  of  tangential  profile  and  phase  angle. 

In  the  error  propagations  for  both  radial  and  decentering  distortion  we  ignored  the 
fact  that  the  coordinates  xp,  yp  of  the  principal  point  are  implicit  in  the  radial  distances  r 
(since  the  radial  distances  are  referred  to  the  calibrated  principal  point).  It  turns  out,  how¬ 
ever,  tha'-  the  errors  in  the  xp,  yP  uave  only  a  second  order  effect  on  the  above  error  propa¬ 
gations  and  thus  need  not  be  considered. 

The  radial  distortion  function  6r  generated  by  the  adjustment  corresponds  to  the  cali¬ 
brated  principal  distance  c.  As  is  shown  in  Brown  (1956)  this  function  con  be  transformed  into 
a  projectively  equivalent  distortion  function  Sf  associated  with  an  arbitrery  specified  principal 
distance  c+Ac.  The  transformation  is  expressed  by 

(56)  s;  =  (1+~)5r  +  fr  , 


which  leads  to  the  alternative  version  of  the  distortion  function 


(57)  s;  =.  KJr  +  KJr3  +  +  K;  ^  +  ... 

in  which 

(58)  K*  Ac/c,  Ki  ^(l+Ac/c)^,  K’  =  (1 +Ac/c)  Ks,  etc. 


One  can  use  this  result  to  force  the  transformed  distortion  function  to  assume  the 
value  of  zero  at  a  specified  radial  distance  r0.  Thus,  if  6ro  denotes  the  value  of  8r  at  r  =  r0, 
it  follows  from  (56)  that  the  choice 


(59)  Ac  =  -(-Ts-)Sr 
ro  °r. 


renders  8^  equal  to  zero  at  r  =  rc. 

Alternatively,  Ac  can  be  chosen  so  that  the  mean  value  of  the  distortion  function 
out  to  a  specified  radial  distance  r0  is  zero.  This  results  in  the  choice 


(60)  Ac 

?n  which 
(51)  w 


wc 

w  +  I  r\ 


_  ri  ,  K, 

~"2~  VT" 


)  • 


Still  another  choice  consists  of  transforming  the  distortion  function  to  minimize  the 
rms  value  of  distortion  out  to  a  specified  radial  distance  r0.  This  is  accomplished  by  choosing 
Ac  cs 


(62)  Ac  =  - 


(s  ±w)  c 
2$  +  w  +  ri 


in  which  w  is  as  above  and 
(63)  s 


-  J  r£r  dr  -  rj  (~  +  -^  r| 


K,  . 

+  -f  r*  + 


•  ) 


The  final  transforma  Mon  to  be  considered  is  one  in  which  the  largest  positive  and 
negative  values  of  distortion  out  to  specified  radial  distance  r0  are  forced  to  have  the  same 
absolute  value.  Here,  the  appropriate  Ac  is  most  readily  established  by  a  two  stage  process. 

The  first  stage  consists  of  choosing  an  initial  value  Ac0  in  accordance  with  either  (60)  or  (62). 
One  then  determines  the  slope  m  of  the  straight  line  that  passes  through  the  origin  and  is  such 
that  the  largest  positive  and  negative  departures  of  the  transformed  distortion  function  from  the 
line  are  of  equal  magnitude.  This  determination  can  be  made  either  graphically  (with  the  aid 
of  a  transparent  ruler)  or  numerically  (by  a  trial  and  error  process).  In  either  case,  the  final 
value  of  Ac  is  given  by 

(63)  Ac  =  Aco-m(c  +  Ac0)  . 

All  of  the  above  transformations  yield  projectiveiy  equivalent  results.  When  Ac/c  is 
small  (as  is  generally  the  case),  the  error  in  the  calibrated  value  of  c  has  no  sensible  effect  on 
the  transformation  and  may  be  safely  ignored.  Also  when  Ac/c  is  small,  the  standard  deviation 
given  by  (47)  may  be  used  in  conjunction  with  the  transformed  distortion  function. 

A  decided  advantage  of  the  analytical  approach  to  camera  calibration  is  that  the  error 
propagation  associated  with  the  end  results  can  be  readily  computed.  By  contrast,  conventional 
laboratory  methods  do  not  lend  themselves  to  precise  statements  concerning  the  accuracies  of  the 
results. 


2.4  Simultaneous  Multiframe  Analytical  Calibration  (SMAC). 

2.4.1  Introduction 

The  results  of  the  preceding  section  are  not  new.  Rather,  they  constitute  a  detailed 
review  of  our  earlier  theory.  The  development  has  been  specifically  constructed  to  facilitate 
the  derivation  of  the  SMAC  calibration.  Before  proceeding,  we  would  emphasize  that  the 
theory  applies  both  to  a  stellar  calibration  and  to  an  aerial  calibration  (i . e . ,  a  calibration 
using  aerial  photographs  taken  over  a  suitable  testing  range).  In  the  stellar  calibration,  the 
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coordinates  X  ,  Y  ,  Z  of  the  came, a  are  inherently  unrecoverable  (and  hence  are  suppressed) 
inasmuch  as  all  control  points  are  sensibly  ot  infinity.  In  the  aerial  calibration,  they  are  carried 
as  unknowns  subject  possibly  to  a  priori  constraints  (as  when  the  precise  position  of  the  aircraft 
is  established  by  a  tra  iking  system).  While  the  elements  of  interior  orientation  (xp,  yp/c)  con 
be  recovered  in  a  stellar  calibration,  they  may  not  be  recoverable  in  ar.  aerial  calibration. 

This  stems  from  the  fact  that  when  all  control  points  are  nearly  coplanar,  a  small  translation 
(6X  ,  6Y  ,  5Z  )  in  object  space  is  projectiveiy  equivalent  to  a  small  translation  (5xp  ,8yp  ,5c) 
in  image  space.  It  follows  that  the  recovery  of  xp,  yp,  c  in  an  aerial  calibration  requires  either 
a  strongly  three  dimensional  distribution  of  control  or  else  precise  external  hacking  of  the  air¬ 
craft.  This  matter  will  be  discussed  more  fully  later.  It  should  be  kept  in  mind  throughout  the 
development  to  follow  that  in  extending  the  single  frame  calibration  to  the  SMAG  calibration, 
we  shall  preserve  the  dual  applicability  of  the  reduction  to  stellar  and  aerial  photographs. 


2.4.2  Formation  of  the  General  Normal  Equations 


We  now  consider  the  genera!  situation  in  which  an  indefinitely  large  number  of  frames 
of  photographed  control  are  to  be  employed  for  camera  calibration.  We  assume  that  a  common 
set  of  interior  projective  parameters  applies  to  ail  frames  and  that  (possibly)  a  different  set  of 
exterior  projective  parameters  applies  to  each  frame.  To  begin,  we  shall  assume  that  only 
three  frames  are  involved  in  the  reduction.  Then  ignoring  momentarily  the  a  priori  constraints 
on  the  interior  projective  parameters,  we  obtain  the  following  three  sets  of  normal  equations  for 
the  independent  single  frame  reductions: 
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If  the  only  avaiioble  information  were  to  consist  of  the  a  priori  constraints  on  the  interior  projective 
parameters,  the  following  system  of  normal  equations  would  apply 

(65)  W  S  =  -Wc  . 

Now  it  has  been  shown  (Brown,  1957)  that  systems  of  normal  equations  formed  from  consistently 
weighted,  independent  observations  and  operating  on  a  common  parametric  vector  are  additive: 
that  is,  they  can  be  summed  to  form  the  system  of  normal  equations  that  would  have  arisen  from 
the  adjustment  of  the  combined  observational  vectors-  This  additive  property  of  independent 
systems  of  normal  equations  cannot  be  applied  directly  to  the  normal  equations  written  above 
because  each  operates  on  a  different  vector.  This  difficulty  is  easily  overcome  by  a  process  of 
zero  augmentation  to  force  each  system  to  operate  on  the  total  parametric  vector  to  be  recovered. 


Accordingly,  equations  (64)  and  (65)  may  be  rewritten  as 
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(66c) 
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Thus  augmented,  each  system  of  normal  equations  now  operates  on  the  same  parametric  vector, 
inasmuch  as  each  is  formed  from  an  independent  set  of  observations,  we  are  thus  in  a  position 
to  apply  the  additive  property  of  independent  sets  of  norma!  equations  to  obtain  the  following 


result: 
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in  which  N  and  c  are  defined  by 


N  =  N,  +  N2  +  N3  , 

C  =  C{  '•  c2  +  c3  . 


This  system  represents  the  normal  equations  resulting  from  a  calibration  bused  on  the  simultaneous 
adjustment  of  plate  coordinates  from  three  arbitrary  frames. 


The  extension  of  the  above  development  to  an  arbitrarily  large  number  of  frames  m  is 
immediate.  The  resulting  system  of  normal  equations  is 
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where  now 

N  -  N|  +  N2  +  *  •  +  Nm 

(71) 

C  -  Cj  +  c2  +  .  •  .  4  cm 

The  above  derivation  of  the  normal  equations  is  designed  to  emphasize  the  close  relationship  that 
exists  between  the  basic  single  frame  reduction  and  the  composite  multiframe  reduction.  An 
alternative  derivation  closely  paralleling  the  procedure  developed  in  Brown  (1958)  is  given  in 
Gambino  (1967). 
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2.4.3  Solution  of  the  Normal  Equations 

The  above  system  of  normal  equations  is  of  order  6m  •!  9  (for  the  general  case)  and 
thus  increases  linearly  with  the  number  of  frames  m.  When  m  is  small,  one  con  solve  the 
system  directly  without  difficulty.  Our  concern,  however,  is  with  the  general  case  in  which 
m  is  without  set  limit.  Here,  a  conventional  direct  solution  becomes  Jmpracticai  after  a  certain 
point  inasmuch  as  the  number  of  computations  is  proportional  to  (6  m  +  9)3 .  Fortunately,  the 
system  of  normal  equations  (68)  has  the  same  general  structure  as  the  system  that  we  had  developed 
in  Brown  1958  for  the  rigorous  adjustment  of  a  general  photogrammetric  net.  There  we  shewed  that 
the  partial  block  diagonaiity  of  the  coefficient  matrix  could  be  exploited  tc  produce  an  efficient 
solution  in  which  compulations  increased  only  linearly  with  the  number  of  unknown  6  vectors. 
Inasmuch  os  the  derivation  of  the  solution  is  given  in  the  above  reference  and  is  further  elaborated 
in  Brown,  Davis,  Johnson  1964,  we  shall  not  repeat  the  derivation  here.  Instead,  we  shah  merely 
adapt  results  pertinent  to  the  problem  at  hand. 


Although  the  general  normal  equations  for  t1  e  SMAC  calibration  are  given  by  (70),  one 
would  not  actually  set  up  this  system  in  the  course  of  the  reduction.  Rather,  one  would  generate 
the  reduced  system  of  normal  equations  resulting  from  the  elimination  of  the  6^'s  from  the  general 
system.  This  entails  the  following  steps.  Starting  with  the  first  frame  (J  =  1),  one  generates  the 
coefficient  matrix  and  constant  column  for  the  single  frame  calibration  according  to  the  develop- 

•  MM*  •  •  •  •  , 

ment  of  the  preceding  section.  This  provides  the  primary  matrices:  N^,  N^,  ,  c^  and 

•  • 

cr 

(72) 


Wj  i'  .  In  terms  of  these,  one  computes  the  following  set  of  auxiliary  matrices 
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As  Sj  and  are  formed  for  Hie  frame,  they  are  added  to  the  sums  of  their  predecessors. 
After  all  m  frames  have  thus  been  processed,  this  produces  the  quantities: 

S  ~  S,  +  S2  +  ...  +  Sn  , 

(76) 

c  ~  c,  +  c2  +  ...  +  cn  . 

The  reduced  system  of  normal  equations,  a  9  x  9  system  involving  only  the  interior  projective 
parameters,  is  then 

(77)  (S  +W)$  =  c  -  We  • 

The  solution  is  given  by 

(78)  6  =  M(c-V/c) 
in  which 

(79)  M  =  (S  +  W)->  . 


Once  the  vector  of  interior  projective  parameters  has  thus  been  obtained,  each  vector  of  exterior 
projective  parameters  can  be  computed  in  turn  from 


(80)  8j  =  (NJ  +  Wj)-'  (c}  -  c})  -  0,6,  (1  =  1,  Z,  ...,  m). 


The  values  obtained  from  (78)  and  (80)  are  added  to  the  original  approximations  to  obtain  improved 
approximat'  'ns.  These  are  then  used  to  initiate  cn  iterative  process  that  is  cycled  to  convergence. 
Computafic  of  observational  residuals  follows  the  procedure  of  Section  2.3-5. 


2.4.4  hrror  Propagation 

Thu  matrix  M  represents  the  covariance  matrix  of  the  interior  projective  parameters. 
Accordingly,  witi.  suitable  reinterpretation,  the  resuits  of  Section  2.3-7  for  error  propagation 
associated  with  a  single  frame  calibration  can  be  applied  to  the  ShAAC  calibiation. 
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The  covariance  matrix  of  the  adjusted  values  of  the  exterior  projective  parameters 
for  the  J^1  frame  is  given  by 

(sn  Mj  =  (Kl,  +  Wj)-1  +  Qj  Mj  Q,?. 

(6,6)  (6,6)  (6,9)  (9,9)  (9,6) 

The  first  term  on  the  right  represents  the  covaiiance  matrix  for  the  limiting  case  in  which  the 
interior  projective  parameters  are  perfectly  known,  and  the  second  term  accounts  for  the 
dilution  of  limiting  accuracies  attributable  to  errors  in  the  adjusted  values  of  the  interior  pro¬ 
jective  parameter.  The  contribution  of  the  second  term  automatically  becomes  suppressed  into 
relative  insignificance  when  the  SMAC  reduction  involves  a  moderate  number  of  frames  each 
containing  a  moderate  number  of  control  points. 

2.4.5  General  Remarks 

The  interior  projective  parameters  resulting  from  a  SMAC  calibration  may  be  considered 
as  the  best  possible  compromise  for  the  se^  of  frames  carried  in  the  reduction.  This  means  that 
errors,  both  random  and  systematic,  arising  from  a  variety  of  sources  to  be  mentioned  below  are 
in  large  measure  averaged  out  insofar  as  their  effects  on  the  interior  projective  parameters  are 
concerned.  In  considering  the  results  of  a  SMAC  calibration,  it  is  convenient  to  classify  systematic 
errors  into  the  following  two  categories: 

0)  transient  systematic  errors, 

(2)  persistent  systematic  errors. 

Transient  systematic  errors  are  those  errors  that  are  systematic  on  a  given  frame  but  are  essentially 
independent  from  frame  to  frame.  Persistent  systematic  errois  arc  those  that  tend  to  remain 
highly  correlated  from  frame  to  frame. 

In  an  aerial  SMAC  calibration  one  finds  it  easy  to  postulate  a  large  number  of  possible 
sources  of  transient  systematic  enor  (particularly,  if  quantitative  estimates  are  not  demanded). 
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A  list  might  include: 

(a)  variations  in  thickness  of  film-emulsion  combination; 

(b)  nonunHorm  dimensional  instability  of  film; 

(c)  broad  refraciive  anomalies  induced  by  atmospheric  turbulence; 

(d)  random  local  failure  of  film  to  conform  precisely  to  platen; 

(e)  angular  motion  of  camera  during  exposure; 

(f)  variations  in  orientation  of  window  and  shockwave  relative  to  camera  because 
of  variations  in  aircraft  attitude; 

(g)  dynamic  deformations  of  camera  induced  by  aircraft  accelerations  or  by  corrective 
accelerations  by  camera  mount; 

(h)  deformations  attributable  to  thermal  imbalances; 

(i)  time  varying  component  of  systematic  errors  in  compaiafor  used  to  measure  film; 

(j)  variations  in  personal  equation  of  measuring  personnel  during  course  of  measuring 
(e.g.,  fatigue  factor). 

No  doubt,  other  potential  contributors  to  transient  systematic  error  could  be  added  to  t'ne  list. 

The  point  is  that  while  the  total  effect  of  such  errors  may  assume  significance  on  a  given  frame, 
the  effect  will  tend  to  be  independent  from  frame  to  frame.  Because  the  combination  of  a  moderate 
number  of  individually  inconsequential,  second  order  systematic  effects  can  amount  to  a  definitely 
significant  first  order  effect,  there  is  little  merit,  in  our  view,  in  employing  in  an  aerial  SMAC 
calibration  more  than  50  or  so  well-distributed  images  on  any  single  frame.  In  many  instances, 
the  reduction  of  500  images  on  a  single  frame  would  produce  only  a  slight  actual  improvement 
over  the  results  obtainable  from  the  reduction  of  50  images.  This  is  because  errors  in  plate  co¬ 
ordinates  ore  not  strictly  independent,  as  one  would  like  to  believe,  but  rather  are  correlated  by 
virtue  of  residual  systematic  effects.  Accordingly,  it  follows  that  redundancy  is  effective  only 
up  to  a  certain  point  on  a  single  frame  and  that  further  improvement  from  redundancy  must  be 
wrought  from  exploitation  of  random  frame-to-frame  vaiiations  of  otherwise  systematic  effects. 
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For  this  rc-oson,  a  SMAC  reduction  of,  say,  20  frames  having  on  average  of  25  control  points 
per  frame  is  vastly  to  be  preferred  over  a  reduction  involving  a  total  of  500  points  on  a  single 
frame. 

Unlike-  the  effect?  of  transient  systematic  errors,  the  effects  of  persistent  systematic 
errors  do  not  average  out  in  a  SiWkC  reduction.  Sources  of  such  error  ir  an  aerial  calibration 
include: 

(a)  departure  of  the  platen  from  a  true  plane; 

(b)  tilt  of  the  platen  at  instant  of  exposure  (pat  tied iy  attributable  to  (a)  ); 

(c)  prism  effect  of  aircraft  window; 

(d)  curvature  of  aircraft  window  induced  by  pressure  and  temperature  differentials 
between  inside  and  outside  of  aircraft; 

(e)  refractive  effects  of  shock  waves  of  aircraft; 

(f)  persistent,  uncompensated  component  of  film  deformation; 

(g'  image  aberrations,  other  than  true  optical  distortion,  that  affect  measured  coordinates 
(e.g.,  in  photographic  measurements,  coma  superimposes  on  classical  radial  distortion 
a  component  that  is  dependent  both  on  radial  distance  and  on  mean  image  diameter); 

(h)  systematic  errors  in  reseau  coordinates  employed  to  compensate  for  film  deformation; 

(i)  systematic  deficiencies  in  image  motion  compensation  (IMC)  system; 

(j)  errors  in  the  given  survey  of  ground  control; 

(k)  uncompensated  systematic  errors  of  comparator  used  to  measure  plates- 

Persistent  systematic  errors  may  be  partially  absorbed  by  the  calibrated  projective  parameters. 

Thus  a  platen  that  happens  to  be  essentially  spherically  concave  or  convex  would  infioduce  an 
effect  that  would  be  completely  absorbed  by  ihe  calibrated  coefficients  of  radial  distortion. 

A  similar  remark  applies  to  the  effeci  of  the  spherical  component  of  the  aircraft  window.  For 
flights  at  a  given  altitude  and  velocity  in  a  given  aircraft,  the  various  layers  bounded  by  shock 
waves  act  to  a  first  approximation  as  very  weak  prisms  and  hence  produce  effects  that  can 
largely  be  absorbed  by  the  calibrated  parameters  of  deccntering  distortion.  A  similar  remark 
applies  to  the  prism  effect  cf  the  aircraft  window.  A  slight  tilt  in  the  camera  platen  at  the 
instant  of  exposure  would  be  projective!/  compensated  by  c  shift  in  the  coordinates  Xp,  >'p 
of  the  principal  point.  In  o  camera  exercising  imceo  motion  compensation,  such  a  shift  would 
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also  compensQie  for  systematic  failure  of  the  fiducials  to  flash  at  the  precise  center  of  exposure 
(for  a  given  V/H  ratio). 

That  portion  of  persistent  systematic  errors  that  is  not  accommodated  by  projective  paiam- 
efers  would  be  reflected  (along  with  random  errors  and  transient  systematic  errors)  in  the  final 
set  of  least  squares  residuals  resulting  from  the  SMAC  calibration.  The  residual  vectors  from  all 
frames  plotted  against  x,y  plate  coordinates  on  a  common  graph  can  be  examined  for  possible 
uncompensated  systematic  effects.  From  numerical  and  statistical  analysis  of  such  residuals,  one 
can  derive  either  error  contour  maps  or  functional  representations  of  residual  systematic  errors. 

These  can  subsequently  be  used,  if  due  precautions  are  taken,  to  supplement  the  calibrated  exterior 
projective  parameters  in  the  correction  of  operational  data.  A  practical  example  of  this  approach 
will  be  given  later. 

From  the  foregoing  considerations,  St  is  clear  that  a  SMAC  calibration  constitutes,  not 
merely  a  lens  calibration,  but  a  total  metric  calibration  of  the  entire  photogrammetric  system 
under  operational  conditions.  To  the  extent  that  all  pertinent  elements  of  the  system  approximate 
in  routine  operations  those  applying  to  the  photography  used  in  the  SMAC  calibration,  the  compen¬ 
sation  for  persistent  systematic  error  afforded  by  the  calibrated  parameters  together  with  error  functions 
derived  from  residuals  is  altogether  desirable.  On  the  other  hand,  such  compensation  would  be 
partially  or  totally  unsound  in  instances  where  key  contributors  to  systematic  error  have  character¬ 
istics  differing  significantly  from  those  prevailing  in  the  SMAC  photography.  For  this  reason, 
attempts  should  be  made  to  isolate,  insofar  as  possible,  the  contributions  of  those  sources  of 
systematic  error  that  are  likely  to  vary  from  one  operation  tc  the  next.  This  matter  will  be  discussed 
more  fully  later. 

3-  APPLICATIONS 

3-1  Introduction 

Having  developed  lire  theorciical  basis  far  SMAC,  we  shall  now  direct  our  attention  toe 
number  of  specific  applications  of  the  method.  These  have  been  selected  to  illustrate  the  power, 
flexibility  ond  universality  of  the  approach.  Most  of  tire  applications  are  illustrated  wiih  actual 
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results,  some  of  which  may  be  of  considerable  significance  to  the  future  of  the  art  of  camera 
calibration.  In  certain  instances,  comparisons  are  made  between  SMAC  results  and  results 
produced  by  conventional  laboratory  methods.  Comparisons  are  also  made  between  stellar 
SMAC  and  aerial  SMAC  calibrations  of  the  same  camera.  The  application  of  SMAC  to 
laboratory  calibrations  is  discussed  and  is  demonstrated  by  SMAC  calibrations  of  Lunar  Orbiter 
cameras.  Consideration  is  also  given  to  applications  of  SMAC  to  the  calibration  of  ballistic 
cameras  and  to  the  analysis  of  their  physical  stability. 

3.2  Calibration  of  KC--6A  Aerial  Cameras 

3*2. 1  Design  and  Execution  of  Aircraft  Tests 

The  KC-6A  Camera  is  the  most  recent  and  most  advanced  of  the  aerial  mapping  cameras 
developed  by  the  U.  S.  Air  Force.  It  employs  the  150mm  f/5*0  Geocon  IV  lens  designed  by 
Dr.  James  Baker.  The  lens  has  a  resolution  of  45  lines/miiiimeter  AWAR  on  Plus  X  film  and  a 
resolution  of  not  less  than  25  lines/miiiimeter  anywhere  within  a  9x9  inch  format.  The  camera 
incorporates  image  motion  compensation  (IMC),  automatic  exposure  control  and  a  platen  reseau 
that  contains  a  set  of  25  points  evenly  spaced  at  2  inch  intervals  and  a  set  of  12  edge  points 
associated  with  the  fixed  fiducials.  It  also  incorporates  an  optica!  linkage  to  the  vertical  gyro 
of  the  Hypernas  li  Inertial  Navigator.  This  permits  rhe  precise  determination  of  the  direction  cf 
the  earnerc  axis  with  respect  to  the  local  vertical.  Further  details  concerning  the  KC-6A  Camera 
are  given  by  Livingstone  (1966). 

A  major  motivation  for  the  development  of  SMAC  wa.  its  application  to  the  flight  testing 
of  the  KC-&A  os  part  of  the  USQ-28  geodetic  mapping  and  surveying  system.  In  view  of  the 
extensive  development  and  demanding  requirements  of  the  KC-6A,  the  decision  was  made  by  the 
U.  S-  A-my  and  Air  Force  to  implement  a  tesi  designed  by  DBA  to  perform  a  definitive  calibra¬ 
tion  of  the  carneio  in  its  operational  environment.  A  desired  by-product  of  the  calibration  was 
to  consist  of  the  delerminoiion  of  the  precise  angular  elements  of  orientation  of  a  series  of 
exposures  for  the  purpose  of  evaluating  the  accuiacies  of  the  vertically  system. 
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The  testing  program  was  based  or.  flights  over  the  McLure  Phofogrammetric  Test  Range, 
a  targeted  range  established  in  Ohio  by  the  US  Coast  and  Geodetic  Survey.  Inasmuch  as  one 
of  the  objectives  of  the  test  was  to  perform  an  operational  calibration  of  elements  of  interior 
orientation,  it  was  necessary  to  provide  extremely  accurate  independent  determinations  of 
exposure  stations.  Because  only  steilar-oriented  ballistic  cameras  could  provide  the  needed 
accuracies,  this  meant  that  the  flight  testing  had  to  be  conducted  at  night.  This  in  turn  meant 
that  the  portion  of  the  range  to  be  used  had  to  be  converted  temporarily  into  a  night  photogram- 
metric  range  through  the  installation  suitable  light  sources  over  selected  targets.  The  layout 
of  the  range,  the  intended  flight  paths  and  the  locations  of  ballistic  camera  tracking  stations 
are  shown  in  Figure  1 . 

In  the  designed  test,  four  flight  paths  cross  the  center  of  the  range  from  the  following 
directions; 


Path  1 
Path  2 
Path  3 
Path  4 


south  to  north, 
east  to  west, 
north  to  south, 
west  to  east. 


At  intervals  of  approximately  two  thousand  feet  throughout  each  crossing,  synchronized  photo¬ 
graphs  are  taken  of  the  range  by  the  pair  of  KC-6A  Cameras  carried  in  the  USQ-28  RC-135  air¬ 
craft.  At  the  center  of  each  exposure,  a  small  strobe  unit  (positioned  close  to  the  prime  camera) 
flashes.  The  five  most  nearly  central  flashes  on  each  flight  path  are  recorded  against  the  stellar 
background  by  at  least  two  of  three  ballistic  cameras  located  at  known  stations  on  the  range 
&ee  Figure  1),  Aircraft  velocity  is  approximately  600  feet  per  second  and  aircraft  altitude  is 
approximately  12,500  feet  above  mean  terrain. 

The  test  program  was  successfully  executed  in  March  1967.  The  actual  test  conformed 
closely  to  the  design.  Army  personnel  made  arrangements  with  local  power  companies  to  provide 
an  electrical  outlet  at  each  of  the  51  selected  targets  and  at  each  of  th’e  three  ballistic  camera  sta¬ 
tions.  The  incomplete  symmetry  in  the  pattern  of  selected  targets  (Figure  1)  is  attributable  to 
unavailability  of  power  in  some  creas  and  to  absence  of  targets  in  others.  Light  sources  con- 
sis.ed  of  500  watt  quartz  iodide  tamps  having  an  output  of  10,000  lumens  and  measuring  1/2  inch 
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in  diameter  by  3|  inches  in  length.  This  type  of  lamp  had  been  selected  from  a  variety 
of  potential  light  sources  on  the  basis  of  preliminary  aircraft  tests.  Each  lamp,  mounted  on  a 
small  shot  bog  as  a  base,  was  carefully  positioned  over  the  center  of  each  tatget. 

A  small  strobe  lamp  was  installed  in  one  of  the  three  camera  windows  at  precisely  known 
offsets  from  the  primary  and  secondary  cameras.  (Offsets  were  41  inches  from  the  primary  camera 
and  71  inches  from  the  secondary  camera.)  The  lamp  was  synchronized  to  flash  at  the  center  of 
each  shutter  opening. 

Ballistic  camera  support  for  the  tests  was  provided  by  DBA  personnel  operating  three 
600mm  f/3.5  ballistic  cameras  leased  from  Space  Systems  Laboratories.  The  camera  shutters 
were  driven  by  DBA  Digital  Programmer  Clocks  synchronized  with  WWV  to  within  one  millisecond. 
Kodak  microflat  plates  coated  with  103F  emulsion  were  employed  for  the  photography. 

Successful  tests  were  conducted  on  the  nights  of  March  24  and  29.  Of  the  four  orthogonal 
passes  for  each  night,  three  turned  out  to  lie  within  one  half  mile  of  the  designated  lines  and  one 
turned  out  to  be  almost  a  mile  off.  Because  of  such  difficulties  of  navigation,  only  eleven  of  the 
twenty  exposure  stations  for  each  night  were  actually  triangulated  by  ballistic  cameras.  This, 
however,  did  not  seriously  compromise  the  reduction,  inasmuch  as  triangulation  of  but  a  single 
exposure  station  would  theoretically  have  been  sufficient  for  the  recovery  of  the  elements  of 
interior  orientation. 

Altogether,  three  KC-6A‘s  provided  photography  over  the  improvised  night  photogrammetric 
test  range.  On  the  first  test,  cameras  005  and  006  were  flown  and  on  the  second,  cameras  006 
and  008  were  flown.  Thus,  camera  006  participated  in  both  tests.  In  the  first  test,  006  occupied 
the  primary  mount  (i.e.,  at  the  station  providing  the  Hypernas  indication  of  vertical ity),  and  in 
the  second  test)  it  occupied  the  secondary  mount.  Unfortunately,  the  Hypernas  verticality  sub¬ 
system  did  not  function  properly  on  either  of  the  two  tests.  Accordingly,  the  desired  evaluation 
of  the  accuracy  of  the  verticality  subsystem  could  not  be  accomplished.  Except  for  this,  the 
flight  test  program  was  a  complete  success. 


Hr 
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3*2.2  Preliminary  Corrections 

DBA  had  full  responsibility  for  the  data  reduction  of  ballistic  camera  plates  and  aerial 
photographs*  All  measurements  were  made  on  a  fully  calibrated  Mann  422G  Comparator  owned 
by  DBA.  The  ballistic  camera  reductions  turned  out  to  be  routine  with  rms  closures  of  triangu- 
lation  of  2  to  3  seconds  of  arc  being  attained.  Such  closures,  propagated  into  triangulated 
cooidinates  of  flashes,  lead  to  one  sigma  accuracies  not  exceeding  0.2  feet  in  X,  Y,  or  Z. 

This  positional  accuracy  is  equivalent  to  better  than  1:60,000  of  the  flying  height  or  to  better 
than  2.5  microns  on  a  given  aerial  photograph.  Inasmuch  as  eleven  of  each  group  of  twenty 
exposure  stations  were  triangulated,  the  net  effect  of  errors  of  ballistic  camera  triangulation 
on  the  determination  of  Xp,  yp,  c  can  be  expected  to  be  well  under  one  micron. 

Before  the  SMAC  reductions  could  be  performed,  a  number  of  small  preliminary  corrections 
had  to  be  applied  to  the  data.  Aside  from  comparator  corrections,  these  consisted  of: 

c  c  c 

(a)  offset  corrections  AX  ,  AY  ,  AZ  to  the  triangulated  positions  of  the  flashes  to 
establish  the  coordinates  of  the  exposure  stations  of  the  aerial  cameras, 

(b)  corrections  to  plate  coordinates  to  account  for  effects  of  atmospheric  refraction, 

(c)  corrections  to  plate  coordinates  to  account  for  effects  of  refraction  of  aircraft 
window, 

(d)  corrections  to  plate  coordinates  to  account  for  effects  of  film  deformation  by 
means  of  measurements  made  of  reseau  images. 

Inasmuch  as  the  flash  lamp  was  located  directly  behind  the  primary  and  secondary  cameras, 
the  direction  cosines  of  the  straight  line  joining  the  two  cameras  and  the  flash  lamp  are  defined 
by 

X  =  sin  H  cos  6 

p  =  cos  H  cos  6 

l>  -  sin  5 

in  which  H  is  aircraft  heading  measured  clockwise  from  north  (the  adopted  direction  of  the 
positive  Y  axis  in  a  plane  tangent  to  the  spheroid  at  the  center  of  the  range)  and  6  is  the  pitch 
angle  of  the  aircraft  (positive  when  nose  is  up).  Precise  values  of  H  and  S  are  provided  by  the 
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navigation  system.  The  offset  corrections  to  be  added  to  the  tricngulated  coordinates  of  a  flash 
ore  given  by 

AX  =  DX 
AY  *=  Dp 
AZ=  Du 

in  which  D  =  3.63  ft.  for  the  primary  camera  and  D  ~  5.88  ft.  for  the  secondary  camera.  The 
error  in  these  corrections  is  unlikely  to  exceed  ±  0.1  feet. 

It  is  readily  shown  by  ray  tracing  based  on  a  flat  earth  model  that  the  corrections  to  be 
applied  or  added  to  the  plate  coordinates  of  a  vertical  photograph  of  flat  terrain  to  account  for 
the  combined  effects  of  atmospheric  refraction  and  'window1  refraction  are  given  by  the  expressions 


8x  -  Gx 
6y  -  Gy 

In  which 

G  =  0  +*¥-) 

wherein 

C"  principal  distance  of  camera, 
x,y—  plate  coordinates  referred  to  principal  point, 

=  Index  of  refraction  of  light  at  ground  level, 
jl  =  index  of  refraction  at  camera. 

The  value  of^  is  computed  from  the  formula 

(jj-1)  106  *=  {77.34  +  )  y+273  ' 

in  which 

P  =  pressure  in  millibars, 

T  =  temperature  in  degrees  Centigrade, 

X  *  wave  length  of  light  in  microns  (approximately  0.55  microns  for  the  middle  ot 
the  visible  ronge). 


When  the  camera  is  in  a  pressurized  compartment  (os  in  the  USQ-28  System),  the  index  p  applies 
to  the  compartment  and  not  to  the  atmosphere  immediately  outside  the  aircraft.  It  follows  that 
when  air  density  in  the  camera  compartment  is  equal  to  tfat  at  ground  level,  p  =  pB  and  the 
result  is  that  atmospheric  refraction  and  window  refraction  perfectly  cancel  each  other.  We 
point  this  out  to  emphasize  the  consequences  of  the  camera  window  in  a  pressurized  aircraft. 
Other  treatments  of  refraction  that  we  have  encountered  (e.g.,  those  in  the  Manual  of  Photo- 
grammetry)  fail  to  distinguish  between  the  refractive  effects  of  pressurized  and  unpressurized 
camera  compartments. 

The  last  of  the  preliminary  corrections,  the  reseau  corrections  for  film  deformation, 
were  applied  to  each  frame  selected  to  be  measured  for  the  SMAC  reduction.  The  entire  set 
of  37  reseau  images  was  measured  on  each  selected  frame  along  with  images  of  fiducials  and 
control  points.  Corrections  for  film  deformation  employed  the  models: 

Ax  =  xr-x  =  a0+ajx+a2y  +  a3  x2  +  ajxy  +  a$y2 

+  a4x3  +  a2x2y+  a8xy2  +  a^y3 

Ay  =  yr-y  =  b0  +  bjx  +  b2y  +  b3x2  +  b4xy  +  bjy2 

4  b^x3  +  br  x2y+  b8xy2  +  b^y3 

in  which 

Ax,  Ay  =  corrections  to  be  added  to  measured  coordinates  x,y; 

xr  ,  yr  “  reseau  coordinates; 

at  ,  bt  =  coefficients  determined  from  least  squares  fit  to  discrepancies  between 
pre-established  coordinates  and  measured  coordinates  of  reseau  images. 

Several  trials  of  the  fitting  process  were  made  to  determine  the  set  of  coefficients  providing 
a  compact  model,  i-e.,  one  containing  no  .inessential  coefficients.  In  boih  the  Ax  and  Ay 
models,  it  was  found  that  the  coefficients  of  x3  and  y3  could  be  consistently  dropped  without 
adversely  affecting  Jhe  quadratic  form  of  the  residuals.  All  other  coefficients  were  found  to 
contribute  significantly  to  the  fit.  Accordingly,  in  the  final  reseau  reductions  only  8  of  the  10 
coefficients  in  each  model  were  exercised. 
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Initial  fits  of  the  models  for  deformation  yielded  rms  eirors  on  the  order  of  5  to  6 
microns  for  each  frame,  in  the  light  of  our  previous  experience,  this  was  not  considered 
to  be  satisfactory.  Upon  examining  the  residuals,  we  found  that  a  few  points  had  large 
residuals  (over  10  microns)  of  nearly  constant  magnitude  on  frame  after  frame.  When  these 
points  were  dropped  from  the  reduction,  rms  errors  were  reduced  to  about  half  their  previous 
values.  This  suggests  that  either  the  pre -calibrated  values  for  the  coordinates  of  these  particular 
points  were  of  marginal  accuracy  or  else  thot  the  individual  projection  units  for  these  points  had 
become  slightly  misaligned  after  calibration. 

In  Table  1  we  have  listed  the  standard  deviations  of  the  reseau  fit  for  each  of  the  twenty 
frames  reduced  from  camera  006  on  each  of  the  two  flights.  All  37  reseau  images  were  used, 
with  the  coordinates  of  the  poor  reseau  points  being  replaced  by  adjusted  values.  On  Test  No.  1 
the  fit  of  the  model  was  excellent  with  typical  rms  errors  being  2.0  microns  in  x  and  2.4  microns 
in  y.  Somewhat  poorer  fits  v/ere  generally  obtained  on  Test  No.  2,  the  typical  rms  errors  being 
2.3  microns  in  x  and  3.6  microns  in  y.  Inasmuch  as  the  y  axis  runs  laterally  across  the  film,  the 
greater  rms  error  in  y  in  both  tests  is  probably  attributable  to  the  occasionally  severe  lateral  defor¬ 
mation  occurring  near  the  edges  of  the  film.  This  is  borne  out  by  the  finding  that  the  largest  y 
residuals  tend  consistently  to  correspond  to  points  closest  to  the  edges  of  the  film.  Residuals  in 
x,  on  the  other  hand,  display  no  such  tenden-y.  When  the  fitting  is  limited  to  the  25  interior 
reseau  points  (spaced  at  2  inch  intervals)  results  on  the  poorer  frames  (e.g.  frames  53,  153,  in 
Table  I)  become  significantly  improved  (to  about  the  2  r  1  -ron  level).  In  applications  to  routine 
aerotriangulation,  it  would  probably  be  best  to  exercise  only  the  interior  reseau  points  in  the 
correction  for  film  deformation  (here,  correction  of  the  fiducials  for  deformation  is  of  no  practical 
consequence,  for  the  recovered  horizontal  coordinates  of  the  exposure  station  provide  effective 
projective  compensation  for  small  errors  in  the  adopted  coordinates  of  the  principal  point).  In 
the  application  to  SMAC,  on  the  other  hand,  the  fiducials  should  define  a  consistent  frame  center 
so  that  the  calibrated  coordir.  'es  of  the  principal  point  are  referred  to  a  common  origin.  It  was 
for  this  reason  that  v.*e  employed  the  entire  set  of  reseau  images  on  each  frame  in  establishing  the 
coefficients  of  the  model  for  film  deformation. 
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TABLE  1.  Results  of  Fitting  of  Models  for  Film  Deformation  of  Frames  from  Camera  00 (> 


TEST  NO.  1,  24 March  1967 

TEST  NO 

.  2,  29  March  1967 

Frame 

No. 

(Microns) 

ar 

(Microns) 

Frame 

No. 

0* 

(Microns) 

°r 

(Microns) 

194 

1.7 

2.0 

118 

1.6 

2.5 

195 

1.2 

2.3 

119 

1.4 

4.3 

196 

1.9 

2.1 

120 

1.8 

4.1 

197 

1.7 

2.1 

121 

2.2 

1.9 

198 

1.5 

2.1 

122 

1.7 

3.8 

222 

1.3 

2.0 

153 

3.9 

6.3 

223 

1.5 

2.5 

154 

1.8 

2.1 

224 

1.9 

2.4 

155 

2.4 

2.6 

225 

2.7 

2.1 

156 

2.2 

2.6 

226 

1.9 

1.9 

157 

1.7 

1.9 

282 

1.7 

2.2 

51 

2.3 

4.0 

283 

2.2 

2.2 

52 

2.0 

2.7 

284 

1.6 

2.1 

53 

3.2 

7.6 

285 

3.3 

4.5 

54 

2.0 

3.1 

286 

1.9 

2.6 

55 

2.5 

3.8 

310 

2.7 

2.3 

86 

2.2 

3.0 

311 

2.4 

2.1 

87 

1.9 

3.4 

312 

1.4 

2.1 

88 

2.0 

2.9 

313 

2.2 

2.2 

90 

2.1 

2.3 

314 

1.8 

2.7 

91 

1.7 

2.8 

Grand  RMS  2.0 

|  Grand  RMS  2.3 

3.6 
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3*2.3  Resuits  of  Aerial  SMAC  Calibration  of  Camera  006 

The  primary  calibration  to  be  performed  from  the  data  gathered  on  the  McLure  test 
was  that  for  Camera  006.  Although,  Cameras  005  and  008  were  also  flown,  they  were  under 
engineering  evaluation,  whereas  Camera  006  had  been  officially  delivered  to  the  Air  Force. 
Accordingly,  our  major  concentration  of  effort  was  on  006  and  it  is  the  results  for  this  camera 
that  we  shall  discuss  in  greatest  detail. 

The  five  most  central  frames  on  each  of  the  four  passes  from  the  two  tests  were  selected 
for  measurement.  All  target  images  (typically  about  30  per  frame),  all  37  reseau  images,  and 
the  images  of  the  primary  and  secondary  fiducials  were  measured  on  each  frame,  double  settings 
being  made  on  each  point.  The  quality  of  all  images  was  excellent.  In  particular,  images  of 
the  targets  throughout  the  format  wore  decidedly  better  than  any  stellar  images  that  we  had 
examined  on  plates  taken  by  commercial  mapping  cameras. 

Each  of  the  two  sets  of  twenty  frames  was  subjected  to  an  independent  ShfaC  reduction 
In  order  to  ascertain  the  consistency  of  the  results  from  one  flight  to  another.  Totals  of  574  images 
and  640  images  were  measured  on  the  first  and  second  sets  of  frames,  respectively.  Trial  runs 
established  that  two  coefficients  of  radial  distortion  (K, ,  Ka)  ond  two  coefficients  of  decentering 
distortion  (Pj ,  P2)  provided  an  adequate  model  (exercise  of  higher  order  coefficients  produced 
essentially  no  reduction  in  the  quadratic  form  of  the  residuals).  Ballistic  camera  positions, 
corrected  for  offset  between  flashing  light  and  camera,  were  exercised  as  apriori  constra«rrs 
in  the  SMAC  reduction  on  the  eleven  frames  in  each  set  that  were  successfully  observed. 

Although  propagation  of  closures  of  ballistic  camera  triangulaiior*  indicated  one  sigma  accuracies 
of  better  then  0.2  feet  in  X,  Y,  and  Z,  the  a  priori  constraints  actually  used  in  the  SMAC  reduc¬ 
tion  were  relaxed  by  a  factor  of  5.  This  was  done  to  avoid  the  possibility  of  overconstraining 
in  view  of  the  consideration  that  errors  in  die  given  survey  of  ballistic  camera  stations  could  con¬ 
ceivably  introduce  a  bias  of  as  much  as  0.5  feet  into  the  triangulations  of  the  flashes. 


The  distortion  functions  resulting  from  the  two  independent  aerial  SMAC  reductions 
are  plotted  in  Figure  2  and  are  tabulated  in  Table  2  (one  sigma  error  bounds  are  indicated 
in  both  presentations).  The  calibrated  values  of  the  interior  projective  parameters  are  also 
tabulated  in  Table  2.  The  rms  error  of  the  plate  measuring  residuals  turned  out  to  be  3.3 
microns  for  the  first  set  of  frames  and  3*7  microns  for  the  second.  The  somewhat  poorer  result 
for  the  second  set  is  probably  attributable  in  part  to  the  somewhat  pocrer  corrections  obtained 
for  film  deformation  (review  Table  1).  The  error  bounds  in  Figure  2  and  the  standard  deviations 
given  in  Table  2  are  based  on  the  rms  errors  of  the  residuals  from  the  respective  adjustments 
(i.e.,  3.3  and  3.7  microns). 

We  see  from  the  plotted  and  tabulated  results  that  the  two  calibrations  are  mutually 
consistent  considering  their  standard  deviations.  The  rms  disrecpancy  between  the  two  radial 
distortion  functions  is  under  one  micron,  and  the  maximum  discrepancy  of  1.8  microns  (at  r=  120mm) 
is  not  out  of  line  with  the  sigmas  of  the  two  curves.  Decentering  distortion  from  both  calibrations 
is  unusually  low,  amounting  to  only  1 .4  and  2.3  microns,  respectively,  at  the  extremeties  of  the 
format.  Although  the  phase  angles  0O  from  the  two  calibrations  differ  significantly  from  each 
other,  this  is  of  no  practical  consequence  inasmuch  as  both  decentering  profiles  are  so  close  to 
zero. 

Distortion  curves  generated  by  the  aerial  SMAC  calibrations  of  Cameras  005  and  008 
are  presented  «n  Appendix  B. 

3.2.4  Discussion  of  Observed  Shift  of  Principal  Point 

An  unexpected  finding  concerns  the  magnitude  of  the  coordinates  of  the  principal  point 
(Table  2).  Both  calibrations  agree  to  within  uncertainties  of  a  few  microns  that  Xp  is  well  over 
200  microns  and  yp  is  about  100  microns.  However,  these  values  are  inconsistent  with  the  labor¬ 
atory  calibration  and  with  the  stellar  calibration,  both  of  which  are  much  closer  to  zero  (these 
will  be  taken  up  later).  The  aerial  SMAC  calibrations  for  Cameras  005  and  008  also  recovered 
large  values  for  Xp,  yp  (namely,  Xp  -  .167  mm,  yp  =  .045  mm  for  005,  and  Xp  =  .175  rr.m, 
yp  =  -.052  mm  for  008).  Naturally,  we  sought  an  explanation  r*?r  such  large  discrepancies  in 
xp,yp.  Several  hypothetical  explanations  were  given  consideration,  namely. 
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10/li  J  Cr  RADIAL  DiSTORliON,  TEST  \  (c  =  i5i.z62mm) 


Figure  2.  Radio!  distortion  curves  and  decentering  profiles  resulting  from  Aerial  SMAC 
calibrations  of  Camera  006  for  two  independent  flight  tests  (accompanied  by  one  sigma 
error  bounds). 
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TABLE  2.  Summary  of  Results  of  Aeriol  5MAC  Calibrations  of  Camera  006 


TEST  NO.  1 

TEST  NO.  2 

Parameter 

Value 

Parameter 

Value 

yp  (mm) 
c  (mm) 

Kj  (mm/mm3) 
Ks(rnm/mmB) 

(mm/mm2) 

(deg) 

.231*. 003 
.104.003 

151.23 1±.002 

(~.274.015)x  10*7 
(  .73Qfc.065)x  10“  * 

(  .706*. 356) x  10“7 
327.9*28.4 

xp  (mm) 
yp  (mm) 
c  (mm) 

Kj  (inm/mm3) 

K8  (mm/mm6 ) 

Jj  (mm/mm2) 

$o  (deg) 

. 222* . 004 
.099-. 004 
151.2274..  002 

(-.214±.015)x  10*7 
(  .509*.068)xl0"i2 

(1.192*.  363)  x  10“7 
193.1*16.4 

POINTS  ON  RADIAL  AND  DECENTERING  DISTORTION  CURVES 


TEST  NO.  1 

TEST  NO. 

2 

r 

«,w 

P, 

r 

K 

(rnm) 

(Microns) 

(Microns) 

(mm) 

(Microns) 

(Microns) 

0 

0 

0 

0 

0 

0 

20 

3.9*0.01 

0.040.01 

20 

3.40.01 

0.05*0.01 

40 

6.6±0.09 

0. 1 1*0.  06 

40 

5.9*0.09 

0.19*0.06 

60 

7.40.27 

0.25±0. 13 

60 

6.6*0.27 

0.43*0.13 

80 

4.9*0.55 

0.45±0.23 

80 

5.1*0.56 

0.76*0.23 

100 

0.5*0.85 

0.71*0.36 

100 

1.7*0.87 

1.19*0.36 

120 

-4.5*1.07 

1.02*0.51 

120 

-2.7*1.08 

1.72*0.52 

140 

-7.2*1. 17 

1.39*0.70 

140 

-6.2*1.18 

2.340.71 

(a)  referred  to  value  of  c  of  151.262mm 

(b)  referred  to  value  of  c  of  151.255mm 
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1 .  Refractive  effects  of  the  shockwave  of  the  aircraft; 

2.  Prism  effect  of  the  camera  window; 

3.  Biases  in  ballistic  camera  triangulaticns; 

4.  Incorrectly  applied  offset  corrections  between  triangulated  flashes  and  cameras; 

5.  Timing  bias  in  occurrence  of  flash  relative  to  exposure  of  frame; 

6.  Tilts  in  the  platen  during  IMC  operation. 

Analysis  eliminated  the  first  four  of  these  hypothesis.  The  first  hypothesis  wos  eliminated  by  three 
considerations:  (a)  'worst  case*  ray  tracing  through  a  prismatic  analog  of  a  shock  wave  indicated 
an  effect  more  than  an  order  of  magnitude  smaller  than  that  observed;  (b)  the  effect,  furthermore, 
would  cause  a  displacement  in  Xp  opposite  in  direction  to  that  observed;  (c)  lateral  symmetry  of 
the  shock  wave  would  preclude  any  effect  on  yp. 

Prism  effect  of  the  window  attributable  to  lack  of  parallelism  of  front  and  rear  surfaces 
would  have  to  be  equivalent  to  0  wedge  angle  of  about  0?1  in  order  to  account  fer  a  200  micron 
displacement  of  principal  point.  This  would  be  about  150  times  greater  than  specifications  for 
the  window.  Moreover,  prism  effect  would  not  explain  the  observed  values  of  yp  which  are  opposite 
sign  for  cameras  005  and  008  relative  to  006. 

Inasmuch  as  the  flying  height  is  about  12,000  feet  above  ground  level,  a  bias  in  the 
triangulated  coordinates  of  the  flashes  of  12000  x  (0.2/150)  =  16  feet  would  be  needed  to  explain 
a  shift  of  200  microns  in  the  principal  point.  This  is  almost  one  hundred  times  greater  than  the 
standard  deviations  of  the  triangulated  coordinates  and  corresponds  to  an  error  of  about  0.5  mm 
Oi*  the  ballistic  camera  plate.  Above  and  beyond  this,  because  of  the  design  of  die  flight  test, 
a  bias  in  ballistic  camera  triangulation,  if  constant  for  all  flight  paths,  would  have  no  ultimate 
effect  on  the  determination  of  the  principal  point.  Thus  a  constant  bias  of  say  16  feet  in  Y  (north) 
would  produce  a  bias  of  plus  200  microns  in  Xp  in  the  separate  reduction  of  photos  taken  on  Hie 
south  to  north  flight  and  o  bias  of  minus  200  microns  in  Xp  in  the  separate  reduct  ion  of  photos  taken 
on  the  north  to  south  flight.  Similarly,  it  would  produce  biases  of  plus  and  minus  200  microns  in 


yp  on  the  separate  reductions  of  west  to  east  and  east  to  west  flights,  respectively.  Accordingly, 
when  all  four  legs  are  reduced  simultaneously,  the  effects  of  a  constant  bias  of  triangulation 
mutually  cancel  one  another  insofar  as  the  net  displacement  of  the  principal  point  is  concerned. 
This,  in  fact,  was  a  guiding  consideration  in  the  design  of  the  flight  test. 

Assuming  Hie  triangulated  positions  of  the  flashes  to  be  correct,  we  find  that  the  misappli¬ 
cation  of  offset  corrections  could  only  partially  explain  the  observed  discrepancies.  If  the  cor¬ 
rections  were  not  applied  at  all,  the  bias  to  be  expected  in  Xp  would  amount  to  plus  43  microns 
when  the  prime  station  is  occupied  and  to  plus  75  microns  when  the  alternate  station  is  occupied. 

If  they  were  applied  with  the  wrong  sign,  the  biases  would  be  doubled  in  magnitude  and  changed 
in  sign  (i.e.,  they  would  become  minus  86  and  minus  150  microns,  respectively).  It  follows  that 
the  observed  values  of  Xp  of  over  plus  200  microns  with  the  camera  in  both  primary  and  alternate 
positions  cannot  be  explained  in  terms  of  any  likely  blunder  in  application  offset  corrections. 
Moreover,  a  thorough  investigation  cf  the  matter  showed  that  the  corrections  were  indeed  correctly 
applied. 

The  flashing  light  is  supposed  to  be  synchronized  to  flash  with  the  exposure  of  the  primary 
fiducials  which  in  turn  is  supposed  to  occur  at  the  center  of  exposure  of  each  frame.  Let  us  postu¬ 
late  that  the  primary  fiducials  were  properly  synchronized  with  the  exposures  but  that  the  flashing 
light  was  not  properly  synchronized  with  the  primary  fiducials.  To  analyze  this  situation  let  us 
assume  perfectly  vertical  photography  at  a  constant  altitude  H  above  a  flat  earth.  Then  we  can 
write  the  projective  equations  in  the  simplified  form 


(86) 


x  -  Xp  =  ~r  [  (X  -  XC)  cos  K  +  (Y  -  Y0)  sin  K  3 


y-yP 


■|j  L-(X-  XC)  sin  k  +  (Y  -  Y°)  cosk  ] 


in  which  K  denotes  the  swing  angle.  If  these  are  viewed  as  observational  equations  for  the 
recovery  of  Xp  and  yp,  we  see  that  errors  SXC,  SYC  in  Hie  coordinates  of  the  exposure  station 
will  result  in  errors  in  Xp  and  yp  of 
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(87) 


6xp  =  C(cos/c)8XC  +  (sinK)6YC], 
8yp  =  -pj  C-(sin/£)5XC  +  (cos  k)  8YC  ]. 


If  we  now  assume  that  errors  In  8XC,  5YC  are  attributable  to  a  timing  offset  8t  between  the 
exposures  of  the  frames  and  the  flashing  of  the  strobe,  we  have 

=  X  St 

=  Y  8t 

"  t  -  t0  -  (time  of  flash)  -  (time  of  frame  exposure). 

if  we  further  assume  the  total  velocity  V  of  the  aircraft  is  constant  and  that  the  positive  x  axis 
of  the  plate  is  always  aligned  with  the  direction  of  flight,  the  X,  Y  components  of  velocity 
become  simply: 


(88) 


SX 
8Y° 
in  which 
8t 


m 


X  “V  cos  K  , 

*  =  V  si  rut. 

Substitution  of  (  89  )  into  (  88  ),  followed  by  substitution  of  the  resulting  expressions  into  (  87  ) 
yields  the  results 


(90) 


8xp  H  6t  , 
6y  =  ° 


These  expressions  are  independent  of  the  direction  of  flight.  It  follows  that  the  exercise  of  direct 
and  reverse  flight  paths  does  not  automatically  cause  the  effects  of  a  constant  bias  in  synchroni¬ 
zation  to  cancel  out  as  is  the  case  with  a  constant  bias  in  position.  It  is  clear  that  such  cancellation 
could  have  been  achieved  from  direct  and  reverse  flight  paths  if  the  camera  w<»re  rotated  180°  in 
swing  angle  between  direct  and  reverse  paths  (on  the  other  hand,  such  a  procedure  would  undo  the 
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the  automatic  cancellation  of  the  effects  of  constant  bias  in  triangulation).  Inasmuch  as 
c«*  150  mm,  H  a  12,000  ft.,  V  «600  ft  ./sec.  for  the  McLure  tests,  it  follows  from  (  90  ) 
that  to  account  for  an  error  of  Sxp  =  0-2  mm,  an  error  of  synchronization  of  6t  -  .028  sec.  is 
needed.  When  image  motion  compensation  (IMC)  is  operative,  as  it  was  In  the  McLure  tests, 
an  error  of  synchronization  projectively  equivalent  to  that  just  considered  is  one  in  which  the 
flashing  light  is  perfectly  synchronized  with  the  exposure,  but  the  primary  fiducials  are  exposed 
.028  sec.  too  late.  A  third  possibility  is  for  the  flashing  light  to  be  perfectly  synchronized  with 
the  primary  fiducials,  which,  in  turn,  are  not  perfectly  synchronized  perfectly  with  the  exposure 
of  the  frame.  In  this  situation,  the  expression  for  fixp  would  hecome 

(?1)  fixp  »  2^  St, 

and  an  error  of  6t  =  .014  sec.  would  be  sufficient  to  explain  a  200  micron  shift  in  Xp. 

From  the  information  available  to  us,  we  are  unable  to  rule  out  any  of  the  three  possibilities 
just  cited,  namely: 

1 .  primary  fiducials  perfectly  synchronized  with  exposure  of  frame,  but  strobe  flash 
occurring  about  28  milliseconds  too  late; 

2.  strobe  flash  perfectly  synchronized  with  exposure  of  frame,  but  primary  fiducials 
exposed  about  28  milliseconds  too  late; 

3.  strobe  flash  and  exposure  of  primary  fiducials  perfectly  synchronized,  but  occurring 
14  milliseconds  after  exposure  of  frame. 

While  each  of  these  could  explain  a  200  micron  shift  in  Xp,  we  have  no  way  of  verifying  which, 
if  any,  is  correct.  We  do,  however,  have  solid  evidence  that  something  is  a.niss  ?n  the  timing 
of  events.  By  measuring  the  distances  between  the  primary  and  secondary  fiducials  on  each  frame 
and  comparing  them  with  their  calibrated  values  (no  IMC),  we  were  able  to  ascertain  the  preci:e 
amount  of  fiim  movement  occurring  between  their  respective  exposures.  This  averaged  93  microns 
on  the  first  flight  test  and  87  microns  on  the  second.  Specifications  for  the  camera  call  for  the 
following  sequence  of  events: 


1 .  well  prior  to  initiation  of  the  shutter  pulse,  the  platen  is  brought  up  to  full  speed 
in  accordance  with  the  V/H  input; 

2.  at  the  initiation  of  the  shutter  pulse,  the  secondary  fiducials  are  exposed  and  the 
shutter  begins  to  open; 

3*  precisely  5  milliseconds  later  at  the  center  of  the  exposure  (which  has  a  constant 
duration  of  10  milliseconds)  the  primary  fiducials  are  exposed; 

4.  5  milliseconds  after  this,  the  shutter  is  fully  closed* 

In  the  McLure  tests,  the  IMC  rate  ideally  should  have  been  about  c  (V/H)  «7. 1  milli- 
meters/second.  Accordingly  the  amount  of  film  movement  in  the  5  milliseconds  between  exposure 
of  secondary  and  primary  fiducials  should  have  been  about  36  microns*  As  has  already  been  pointed 
out,  it  was  in  fact  87  and  93  microns,  on  the  average,  for  the  two  tests.  This  suggests  that  either 
the  IMC  rate  was  in  error  by  a  factor  of  about  2*5  or  else  the  interval  between  exposures  of  secon¬ 
dary  and  primary  fiducials  was  about  2*5  times  longer  than  the  5  milliseconds  specified.  Inasmuch 
cs  the  images  of  the  ground  lights  were  perfectly  round  and  not  the  least  elongated,  we  may  con¬ 
clude  that  the  IMC  rate  was  not  ot  fault  and  that  the  time  interval  between  exp^ures  of  secondary 
and  primary  fiducials  therefore  must  have  been  about  13  milliseconds.  The  finding  that  the  actuc' 
interval  between  one  pair  of  events  is  definitely  incorrect,  lends  support  to  the  hypothesis  that  a 
deficiency  in  synchronization  is  a  iikely  explanation  for  the  observed  displacement  of  the  principal 
point  in  Xp.  Recall  that  a  delay  of  14  milliseconds  in  the  exposure  of  the  primary  fiducials  relative 
to  the  exposure  of  the  frame  would  explain  a  200  micron  shift  in  Xp  if  the  flashing  light  were  truly 
synchronized  with  the  primary  fiducial  (as  it  is  supposed  to  be).  In  view  of  the  similar  interval 
(13  milliseconds)  beh  >n  exposure  of  secondary  and  primary  fiducials,  a  natural  conjecture  is  that 
it  war  the  exposure  o  e  secondary  fiducials  that  actually  coincided  with  the  center  of  the  frame 
exposure.  Whether  or  not  this  is  the  cose,  it  is  clear  that  a  thorough  investigation  should  be  made 
into  the  timing  and  synchronization  of  events  occurring  throughout  the  exposure  cycle  of  the  KC-6A* 

*  After  the  above  was  written  an  oscilloscope  test  performed  on  Camera  006  by  the  1370  Photo 
Squadron  showed  that  the  flashing  light  actually  occurred  between  25  and  30  milliseconds 
after  the  midpoint  of  the  exposure.  Thus  the  first  hypothesis  forwarded  above  is  now  confirmed. 
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An  error  in  synchronization  does  not  serve  to  explain  the  observed  displacement  of 
about  100  microns  in  yp.  If  the  platen  should  somehow  become  tilted  laterally  by  about 
2  minutes  of  arc  during  the  IMC  movement,  such  a  displacement  in  y^  wou'1  cccur.  The 
effects  of  such  a  tilt  on  focus  would,  at  worst,  be  only  barely  noticecfble.  If,  for  instance, 
the  edge  of  the  frame  were  75  microns  from  the  plane  of  best  focus,  the  resulting  'circle  of 
confusion1  of  an  image  patch  near  the  edge  would  have  a  diameter  of  only  15  microns  (for 
the  camera  aperture  at  f/5).  This  would  nor  perceptibly  degrade  the  images  of  the  ground 
targets,  for  these  are  about  50  microns  in  diameter.  For  this  reason,  examination  of  image 
quality  does  not  provide  a  means  for  either  confirming  or  denying  the  existence  of  a  tilt 
of  the  platen  as  small  as  two  minutes  of  arc.  Although  we  are  unable  at  this  time  to  suggest 
any  other  plausible  explanation  for  the  displacement  in  yp,  this  naturally  does  not  rule  out 
the  existence  of  alternative  explanations. 

Considering  the  possible  implication  of  IMC  in  the  shift  of  the  principal  point,  we 
would  urge  that  in  any  future  aerial  SMAC  calibrations  of  the  KC-6A,  the  basic  flight  test 
pattern  be  flown  both  with  and  without  the  operation  of  IMC.  This  would  settle  many 
matters  that  are  now  subject  to  conjecture. 

Just  what  is  the  precise  photogrammetric  significance  of  the  elements  of  interior 
orientation?  In  view  of  our  findings  concerning  the  principal  point,  consideration  of  this 
matter  is  clearly  in  order.  In  conventional  aerotriangulation  over  flat  terrain  rather  large 
errors  in  Xp,  yp  and  c  are  readily  tolerated  inasmuch  as  almost  perfect  projective  compensation 
is  provided  by  the  recovered  coordinates  Xc  ,YC  ,ZC  of  the  exposure  sta’.on.  Here,  the  central 
problem  reduces  to  one  of  interpolation  within  a  framework  of  known  control  points  to  establish 
the  coordinates  of  other  points  of  interest.  Thus  the  process,  being  essentially  interpolative  in 
character,  is  not  normally  compromised  by  projective  tradeoffs  of  errors  in  Xp,  yp  and  c 
insofar  as  the  desired  eiH  results  are  concerned.  However,  this  comfortable  situation  no 
longer  applies  as  soon  as  the  photogrammetric  process  becomes  extrapolative,  as  is  the  case 
in  operations  exercising  the  full  capability  of  the  USQ-28  System.  Here,  external  sensors 
(Shiran,  Hypernas,  Terrain  Piofile  Recorder)  provide  sharp  constraints  on  elements  of  exterior 
orientation.  Theoretically,  this  permits  one  to  perform  aerotriangulation  without  recourse  to 
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ground  control.  However,  since  tight  constraints  on  elements  of  exterior  orientation 
sharply  limits  ihe  effectiveness  of  projective  compensation  for  defects  in  elements  of 
interior  orientation,  it  follows  that  far  higher  accuracies  are  required  of  elements  of 
interior  orientation  when  such  constraints  are  brought  to  bear.  While  errors  of  TOO  or 
200  microns  in  Xp,  yp,  c  are  of  no  practical  consequence  in  conventional  aerotrianguiation, 
such  errors  would  be  intolerable  when  accurate  external  sensors  are  exercised.  Full 
realization  of  the  potential  of  the  USQ-28  System  demands  that  errors  in  xn,  y  and  c  be 
suppressed  to  under  10  microns.  Hence,  our  concern  with  the  results  of  aerie!  SMAC  with 
regard  to  the  principal  point;  the  200  micron  shift  in  xp  would,  if  real  and  unaccounted  for, 
introduce  an  error  in  planimetry  of  H/750. 

Another  situation  in  which  errors  in  elements  of  interior  orientation  have  a  first  order 
effect  merits  discussion.  This  is  in  photogrammetric  resection  wherein  the  objective  is  the 
accurate  recover/  of  the  coordinates  Xc  ,YC,ZC  of  the  exposure  station.  This  approach  has 
been  used,  for  example,  for  the  evaluation  of  tracking  accuracies  of  Hiran  and  Shiran  from 
photography  taken  over  the  Phoenix  Test  Range.  For  H  =  30,000  feet,  an  error  of  200  microns 
in  principal  point  would  introduce  an  error  of  H/750  =  40  feet  in  the  horizontal  position  of  the 
camera  as  obtained  by  resection.  Such  an  error  would  clearly  invalidate  the  evaluation  of  any 
accurate  tracking  system.  It  is  to  be  emphasized  thet  whenever  the  coordinates  Xc ,  Yc  ,ZC 
recovered  from  a  photogrammetric  reduction  are  to  be  interpreted  as  having  physical 
significance,  the  errors  in  Xp,  yp,  c  assume  paramount  importance.  This  fact  seems  not  to 
have  been  given  due  consideration  in  certain  flight  testing  programs  which  may  have  been 
significantly  compromised  as  a  result. 


I  i 


3.2.5  Analysis  of  Residuals 

We  have  yet  to  consider  ir.  detail  the  observational  residuals  obtained  from  the  SMAC 
reductions.  These  should  ideally  be  completely  random  and  should  have  a  standard  deviation 
consistent  with  the  2  to  3  microns  typical  of  the  resea-j  reduction.  As  we  have  seen,  the 
standard  deviations  of  the  residuals  from  the  two  reductions  turned  out  to  be  3.3  and  3.7 
microns,  respectively,  or  somewhat  higher  than  desirable.  This  suggests  the  presence  of 
some  degree  of  unmodeled  systematic  error.  To  aid  in  the  examination  of  the  residuals. 
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we  have  provided  in  Figures  3  end  4  graphs  of  the  residual  vectors  fi  :~t  the  two  SfAAC 
reductions.  In  these  figures,  all  residuals  from  the  20  frames  carried  in  each  5MAC 
reduction  are  plotted  on  a  common  graph  against  the  plate  coordinates  of  the  Images, 

The  resulting  high  density  of  residuals  helps  to  isolate  any  localized  systematic  tendencies. 

We  see  from  Figures  3  and  4  that  definite  systematic  patterns  do  exist,  psrtlc  larly 
near  the  lower  left  hand  and  upper  right  hand  corners  as  well  os  near  the  right  hand  edge 
between  y  =  20  and  y  ~  70  millimeters.  A  fairly  good  correlation  Is  seen  to  exist  between 
the  more  pronounced  systematic  tendencies  in  both  figures.  Indeed,  when  both  figures 
are  superimposed,  os  is  done  in  Figure  5,  systematic  patterns  become  even  more  clearly 
defined. 

In  examining  the  tabulations  of  residuals,  we  noticed  that  the  residuals  for  certain 
control  points  appeared  to  be  consistently  larger  than  average  In  magnitude  and  to  be 
systematic  in  direction  within  a  given  run.  To  illustrate  this,  we  have  listed  in  Table  3 
the  residuals  for  point  no.  39  for  all  frames  on  both  tests.  It  is  clear  from  the  fable  that  the 
residuals  from  the  two  fesfs  are  strongly  correlated.  The  pattern  of  the  residuals  corresponds 
closely  to  that  to  be  expected  from  an  error  in  the  survey  of  the  control  point.  Thus 
significant  negative  errors  in  the  given  X  (east)  and  Y  (north)  coordinates  of  a  target  would 
tend  to  generate  positive  x  residuals  and  positive  y  residuals  on  west  to  east  runs;  positive  y 
residuals  and  negative  x  residuals  on  south  to  north  runs;  negative  x  residuals  end  negative  y 
residuals  on  east  to  west  runs;  and  negative  y  residuals  and  positive  x  residuals  on  north  to 
south  runs.  This  is  indeed  the  general  tendency  of  the  residuals  in  Table  3,  with  the  effects 
of  a  possible  error  in  X  predominating. 

In  view  of  the  fact  that  residuals  from  several  other  control  points  also  indicated  the 
presence  of  significant  survey'  error,  we  decided  to  examine  more  closely  the  matter  of  survey 
accuracies.  From  ETL  contacts  with  USCGS  we  had  been  informed  that  the  McLure  survey 
was  considered  to  be  a  third  order  survey  exercising  second  order  procedures.  Although  this 
was  deemed  tc  be  decidedly  marginal  for  the  SMAC  c jdibrafions,  we  reasoned  that  ultimate 
effects  of  survey  errors  on  interior  projective  parameters  would  be  largely  averaged  out  by 
virtue  of  the  orthogonal  pattern  of  direct  and  reverse  flights.  Inasmuch  as  the  basic  3MAC 
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Figure  3,  Composite  plot  of  residual  vectors  from  20  frames  carried  in  Aerial  SMAC 
calibration  of  Camera  006,  Flight  Test  No.  1  (nns  vaiue  of  residuals ”3.3  microns). 
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Figure  5.  Combined  plot  of  residual  vectois  in  Figures  3  and  4. 
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TA8LE  3.  Plate  residuals  for  Control  Point  No.  39  on  each  frame  from  separate 
SMAC  reductions  of  two  flight  tests  (s)  with  Grigs  no!  survey  -and  (h)  after  adjustment 
of  survey  (numbers  in  parenthesis) 


Run 

TEST  1 

TEST  2 

Frame 

4  vs 

(Microns) 

vr 

(Microns) 

Frame 

v* 

(Microns) 

vr 

Microns) 

194 

3.9(  2.2) 

3.2{-4,Q) 

51 

2.8(  0.3) 

8.5(-2.0) 

195 

2.7(  1.4) 

2.1(-5.0) 
7.9(  0.7) 

52 

1.4(-0.7) 

7.6(  1.2) 

W-*E 

196 

1.5(  0.6) 

53 

5.6(  4.0) 

8.5(  2.0) 

197 

-  0.6(-1.0) 

5.5(-1.7) 

54 

“1.3(“2.6) 

6. 2(~0.2) 

198 

0.9(-0 .9) 

0.6(-6.5) 

55 

1.5(  0.6) 

5. 7  (-0.7) 

222 

-  5.4(-0.8) 

-0.8(~2, 1) 

86 

-5.4(  0.8) 

0.5(-1.Z) 

223 

-ll.l(-4.5) 

-1.9(»3.2) 

87 

-6.5(  0.0) 

-  0.5(  0.0) 

N-+S 

224 

-  3.9(~2. 1) 

2. 1(  0.8) 

88 

~8.4(-1.5) 

-  0.5(-2.2) 

225 

-  6.2(  1.0) 

-0.6(-1.9) 

90 

-6.2(  1.4) 

1.6(  0.0) 

226 

-?1.9(-4.3> 

-0.5(-1.8) 

91 

~7.5(  0.6) 

-  0.9(-2.5) 

282 

-  4.2(-2.9) 

**5.6(  0.3) 

118 

-3.4(-2.4) 

-  5.6(  1.4) 

283 

-  4.8(-3.0) 

-7 , 8{- 1 . 9) 

119 

-1.7(-0.3) 

-  9.9(-2.8) 

E  -*W 

284 

-  6.5(-4.4) 

-4.3(  1.6) 

120 

-2.0(-0.2) 

-11.5(-4.4) 

285 

-  2.0(  0.3) 

-6.2(-0.3) 

121 

-1.0(  1.0) 

-  6. 1{  1.1) 

286 

-*  6.4(-3.8) 

-5.2(  0.0) 

122 

-2,5{  0.0) 

-  4.6(  2.5) 

310 

8.2(  0.8) 

-0,2(  0.8) 

153 

8.7(  0.8) 

4. 1(  5.6) 

311 

7.4(  0.4) 

0.2(-0.4) 

154 

7.Q(-0.4) 

1.2(  2.7) 

S  -*N 

312 

6.7(  0.0) 

-5.2(-4.2) 

155 

7.9(  0.9) 

-  4.6(-3. 1) 

313 

6.9(  0.7) 

-4.5(-3,4) 

156 

6.7(  0,1) 

-  G,4(  1.0) 

314 

4 . 8  (—  1  -  4) 

-0.7(  0.4) 

_ 

157 

6.5(-0.3) 

-  2,7(  2.4) 

Grand  RMS  for  point  39  before  adjustment  of  survey:  5,4  microns 
Grand  RMS  for  point  39  after  adjustment  of  survey;  2.2  microns 
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reduction  does  not  consider  survey  error'.,  it  follows  that  their  influence  will  be  reflected, 
os  we  have  seen,  in  the  plate  measuring  residuals.  Inasmuch  as  this  influence  does  indeed 
assume  significance  in  the  residuals  for  several  points,  we  decided  to  adjust  the  survey  to 
mini.nize  5ts  effects.  This  would  most  pioperly  be  done  through  a  rigorous  extension  of 
SMAC  to  consider  errors  in  the  given  control,  os  is  discussed  in  a  later  section.  Since 
such  an  approach  would  be  beyond  the  scope  of  the  present  undertaking,  we  adopted  a 
compromise  along  the  following  lines.  The  projective  parameters  recovered  from  the 
SMAC  reductions  were  considered  to  be  flawless  for  each  frame.  The  coordinates  of  all 
control  points  were  then  regarded  as  unknown  and  the  coordinates  of  a  given  control  point 
were  established  by  spatial  trianguiation  of  rays  from  all  frames  containing  the  point.  As 
many  as  40  rays  were  involved  in  the  trianguiation  of  the  more  central  points,  inasmuch 
as  these  were  recorded  on  c(l  40  frames.  The  corrections  to  the  survey  resulting  from  this 
process  are  listed  in  Table  4  along  with  their  standard  deviations. 

The  nns  values  of  the  corrections  for  those  points  appearing  on  least  20  frames  turn 
out  to  be  0.20  ft.  in  X,  .  17  ft.  in  Y,  and  0.22  ft.  in  Z,  respectively.  Accordingly, 
these  values  were  adopted  as  being  characteristic  of  the  given  sutvey.  The  revised 
corrections  listed  in  the  final  three  columns  of  Table  4  are  derived  from  the  v  sighted 
averages  of  the  original  values  and  the  corresponding  values  resulting  from  the  adjustment. 

Thus  6X,  $Y,  fiZ  are  given  by 

AX 

AY 

AZ 

In  which  0*/  <7y,  Oz  are  the  sigmas  listed  in  Table  4.  By  virtue  of  the  weighting,  the  revised 
coriections  for  those  points  that  are  poorly  determined  in  adjustment  become  close  to  zero,  where¬ 
as  the  revised  corrections  for  those  points  that  are  strongly  determined  become  only  slightly  altered. 
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TABLE  4.  Corrections  to  Survey  of  McLure  Range  Derived  fiom  Aerial  SMAC  Residuals 


Pt. 

No. 

Original  Corrections  to 
Survey  from  Aerotriang . 

Standard  Deviations  of 
Original  Corrections 

Revised  Corrections  from 
Weighted  Averages 

No. 

Frames 

AX 

hL 

0^2 

6X 

6Y 

6Z 

(Ft) 

(Ft) 

(Ft) 

(Ft) 

(Ft) 

(Ft) 

(Ft) 

(Ft) 

(Ft) 

1 

2 

0.66 

0.35 

0.52 

1.40 

1.25 

2.32 

0.01 

0.00 

0.00 

2 

9 

0.33 

-0.15 

0.58 

0.29 

0.23 

0.46 

0.11 

-0.04 

0.11 

3 

14 

0.22 

-0.10 

0.48 

0.21 

0.19 

0.35 

0.10 

-0.04 

0. 14 

4 

20 

-0.25 

0.31 

-0.23 

0.13 

0.15 

0.26 

-0.18 

0.15 

-0.10 

5 

21 

-0.20 

-0.03 

-0. 15 

0.10 

0.14 

0  25 

-0.16 

-0.02 

-0.07 

6 

22 

-0.22 

0.11 

-0.21 

0.08 

0.13 

0.22 

-0.19 

0.06 

-0.11 

7 

22 

-0.17 

-0.C2 

0.05 

0.06 

0,13 

0.22 

-0.15 

-0.01 

0.03 

8 

21 

-0.19 

-0.07 

-0.07 

0.05 

0.13 

0.22 

-0.18 

-0.04 

-0.04 

9 

23 

-0.07 

-0.13 

0.18 

0.05 

0.12 

0.21 

-0.06 

-0.08 

0.10 

10 

21 

-0. 17 

-0.35 

0.56 

0.08 

0.14 

0.24 

-0. 15 

-0.19 

0.26 

il 

10 

-0.54 

-0.67 

0.91 

0.25 

0.25 

0.45 

-0.21 

-0.18 

0.18 

12 

7 

0.22 

0.37 

-0.94 

0.37 

0.35 

0.62 

0.05 

0.06 

-0.11 

13 

2 

0.92 

1.04 

-1.40 

1.03 

1.07 

2.00 

0.03 

0.02 

-0.02 

14 

2 

0.02 

-0.01 

-0.28 

1.24 

0.25 

2.00 

0.00 

0.00 

0.00 

15 

10 

0.41 

0.02 

-0.91 

0.21 

0.11 

0.34 

0.20 

0.01 

-0.28 

16 

34 

0.02 

-0.04 

0.02 

0.08 

0.06 

0.15 

0.02 

-0.03 

0.01 

17 

27 

-0.17 

-0.08 

0.26 

0.09 

0.09 

0.18 

-0.14 

-0.06 

0.16 

18 

40 

0.03 

0.04 

-0.11 

0.04 

0.05 

0.12 

0.03 

0.04 

-0.08 

19 

40 

0.04 

0.15 

0.11 

0.04 

0.05 

0.12 

0.04 

0.14 

0.08 

20 

40 

0.14 

0.15 

0.16 

0.04 

0.05 

0.12 

0.13 

0. 14 

0. 12 

21 

40 

-0.08 

0.01 

0.01 

0.05 

0.05 

0.12 

-0.08 

0.01 

0.01 

22 

30 

0.09 

-0.08 

-0.29 

0.07 

0.08 

0.17 

0.08 

-0.06 

-0.18 

23 

39 

0.05 

-0.05 

-0.30 

0.06 

0.05 

0.13 

0.05 

-0.05 

-0.22 

24 

38 

0.12 

-0.04 

-0.21 

0.07 

0.06 

0.13 

0.11 

-0.03 

-0.16 

25 

14 

-0.03 

0.45 

-0.09 

0.18 

0.12 

0.28 

-0.02 

0.27 

-0.04 
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TABLE  4  (Continued) 


Original  Corrections  to 

Standard  Deviations  of 

Revised  Corrections  from 

Pt. 

No. 

Survey  from  Aerotriang . 

Original  Corrections 

Weighted  Averages 

No. 

Frames 

AX 

ay 

&  i 

QAy 

z 

6X 

6Y 

6Z 

(Ft) 

(Ft) 

(Ft)  1 

(Ft) 

(Ft) 

(Ft) 

(Ft) 

(Ft) 

26 

21 

0.38 

0.09 

0.20 

0.12 

0.08 

0.  i8 

0.28 

0.07 

0.12 

27 

5 

0.36 

-0.16 

-0.02 

0.27 

0.17 

0.45 

0. 13 

-0.07 

0.00 

28 

4 

0.42 

-0.36 

0.33 

0.34 

0.22 

0.54 

0.11 

-0.11 

0.05 

29 

20 

-0.18 

0.18 

-0.09 

0.13 

0.06 

0.19 

-0.13 

0.16 

-0.05 

30 

7 

0.15 

-0.32 

0.33 

0.24 

0.09 

0.40 

0.06 

-0.24 

0.08 

31 

21 

-0.39 

-0. 15 

-0.22 

0.12 

0.06 

0.18 

-0.29 

-0.13 

-0.13 

32 

27 

-0.32 

-0.25 

-0.30 

0.09 

0.05 

0.16 

-0.27 

-0.22 

-0.20 

33 

36 

-0.06 

-0.17 

-0. 12 

0.07 

0.04 

0.13 

-0.05 

-0.16 

-0.09 

34 

20 

-0.00 

-0.06 

-0.44 

0.07 

0.07 

0.20 

-0.00 

-0.05 

-0.25 

35 

38 

-0.03 

-0.14 

-0.10 

0.05 

0.04 

0.13 

-0.03 

-0.13 

-0.07 

36 

39 

-0.C9 

0.10 

-0.20 

0.05 

0.04 

0.12 

-0.08 

0.09 

0.18 

37 

40 

0.01 

-0.27 

-0.12 

0.05 

0.04 

0.13 

0.01 

-0.25 

0.16 

33  ! 

40 

-0.09 

-0.24 

0.00 

0.08 

0.04 

0.15 

-0.09 

-0.22 

-0.08 

39  | 

40 

-0.08 

0.52 

-0.09 

0.08 

0.05 

0.15 

-0.08 

0.48 

-0.16 

40 

39 

0.21 

0.08 

0.33 

0.16 

0.07 

0.24 

0.20 

0.07 

-0.09 

41 

33 

0.13 

0.01 

0.00 

0.08 

0.04 

0.15 

0.11 

0.01 

0.00 

42 

33 

0.23 

-0.13 

-0,09 

0.08 

0.05 

0.15 

0.20 

-0.12 

-0.06 

43 

16 

0.10 

-0.05 

0.33 

0.16 

0.07 

0.24 

0.06 

-0.04 

0.15 

44 

1 

- 

- 

- 

- 

- 

0.00 

0.00 

0.00 

45 

22 

-0.38 

0.14 

0.54 

0.11 

0.12 

0.21 

-0.29 

0.08 

0.29 

46 

30 

-0.13 

0.23 

0.09 

0.06 

0.09 

0.16 

-0.12 

0.17 

0,06 

47 

33 

0.12 

0.14 

-0.04 

0.04 

0.08 

0.14 

0.12 

0.11 

-0.03 

48 

28 

-0.22 

-0.C4 

-0.01 

0.05 

0.10 

0.17 

-0.21 

-0.03 

0.00 

49 

30 

-0.14 

-0.21 

-0.06 

0.06 

0.09 

0. 16 

-0.13 

-0.15 

-0.04 

50 

16 

-0.06 

0.32 

-0.20 

0.16 

O.lo 

0.26 

0.04 

0.18 

-0.08 

51 

25 

0.02 

-0.02 

-0.34 

0.05 

0.10 

0.20 

0.02 

-0.01 

-0.19 

52 

20 

0.06 

-0.18 

-0.27 

0.05 

0.07 

0.20 

0.06 

-0.15 

-0.15 

-^-r,-,„  '  ^  -rr“-* 

f-  -  -■ 
f  --: 

I 

* 

I  j 
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It  will  be  noted  from  Table  4  that  the  stronger  triangulations  produced  proportional 
accuracies  on  the  order  of  1  part  in  300,000  of  the  flying  height  in  the  horizontal 
coordinates  X  and  Y  and  1  part  in  100,000  of  the  fly'tjg  height  in  the  vertical  coordinate 
Z.  Such  a  degree  of  accuracy  is  unprecidented  in  aerotriangulation,  but  then  so  also  is 
the  experiment  itself.  The  results  indicate  that  the  unconventional  ‘cloyerleaf1  scheme 
of  aerotriangulation  can  maintain  first  order  accuracies  in  the  process  of  densification  of  a 
sufficiently  dense  first  order  network. 

When  the  revised  survey  was  exercised  ?n  a  repetition  of  the  two  SMAC  calibrations, 
the  resulting  inferior  projective  parameters  were  not  sensibly  altered  from  their  former  values. 

On  the  other  bond,  the  combined  set  of  measuring  residuals  from  both  tests  was  moderately 
improved,  the  rms  error  being  reduced  from  3.6  microns  to  2.9  microns.  As  can  be  seen 
from  the  plotted  results  in  Figure  6,  the  improvement  is  especially  great  over  the  central 
two  thirds  of  the  plate;  here,  the  residuals  are  almost  perfectly  random.  Nonetheless,  . 

systematic  tendencies  do  persist  near  the  edges  of  the  format.  Indeed,  now  that  the  effects 
of  survey  errors  are  diminished,  persistent  systematic  components  are  somewhat  more  clearly 
defined  than  they  were  previously  (Figure  5), 


As  an  example  of  the  improvement  in  the  residuals  wrought  by  the  adjusted  survey, 
v/e  have  also  provided  in  Table  3  the  plate  residuals  for  Point  No.  3?  following  the 
application  of  survey  corrections  in  the  revised  SMAC  reductions.  The  rms  value  of  the 
original  SMAC  residuals  for  Point  No.  39  is  5.4  microns;  after  the  application  of  survey 
corrections  the  rms  value  is  reduced  to  2.2  microns. 


3.2.6  Empirical  Modeling  of  Residual  Systematic  Errors 

While  we  were  unable  to  isolate  the  causes  of  the  observed  residual  systematic  error, 
we  did  find  that  systematic  components  Ax,  Ay  of  the  residual  vectors  could  be  adequately 
represented  by  general  fifth  order  polynomials  of  the  form; 
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figure  6.  Revised  version  of  Figure  5  offer  adjustment  of  Me  Lure  survey  (rms  value  of 
residuals  =2*9  microns). 


The  coefficients  of  Ofjj  ,  /3jj  of  these  error  functions  were  establi*hed  by  means  of  a  least 
squares  fit  to  the  residuals  ?n  Figure  6.  The  rms  error  of  the  revised  residuals  resulting  from  the 
fit  turned  out  to  be  2.3  microns,  a  substantia!  improvement  over  the  2.9  microns  prior  to  the  fit. 
Contour  maps  of  the  fitted  error  functions  for  Ax  and  Ay  are  shown  in  Figures  7  end  8.  Ideally, 
the  error  functions  could  be  regarded  as  empirically  established  corrections  accounting  for  those 
systematic  errors  that  are  inadequately  modeled  in  the  SMAC  calibration.  However,  before  they 
can  be  applied  with  confidence  as  standard  corrections,  their  general  validity  must  be  established 
in  a  variety  of  situations.  It  remains  to  be  determined,  for  example,  to  what  extent  the  error 
functions  remain  valid  when  such  factors  as  the  following  are  altered: 

(a)  flying  altitude  and  velocity; 

(b)  calibrated  coordinates  of  platen  reseau; 

(c)  comparator  used  for  measuring; 

(d)  exercise  of  image  motion  compensation; 

(e)  film  magazine  (and  hence  platen)  used; 

(f)  photo  processing  of  film. 

Until  such  time  as  appropriate  supplementary  testing  can  be  performed,  the  error  functions  deriyed 
from  the  SMAC  residuals  must  be  regarded  as  having  no  more  than  provisioned  validity.  The  fact 
that  their  application  does  reduce  the  rms  of  the  revised  residuals  to  2.3  microns  in  the  present 
instance  is  clearly  a  strong  inducement  for  a  thorough  investigation  of  the  soundness  of  such 
empirically  derived  error  functions.  Thus,  in  conjunction  with  possible  future  testing  programs, 
specific  attention  should  be  directed  toward  deliberate  variation  of  those  factors  that  could  con¬ 
ceivably  influence  the  results.  By  so  doing,  one  would  be  better  able  to  isolate  specific  causes 
of  residual  systematic  error. 

It  should  be  appreciated  that  when  only  a  modest  number  of  control  points  appear  on  the 
typical  frame,  empirical  modeling  of  residual  systematic  errors  becomes  practicable  solely  by 
virtue  of  the  large  sample  of  residuals  that  can  nonetheless  be  generated  by  a  SMAC  reduction 
of  a  moderately  large  number  of  frames. 


figure  7.  feocontours  of  empirical  corrections  to  x  coordinates  of  frames  exposed  by  Camera  006 
(contour  interval  =  1  micron). 


While  we  believe  that  the  ultimate  refinement  of  the  process  of  camera  calibration 
may  well  consist  of  the  empirical  modeling  of  the  residual  systematic  errors  remaining  after 
the  execution  of  a  comprehensive  SMAC  reduction,  this  is  hardly  a  matter  for  complacence. 

As  long  as  significant  residual  systematic  errors  do  persist,  their  source  should  be  a  matter  of 
continuing  concern,  and  one  should  not  remain  content  until  satisfactory  physical  explanations 
of  such  errors  are  established.  Physical  explanations  can  be  expected  to  lead  either  to  imple¬ 
mentation  of  physical  corrections,  which  is  desirable  when  practicable,  or  to  mathematical 
modeling  of  the  physical  process,  which  is  far  to  be  preferred  over  blind  empirical  modeling. 

At  best,  empirical  modeling  should  be  viewed  as  a  crutch  for  ones  ignorance,  rather  than  as  a 
panacea. 

3.2.7  Empirical  Determination  of  Weighting  Functions  of  Plate  Coordinates 

# 

It  Is  generally  appreciated  that  effective  accuracies  of  plate  coordinates  vary  with  radial 
distance-  Only  a  minor  part  of  such  variation  can,  in  general,  be  attributed  to  shortcomings  of 
the  plate  measuring  comparator.  On  the  other  hand,  much  of  it  can  presumably  be  attributed  to 
the  variations  in  image  quality  resulting  from  decreasing  resolution  of  the  lens  with  increasing 
radial  distance.  Other  major  contributions  might  well  stem  from  departures  from  flatness  of  the 
photographic  emulsion  at  the  instant  of  exposure  end  from  uncompensated  deformation  of  the  film. 

Moren  (1965)  has  reported  an  analysis  of  the  pooled  residuals  obtained  from  a  series  of 
camera  calibrations  performed  from  photography  of  the  Oland  test  field.  The  calibrations  were 
accomplished  by  the  grid  method  developed  by  Hallert  (1954).  Moren  found  that  the  dependence 
of  the  standard  deviation  of  image  coordinates  on  radial  distance  r  is  described  by  the  formula; 

(95)  a  -•=  2.-1  +  .  053  r  +  .  00023  r2 

in  which  <j  is  in  microns  and  r  is  In  millimeters  Thus  a  varies  from  about  2  microns  in  the  center 
of  the  frame  to  about  10  microns  at  r  =  100  mm,  and  to  about  15  microns  in  the  corners. 
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We  considered  it  pertinent  to  our  Investigation  to  follow  Moren‘s  example  by  performing 
a  similar  analysis  of  the  pooled  SMAC  residuals.  Unlike  Moreri,  we  decided  to  perforin  separate 
analyses  of  the  redial  and  tangential  components  of  the  residual  vectors.  Accordingly,  we  divided 
the  format  out  to  r  -  120  mm  Into  five  circular  zones  of’  approximately  equal  area  arid  then  com¬ 
puted  the  rms  errors  of  the  radial  and  tangential  components  of  the  residuals  in  Figure  6  for  each 
of  the  resulting  zones.  We  also  computed  the  coefficient  of  correlation  of  the  radial  and  tangential 
components  for  each  zone.  The  resulting  values  are  listed  in  Table  5  and  are  plotted  in  Figure  9. 


Table  5.  Statistical  Analysis  of  Residuals  Plotted  in  Figure  6. 


Interval 

Mean  Radius 

r 

No.  Residual 
Vectors 

RMS  of  Radial 
Components:  Or 

RMS  of  Tangential 
Components:  o> 

Correlation 

Coeff: 

0  -  53mm 

26.5mm 

269 

2.63p 

2.25p 

+.012 

53-  75 

64.0 

251 

3.10 

2.09 

-.035 

75-  92 

83.5 

217 

3.48 

2.10 

-.054 

92  -  106 

99.0 

208 

3-60 

2.08 

-.117 

106-120 

113.0 

177 

3.89 

2.66 

-.032 

5120 

131.2 

88 

4.23 

3.13 

-.174 

*None  of  the  Pj+  differs  significantly  from  zero  at  the  95%  level  of  confidence. 


interestingly  enough,  we  find  that  the  sigmas  of  the  tangential  components  actually  decrease 
slightly  with  radial  distance  our  to  r  =  100mm.  Moreover,  they  are  significantly  smaller  throughout 
the  format  than  the  sigmas  of  the  radial  components.  This  suggests  that  departures  from  flatness  of 
the  surface  of  the  emulsion  may  be  of  greater  importance  to  our  results  than  are  variations  in  photo¬ 
graphic  resolution.  We  find  the  variation  in  the  standard  deviations  of  the  image  coordinates  to  be 
much  less  pronounced  throughout  the  format  than  that  experienced  by  Moren.  This  may  be  an  indi¬ 
cation  of  superior  optical  performance  by  the  Geocon  IV  Lens.  However,  we  suspect  that  Moren‘s 
results  may  reflect,  to  seme  extent,  the  influence  of  decentering  distortion  inasmuch  as  th»s  was 
not  explicitly  accounted  for  in  hit  calibrations. 
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!t  wM(  be  noted  that  we  chose  fo  employ  even  powered  polynomials  to  represent  crr 
and  $  In  Figuie  9.  This  is  strictly  a  matter  of  preference  stemming  from  intuitive  optica! 
considerations.  Two  terms  of  the  expansion  were  found  adequate  to  fit  the  observed  values 
of  0r  whereas  thne®  were  needed  for  . 

Had  *ve  used  the  residuals  remaining  after  the  fir  of  the  empirical  functions  £y 
discussed  in  the  preceding  section,  f,he  resulting  values  of  gf  and  g^  would  have  exhibited 
even  less  variation  than  that  encountered  above.  Indeed,  should  the  empirical  corrections 
prove  to  have  general  validity  for  the  camera  In  question  (KC-6A  No.  006),  one  could 
safely  regard  the  stanclard  deviations  of  the  image  coordinates  as  being  essentially  invariant 
throughout  the  format.  However,  pending  clarification  of  this  matter,  vve  consider  it  more 
prudent  to  adopt  the  results  of  Table  5  and  Figure  9. 

The  practical  application  of  the  foregoing  results  to  analytical  cerotriangulation 
performed  from  photography  taken  by  Camera  006  remains  to  be  considered.  Even  though 
Table  5  indicates  that  radial  and  tangential  components  are  essentially  uncorrelated,  this  is 
not  necessarily  true  of  the  x,y  components.  Inasmuch  *  components  of  error  are  related 
to  radial  and  tangential  components  by 


y3  crfs  +  *3  ffi8  -~y(ars  ~  at?  ) 
-xy(ar3  ~  at3)  x2  <7r3  +  y2  fft3 


For  the  special  case  in  which  ar  =  at  =  o,  one  has  Oy  =  oy  =  a  and  ax y  -  0 .  In  the  general 
case  in  which  ar  /  a^,  axy  =  0  only  when  x  or  y  is  equal  to  zero.  When  x  =  y  the  coefficients 
of  correlation  between  x  and  y  is  given  by  fay  ~  {of  -  of)/  {of  +  e^). 

Several  of  the  values  of  ar  and  G>  in  Figure  9  lead  to  coefficients  of  correlation  approaching  0.5, 
In  view  of  this,  it  is  clear  that  the  adjustment  of  photogrammetric  observations  in  the  process  of 
analytical  aerotriangulation  should  admit  the  possibility  of  correlated  x,y  coordinates,  as  is  done 
in  our  general  theory  (Brown,  1958;  Brown,  Davis,  Johnson,  1964).  The  covariance  matrices 
required  for  the  adjustment  can  be  computed  by  means  of  equation  (97)  in  which  ffr  and  0^  are 
defined  by  the  polynomial  expressions  given  in  Figure  9. 

The  weighting  functions  generated  from  an  cnalysis  of  SMAC  residuals  constitutes  a  most 
important  by-product  of  the  calibration  of  a  given  camera.  We  believe  that  in  the  future,  the 
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determination  of  such  weighting  functions  should  be  required  as  one  of  the  standard  outputs  of 
the  process  of  calibration. 

3.2-8  Stellar  SMAC  Calibration  of  KC-6A  No.  006  and  No.  008. 

One  of  the  objectives  of  the  testing  program  was  to  compare  the  results  obtained  from  an 
operational  aerial  calibration  with  those  obtained  from  a  stellar  calibration.  Because  of  delays 
encountered  at  Wright-Fatterson  Air  Force  Base  in  the  design  and  fabrication  of  a  suitable  plate 
back  for  the  camera,  the  stellar  shots  were  nor  made  until  August  1967,  or  about  five  months 
offer  the  completion  of  the  aerial  tests. 

For  the  stellar  photographs  it  was  necessary  to  employ  12  by  12  inch  plates  in  order  to 
accommodate  the  three  pads  that  define  the  focal  plane  of  the  camera.  Kodak  microflat  plates 
(6  mm  thick)  coated  with  1 03— F  Spectroscopic  emulsion  were  used  for  the  photographic  recording. 

In  conventior-i  stellar  calibrations  of  a  camera,  great  pains  must  be  taken  to  insure 
absolute  Stability  of  the  camera  throughout  the  period  of  recording.  In  addition,  each  exposure 
must  be  precisely  timed  in  order  to  account  for  the  rotation  of  the  earth  between  exposures. 


With  cameras  having  the  fairly  small  aperture  typical  of  aerial  mapping  camera$/  one  usually 
has  to  make  periodic  exposures  over  a  period  of  an  hour  or  so  in  order  to  obtain  the  number 
(200-300)  and  distribution  of  stellar  images  needed  for  a  comprehensive  calibration.  Unfortunately, 
ones  success  in  maintaining  stability  and  precise  timing  is  not  ordinarily  known  until  after  the 
reduction  has  been  completed.  When  inconsistencies  are  then  found  to  result  from  deficiencies 
in  stability  or  timing,  the  test  will  usually  have  to  be  repeated.  Because  of  such  demanding  re¬ 
quirements,  conventional  stellar  calibration  of  mapping  cameras  has  been  too  expensive  for  general 
application.  Although  expense  was  not  a  serious  factor  in  the  calibration  of  the  KC-6A  cameras, 
their  limited  availability  made  ihe  execution  of  a  conventional  stellar  calibration  too  risky  to  con¬ 
sider  in  view  of  the  alternative  afforded  by  a  stellar  SMAC  calibration.  Inasmuch  as  each  stellar 
exposure  can  be  regarded  as  a  separate  frame,  a  stellar  SMAC  calibration  imposes  no  requirements 
for  stability.  Indeed,  if  desired,  the  camera  could  be  deliberately  reoriented  between  exposures. 
Likewise,  stellar  SMAC  imposes  no  requirements  for  the  precise  timing  of  stellar  exposures.  This'is 
because  an  error  in  the  assigned  time  of  an  exposure  is  projectively  equivalent  to  an  error  in  the 
instantaneous  right  ascension  of  the  camera  axis*  Accordingly,  the  angular  elements  of  orientation 
for  a  given  frame  can  fully  accommodate  timing  errors.  Indeed,  were  it  not  necessary  to  correct 
for  the  effects  of  atmospheric  refraction,  one  could  forego  all  need  for  timing  of  exposures.  As  it 
is,  one  must  know  the  approximate  zenith  distance  of  each  star  at  the  instant  of  exposure  in  order 
to  correct  for  refraction.  Since  an  accuracy  of  0?5  is  altogether  adequate  for  this  purpose,  it  is 
more  than  sufficient  if  times  of  exposures  are  known  to  an  accuracy  of  one  minute.  Accordingly, 
a  good  wristwatch  is  an  adequate  chronometer  for  a  stellar  SMAC  calibration.  Because  timing 
requirements  are  so  coarse,  the  special,  electrically  timed  foreshutter  needed  in  a  conventional 
stellar  calibration  can  be  replaced  in  a  SMAC  calibration  with  a  hand-operated  piece  of  darkened 
cardboard.  Because  stability  is  of  no  concern,  a  SMAC  calibration  demands  nothing  in  the  way  of 
special  facilities  such  as  massive  concrete  pedestals  enclosed  in  observing  domes.  Thus,  a  stellar 
SMAC  calibration  not  only  has  the  advantage  of  rendering  the  results  totally  immune  to  deficiencies 
in  stability  and  timing,  but  it  also  vastly  expedites  and  simplifies  the  field  operation. 


Zenithal  stellar  plates  were  obtained  from  cameras  006  and  008.  (Camera  005  was 
unavailable  at  the  time.)  Stellar  exposures  were  made  at  ten  minute  Intervals  over  a  period 
of  50  minutes.  About  80  star  trails  were  selected  for  measurement  on  the  plate  from  C06  and 
about  70  from  008.  Up  to  6  Images  were  measured  on  each  selected  trail,  leading  to  totals 
of  436  measured  images  for  006  and  385  measured  images  for  008.  The  calibrated  values  of  the 
interior  projective  parameters  resulting  from  the  various  SMAC  reductions  are  given  in  Table  6 
and  the  calibrated  distortion  functions  are  plotted  in  Figure  10.  For  future  reference.  Table  6 
also  includes  corresponding  results  from  a  laboratory  calibration  performed  by  Fairchild. 

In  comparing  the  results  of  the  stellar  calibration  for  006  with  the  results  of  the  aerial 
calibration  we  find: 

(1)  the  calibrated  principal  distance  from  stellar  SMAC  is  about  30  microns  shorter 

than  that  from  aerial  SMAC;  • 

(2)  x  and  y  from  stellar  SMAC  are  much  smaller  than  the  values  from  aerial  SMAC; 

P  P 

(3)  the  radial  distortion  curve  from  stellar  SMAC  is  in  excellent  agreement  with  the 
mean  of  the  two  aerial  SMAC  curves; 

(4)  the  profile  function  of  decentering  distortion  is  not  in  good  agreement  with  the 
corresponding  profile  functions  from  aerial  SMAC. 

We  believe  that  (1)  can  be  considered  a  measure  of  the  actual  physical  change  resulting  from 
the  orientation  of  the  camera  which  was  pointed  downward  for  the  aerial  SMAC  exposure  and 
upward  for  the  stellar  SMAC  exposures.  The  sign  of  the  difference  between  the  principal  distances 
agrees  with  what  one  would  expect  from  the  influence  of  gravity  on  the  position  of  the  lens  relative 
to  the  focal  plane. 

The  fact  that  x  ,  y  from  stellar  SMAC  fails  to  confirm  the  large  values  resulting  from 

li 

aerial  SMAC  lends  confirmation  to  the  hypothesis  that  the  latter  are  compromised  (most  likely)  by 
errors  of  synchronization . 

While  the  excellent  agreement  between  the  radial  distortion  curves  is  assuring,  this  is  offset 
by  the  rather  poor  agreement  between  decentering  profiles  which  is  too  great  to  be  accidental. 

WJe  suspect  that  an  actual  change  in  decentering  of  the  lens  elements  may  well  have  accompanied 
the  opposing  changes  in  orientation.  If  this  should  indeed  prove  to  be  the  case,  the  validity  of 
the  stellar  method  for  the  calibration  of  mapping  cameras  would  be  significantly  compromised. 
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TABLE  6.  Elements  of  Interior  Orientation  of  KC-6A  Cameras  Resulting  from 
Different  Calibrations 


Type  of 

Parameter 

camera 

Calibration 

005 

006 

008 

Aerial  SMAC 

xp  (mm) 

.167+. 004 

.231*. 003 

— 

(Test  No.  1) 

yp  (mm) 

-.045+. 004 

.  104+.003 

- 

c  (mm) 

152.539i.002 

151. 23 1+. 002 

- 

Aerial  SMAC 

xp  (mm) 

- 

.222+. 004 

.175+, 003 

(Test  No.  2) 

"p  (mm) 

- 

. 09?± . 004 

-.052+. 003 

c  (mm) 

- 

151. 227+.  002 

154. 626+. 002 

Stellar  SMAC 

Xp  (mm) 

- 

-.035+.  002 

-.033+. 003 

yv  (mm) 

_ 

-.017+. 002 

.00Q+.003 

c  (mm) 

- 

151.200+.001 

154.582i.002 

*  Laboratory 

xp  (mm) 

-.011 

.015 

-.009 

(Fairchild 

yp  (mm) 

.020 

-.013 

.014 

multicollimator) 

c  (mm) 

152.565 

151.220 

154.626 

*Standard  deviations  of  laboratory  results  are  not  given  in  calibration  certificate. 


10u  T  6.  RADIAL  DISTORTION,  CAMERA  006  (c  =  151 .234mm) 
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Figure  10.  Radial  distortion  curves  and  decentering  profiles  resulting  from  Stellar  SMAC 
calibrations  of  Cameras  006  and  008  (accompanied  by  one  sigina  error  bounds), 
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A  possible  way  around  this  difficulty  would  be  to  use  an  optically  flat  mirror  or  large  pool  of 
mercury  to  permit  exposures  about  the  zenith  to  he  made  with  the  camera  pointed  downwards. 

We  would  recommend  that  a  future  investigation  be  undertaken  along  such  linos  in  order  to 
establish  the  precise  influence  of  the  orientation  of  the  camera  on  the  interior  projective  para- 
meters. 

The  agreement  between  the  stellar  SMAC  calibration  of  radial  distortion  for  008 
(Figure  10)  and  the  aerial  SMAC-  calibration  (Appendix  B)  is  not  as  good  as  that  between  the 
corresponding  calibrations  for  006.  This  could  well  be  attributable  to  the  effects  of  curvature 
of  the  surface  of  the  photographic  plate,  a  topic  taken  up  below. 

3.2.9  Analyses  of  Stellar  S MAC  Residuals 

The  residual  vectors  resulting  from  the  stellar  SMAC  calibrations  of  cameras  006  and  008 
are  plotted  in  Figures  Hand  12.  A  grid  having  special  vectors  associated  with  each  intersection 
has  been  superimposed  on  each  figure.  These  are  for  future  consideration  and  are  to  be  ignored 

• 

for  the  time  being,  in  both  figures,  those  images  on  common  star  trails  are  joined  by  straight 
lines. 

Even  a  casual  examination  of  the  residual  vectors  in  Figures  11  and  12  reveals  the  presence 
of  pronounced  systematic  tendencies  (especially  so  in  Figure  12).  The  rms  errors  of  the  residual 
vectors  of  3.5  and  4.4  microns  respectively  are  unduly  large  for  a  stellar  calibration.  When  the 
residual  vectors  are  analyzed  in  terms  of  radial  and  tangential  components,  we  find  the  radial 
components  to  have  rms  values  of  3-0  and  4.0  microns,  respectively,  and  the  tangential  components 
to  have  rms  values  of  1 .8  and  1 .9  microns.  This  immediately  suggests  that  departures  of  the  photo¬ 
graphic  surface  from  a  true  plane  might  be  at  the  root  of  the  difficulty.  However,  the  specifica¬ 
tions  for  Kodak  microflat  glass  plates  states:  "the  front  surface  of  microflat  glass,  on  which  the 
emulsion  is  coated,  does  not  depart  by  more  than  0.00002  inch  per  linear  inch  from  a  true  plane." 
This  would  permit  a  departure  from  flatness  of,  at  most,  about  +  3  microns  across  the  diagonal 
of  the  9x9  inch  format,  which  is  not  nearly  enough  to  account  for  the  observed  effect.  In  looking 
more  deeply  into  the  matter,  it  occurred  to  us  that  the  above  specification  actually  applies  to 
uncoaied  plates.  This  consideration  brought  to  our  recollection  an  article  in  the  Astronomical 
Journal  on  the  subject  of  the  effects  on  flatness  of  the  photographic  emulsion.  Upon  searching,  we 


Figure  12,  Residual  vectors  from  Stellar  SMAC  calibration  of  Comera  008  (rms  value  of 
residuals  =  4.6  microns). 
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found  the  desired  article.  It  was  in  the  1942  Volume  of  the  Journal  and  was  entitled,  "The 
Concavity  of  Photographic  Pictes",  by  Professor  Frank  Schiesinger  of  Yale  University  Observa¬ 
tory.  Because  of  its  pertinence  and  brevity,  v/e  quote  it  Here  in  full: 


"The  writer  has  several  times  itated  that  the  shrinking  of  the  gelatins  in  drying  caused 
photographic  plates  to  be  concave  on  the  coated  side.  Doubts  of  the  correctness  of  this 
statement  having  been  expressed,  the  matter  now  has  been  definitely  tested  by  measuring 
the  curvature  of  the  identical  plates  before  and  after  they  had  been  coated.  This  was 
dene  by  means  of  the  simple  spheromete-r  described  In  the  introduction  to  Volume  9  of  the 
Yale  Observatory  Transactions.  The  plates,  76  in  number,  are  of  plate  glass:  they  are 
43cm  square. by  5  or  6  mm  thick.  In  every  case  they  were  found  to  be  more  concave  on 
the  coated  side  than  they  were  before  coating.  The  change  varied  from  0.004  mm  to 
Q.HG  mm,  the  overage  being  0.046mm.  " 

"This  concavity  would  give  rise  to  an  appreciable  error  in  star  positions  near  the  edges  of 
the  plates  if  we  did  nor  flatten  out  the  plates  when  they  were  being  exposed. " 


"While  this  experiment  proves  that  the  shrinking  o*  the  gelatine  on  plate-glass  causes  their 
concavity,  the  conclusion  does  not  necessarily  apply  to  ordinary  glass.  However,  by 
measuring  14  such  plates  in  the  16.5  by  2.15  cm  size,  about  1.6mm  thick,  we  found  that 
they  were  concave  on  the  coated  side  by  an  average  of  0.023 mm.  Allowing  for  the  differ¬ 
ence  In  size  and  thickness  between  the  plate-glass  and  these  plates,  it  looks  as  though  the 
concavity  In  the  latter  case  were  likewise  due  to  the  contraction  of  the  gelatine." 


"1  am  indebted  to  Mr.  HOLLANDER  for  making  the  necessary  measurement. " 


In  view  of  Professor  Schlesinger’s  findings,  we  deemed  it  advisable  to  investigate  the  matter 
of  the  flatness  of  the  actual  plates  used  in  the  stellar  SMAC  reductions.  On  the  glass  side  of  each 
plate,  we  established  a  grid  array  of  25  points  evenly  spaced  over  the  9x9  inch  area  corresponding 
to  the  photographic  format.  The  plate  was  then  placed,  glass  side  a‘own,  on  a  Brown  and  Sharpe 
24"  x  30"  granite  surfece  plate  accurate  to  ±  0.6  microns.  A  Brown  and  Sharpe  dial  indicator 
was  employed  to  measure  relative  spot  elevations  of  the  points  on  the  emufsion  coiresponding  to 
the  grid.  The  nvr'  ;ienf$  were  recorded  to  the  nearest  micron  (or,  more  precisely,  to  the 
nearest  .00005  inen),  with  double  settings  being  made  on  each  point,  but  not  In  succession. 
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The  rms  of  the  means  of  the  settings  was  found  to  be  1 .5  microns.  A  best  fitting  least  squares 
plane  was  established  from  the  measurements.  The  departures  of  the  spot  elevations  from  Hie 
best  fitting  plane  were  found  to  be  entirely  regular  and  were  used  to  generate  horizontal  and 
vertical  profiles  which,  in  turn,  were  used  to  generate  contour  maps.  The  contour  maps  of  the 
Surfaces  of  the  hvo  plates  relative  to  best  fitting  planes  are  presented  in  Figures  13a  and  13b. 

We  note  from  the  contour  maps  that  ihe  surfaces  are  not  concave,  os  reported  by 
Schiesingcr,  but  rather  are  roughly  saddle  shaped,  in  addition,  we  see  that  the  gradients  of  the 
surfaces  in  certain  areas  become  up  to  fen  times  greafei  than  that  allowed  by  the  specifications  for 
micro  flat  glass.  From  this  we  may  conclude  that  either  the  plates  are  significantly  deformed  by 
the  stresses  induced  by  the  drying  of  the  photographic  emulsion,  or  else  the  thickness  of  the 
emulsion  itself  is  far  from  uniform.  To  check  the  latter,  we  carefully  scaped  away  a  small  patch 
of  emulsion  in  the  vicinity  of  each  grid  point  and  then  took  supplementary  spot  readings  to  ascertain 
the  thickness  of  the  emulsion  We  found  the  variation  in  thickness  to  be  negligible,  thus  confirming 
that  it  is  ihe  piafe  itself  the  f-  deformed. 

Going  a  step  further,  we  derived  spot  elevations  Az  from  the  contour  maps  for  81  points 
evenly  spaced  over  the  format.  For  each  point  we  then  computed 

(98)  6r  =  -  Az 
c 

in  which  r  Is  the  radial  distance  to  the  point  and  c  is  the  focal  length  of  the  camera.  The 
quantity  6r  defines  the  contribution  to  the  radial  residua i  to  be  expected  from  departure  of 
the  photographic  surface  from  the  best  fitting  plane.  The  resulting  radial  vectors  at  the  grid 
intersections  are  superimposed  on  the  residua!  plots  of  Figures  11  and  12.  In  Figure  11  the 
correlation  between  the  radial  vectors  on  the  grid  and  the  systematic  components  of  the  SMAC 
residuals  is  seen  to  be  remarkably  good  and,  in  our  view,  is  clearly  too  strong  to  be  accidental, 
the  correlation  between  the  two  sets  of  vectors  in  Figure  12,  though  not  quite  as  pronounced  as 
in  Figure  11,  is  nonetheless  strong. 
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In  both  cases,  but  particularly  in  Figure  12,  the  correlation  in  direction  is  more  pronounced 
than  the  correlation  in  magnitude.  For  the  most  part,  the  magnitudes  of  the  grid  residuals 
tend  to  be  somewhat  smaller  than  those  of  their  counterparts  from  the  SMAC  reductions.  This 
suggests  that  the  deformation  of  the  plate  following  photographic  processing  may  actually  be 
somewhat  less  severe  than  that  of  the  undeveloped  plate  (quite  possibly,  stresses  in  the  emulsion 
become  partially  relieved  as  a  result  of  processing).  We  recommend  that  this  matter  be  taken 
up  in  future  investigations.  With  the  aid  of  air  gauging  techniques  it  should  be  possible  to  work 
out  a  practical  scheme  for  measuring  the  surface  of  undeveloped  plates  in  total  darkness. 


Even  though  modern  emulsions  are  likely  to  be  appreciably  thinner  than. those  employed 
in  Schlesinger's  day,  our  results  suggest  that  his  findings  continue  io  have  pertinence.  Even 
with  6  mm  thick  plates  ,  one  is  not  assured  that  the  flatness  prior  to  coating  wili  be  maintained 

* 

after  coating.  While  the  departures  from  flatness  that  we  have  observed  are  admittedly  small 
by  customary  standards  (rms  values  of  6.2  and  7.8  microns  for  <!. ,  two  plates  measured),  they 
are,  as  we  have  seen,  enough  to  contaminate  our  results  to  a  significant  degree.  A  brute  force 
solution  to  the  problem  would  lie  in  the  use  of  extra  thick  plates  (perhaps  a  centimeter  or  more). 
Alternatively,  satisfactory  results  might  be  obtainable  if  the  plates  were  coated  on  both  sides  in 
order  to  balance  the  contractive  stresses  of  the  emulsion.  Still  another  solution,  although  an 
operationally  awkward  one,  would  consist  of  employing  a  piate  fiattener  prior  to  exposure,  as 
was  practiced  by  Schlesinger.  Such  a  procedure  has  also  been  described  by  Carman  (1968). 
Finally,  if  it  were  established  that  the  free  state  figure  of  the  plate  is  not  significantly  disturbed 
when  it  is  pressed  against  the  focal  reference,  one  could  measure  the  surface  of  the  exposed  plate 
before  development  to  generate  appropriate  corrections  to  be  applied  prior  to  the  adjustment. 
Once  perfected,  such  a  procedure  would  not  necessarily  be  excessively  burdensome. 

For  future  investigations  of  the  stellar  SMAC  calibration  of  wide  angle  mapping  cameras 
we  recommend  that  plates  be  abandoned  altogether.  Even  if  perfectly  flat,  they  constitute  an 
unnatural  intrusion  Info  the  system.  Better  would  be  the  use  of  thick  based  (preferably,  7mil) 
polyester  film  in  conjunction  with  the  actual  magazine  and  platten  associated  with  the  camera. 
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Jf  a  totai  of  20  or  more  frames  (in  the  usual  sense  of  the  word),  each  containing  about  50 
usable  images,  were  measured,  the  final  results  would  not  be  compromised  by  the  vagaries  of 
any  one  frame  (a  distinct  drawback  to  conventional  calibrations  employing  film).  Exposures 
could  be  made  with  the  camera  pointed  upward,  as  is  customarily  done  in  stellar  calibrations, 
and  with  the  camera  pointed  downward  (using  the  mirror  method  alluded  to  earlier). 

Separate  SMAC  reductions  performed  on  the  two  sets  would  establish,  in  a  definitive  manner, 
the  influence  of  extremes  of  gravity  on  the  metric  properties  of  the  camera.  Fresh  comparisons 
should  also  be  made  with  aerial  SMAC  calibrations,  preferably  based  on  photography  taken  over 
the  newly  established  Casa  Grande  Range.  Finally,  a  future  program  should  also  consider  SMAC 
calibrations  based  on  the  reduction  of  a  moderate  number  of  frames  (again,  at  least  20)  of  a 
standard  bank  of  collimators  (the  application  of  SMAC  to  laboratory  calibrations  will  be  taken  up 
later).  From  the  large  sets  of  residuals  generated  by  each  of  the  various  calibrations,  empirical 
weighting  functions  could  be  generated  and  compared,  as  could  empirical  corrections  for  residual 
systematic  errors.  Such  an  exercise  would  determine  the  relative  merits  of  the  different  approaches 
and  would  help  to  establish  the  confidence  to  be  accorded  the  various  empirical  functions  derived 
from  the  residuals.  In  Section  4  specific  recommendations  are  made  for  future  studies. 

3.2. 10  Comparisons  with  Laboratory  Calibrations 

AH  three  cameras  subjected  to  the  Aerial  SMAC  calibration  had  previously  been  cali¬ 
brated  by  Fairchild  on  a  multiccllimator  bank.  This  provided  an  opportunity  for  further  intercom- 
parison  of  results  frcm  the  different  methods  of  calibration.  In  Figure  14  we  have  plotted 
the  radial  distortion  curves  produced  by  Fairchild  and  the  corresponding  curves  resulting  from 
Aerial  SMAC  calibrations.  The  signs  of  the  original  Fairchild  curves  have  been  reversed  to  make 
the  corrections  conform  to  our  sign  convention.  Each  Fairchild  curve  represents  the  average  of 
the  four  curves  generated  by  the  four  semi-diago  n{$  of  the  fe-mat,  and  each  curve  is  accompanied 
by  an  upper  and  lower  envelope.  The  envelopes  define  boundaries  that  enclose  all  four  individual 
distortion  curves  and  thus  provide  a  measure  of  dispersion. 
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Figure  14.  Comparison  of  Aerial  SMAC  radial  distortion  curves  with  Fairchild  Multicoilimator 
curves  (broken  curves  represent  envelopes  enclosing  individual  multicoilimator  curves  along 
different  semi -diagonals). 


The  Aerial  SMAC  curve  for  camera  005  is  in  generally  good  agreement  with  the 
Fairchild  curve  and  lies  entirely  within  the  Fairchild  envelopes.  The  rms  departure  between 
the  two  curves  is  2.4  microns,  and  maximum  departure  is  about  5  microns. 

The  Aerial  SMAG  curves  for  cameras  C06  and  008  lie,  for  the  most  part,  outside  the 
Fairchild  envelopes.  Their  rms  departures  from  the  Fairchild  curves  are  4.6  and  3.3  microns, 
respectively,  and  their  maximum  departures  are  10.5  and  5.5  microns.  While  such  agreement  is 
not  especially  good,  neither  is  it  particularly  poor  by  normal  standards  of  camera  calibration. 

It  is  altogether  likely  that  the  Fairchild  calibrations  were  compromised,  as  were  our  stellar 
calibrations,  by  unflatness  of  the  photographic  plates  employed  (the  plates  used  by  Fairchild 
were  of  the  same  type  as  those  used  in  our  stellar  calibrations). 

The  values  of  the  elements  of  interior  orientation  resulting  from  the  various  methods  of 
cal  ibration  are  listed  in  Table  6.  If  corrections  for  gravity  of  30  microns  are  applied  1o  the 
principal  distances  resulting  from  Stellar  SMAC,  the  departures  of  the  various  calibrated 
principal  distances  for  cameras  006  and  008  all  become  within  10  microns  of  the  mean .  The 
departures  for  camera  005  are  within  13  microns  of  the  mean.  Such  agreement  is  consistent 
with  the  accuracies  generally  attributed  to  laboratory  calibrations. 

For  reasons  discussed  earlier,  the  results  from  Aerial  SMAC  for  principal  point  are 

partially  suspect  and  require  investigation  to  establish  their  precise  significance.  The  results 

from  Stellar  SMAC  and  Fairchild  for  principal  point  of  cameras  006  and  008  agree  to  within  50 

microns  and  22  microns,  respectively  for  x  and  to  within  14  microns  for  y  .  Except  for  the 

P  P 

50  micron  discrepancy,  this  agreement  is  within  the  range  to  be  expected  from  the  accuracies  of 
the  laboratory  method.  Here  again,  lack  of  flatness  of  the  plates  may  have  influenced  the  results. 

3.2.11  General  Conclusions  Regarding  KC-6A  Calibrations 

The  application  of  SMAC  to  the  calibration  of  KC-6A  cameras  has  not  only  demonstrated 
the  feasibility  of  the  approach,  but  also  has  shown  that  SMAC  can  produce  a  total  systems  calibiation 
to  a  degree  previously  unattainable.  Indeed,  it  has  become  clear  from  our  results  that  a  more 
pragmatic  view  of  calibration  is  In  order.  The  interior  projective  parameters  derived  from 
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a  mo  derate  number  of  frames,  are  essentially  unaffected  by  transient  systematic  errors  that 
would  be  of  consequence  to  a  calibration  based  on  any  single  frame.  Pro{ecttvely  unmodeled, 
persistent  systematic  errors  can  assume  prominence  in  the  large  sample  of  residuals  produced  by 
a  SMAC  reduclion  of  a  moderate  number  of  frames  each  containing  e  moderate  number  of  images. 
Thus  empirical  modeling  of  residual  systematic  error  can  emerge  as  a  practical  refinement  of  the 
SMAC  calibration.  So  too,  can  the  derivation  of  radial  and  tangential  weighting  functions. 

In  view  of  this  expanded  scope  of  the  process  of  calibration,  one  should  not  strive  to  suppress 
certain  sources  of  systematic  error  as  is  customarily  done.  Quite  the  contrary,  one  should 
deliberately  exercise  sources  of  persistent,  unknown  systematic  error  in  order  that  their 
influence  may  be  revealed.  For  this  reason,  we  now  appreciate  that  in  the  calibration  of 
mapping  cameras,  plates  should  not  be  substituted  for  film.  The  magazine  and  film  platen  are 
integral  parts  of  the  mapping  camera  and  their  contribution  to  systematic  errors  should  not  be, 
end  need  not  be,  slighted.  Similarly,  if  empirical  weighting  functions  are  to  have  validity 
for  general  operational  applications,  they  must  be  derived  from  observational  material 
appropriate  to  such  applications.  Thus  if  a  camera  does  not  employ  a  reseau,  corrections  for 
film  deformation  for  that  particular  camera  should  be  limited  to  what  can  be  obtained  from 
measurements  of  fiducial  marks.  Calibration  should  be  as  consistent  as  possible  with  operation. 

We  shall  postpone  making  specific  recommendations  for  future  lines  of  investigation 
into  the  metric  properties  of  the  KC-6A  cameras  until  after  we  have  explored  still  other 
applications  of  SMAC. 
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3.3  Application  of  SMAC  to  Calibration  of  Lunar  Orbiter  Cameras 


3.3.1  Background 

Having  reviewed  applications  of  SMAC  to  aerial  and  stellar  calibrations,  we  shall 
now  direct  our  attention  to  the  application  of  SMAC  to  laboratory  calibrations  of  three,  high 
resolution  Lunar  Orbiter  Cameras.  These  cameras  have  focal  lengths  of  24  inches  and  subtend 
5°  x  20°  formats  on  70  mm  film.  The  calibrations  were  performed  at  Cape  Kennedy  for  Boeing 
Company,  NASA's  prime  contractor  for  the  Lunar  Orbiter  Project.  A  key  requirement  of  the 
calibration  was  that  the  operational  environment  of  the  camera  had  to  be  simulated.  This  meant 
that  the  camera  had  to  be  operated  from  within  its  space  capsule  which,  in  turn,  had  to  be  placed 
in  a  vacuum  chamber.  These  constraints  ruled  out  any  conventional  approach  to  calibration  (e.g., 
goniometer  or  standard  collimator  bonk).  Although  the  vacuum  chamber  itself  was  fairly  small  and 
had  a  high  quality  optical  window,  it  was  not  permissible  to  remove  the  chamber  from  the  clean 
room  at  Cnpe  Kennedy.  This  ruled  out  the  possibility  of  a  stellar  calibration.  Moreover,  less  than 
a  month  was  available  from  the  approval  by  NASA  to  proceed  and  the  cutoff  date  dictated  by  the 
launch  schedule.  Only  a  minimal  laboratory  set  up  could  therefore  be  considered.  Accordingly, 
NASA  and  Boeing  decided  to  settle  on  a  calibration  leading  primarily  to  a  precise  determination 
of  focal  length;  all  other  potentially  obtainable  results  were  to  be  considered  to  be  of  secondary 
importance.  This  guideline  led  to  an  experimental  set  up  limited  to  photographs  of  a  single  pair 
of  collimators  subtending  a  precisely  known  angle  of  nominally  10°  (or  about  half  the  horizontal 
field). 

Inasmuch  as  up  to  eighteen  frames  could  be  allotted  to  each  camera  for  the  calibration, 
we  suggested  a  SMAC  calibration  wherein  the  pair  of  collimators,  mounted  on  a  horizontal  rotary 
toble,  would  be  rotated  about  one  half  degree  between  exposures  on  successive  frames  so  that  the 
pair  of  collimator  images  would  proceed  from  one  extremity  of  the  horizontal  field  to  opposite 
extremity.  Thus  if  all  frames  were  superimposed  to  produce  a  single  composite  frame,  the  eighteen 
pairs  of  collimator  images  would  be  fairly  uniformily  distributed  across  the  format  with  all  images 
very  nearly  lying  on  a  common  straight  line  passing  through  the  center  of  the  frame  in  the  long 


(or  x)  direction  of  the  format.  Clearly,  a  more  desirable  distribution  would  have  been  one  with 
images  arrayed  across  both  diagonals  of  the  format.  Unfortunately  this  could  not  be  considered, 
in  view  of  the  time  available,  because  of  the  severe  physical  constraints  imposed  on  the  operation. 

A  computer  simulation  of  the  SMAC  reduction  appropriate  to  proposed  experiment  showed 
that  focal  length  could  theoretically  be  recovered  to  an  accuracy  (one  sigma)  of  about  30  microns. 
Thus,  the  desired  result  was  potentially  attainable.  The  simulation  also  showed  that  coefficients 
of  radial  distortion  (Kj ,  K2)  and  one  of  the  coordinates  (xp)  of  the  principal  point  could  be 
recovered  to  meaningful  accuracies.  The  second  coordinate  of  the  principal  point  yp  was  found 
to  be  inherently  unrecoverable  because  the  distribution  of  control  points  produced  no  variation 
in  the  y  coordinates  of  the  images.  The  recovery  of  decentering  distortion  could  not  be  effected 
because  it  was  found  to  require  a  minimum  of  three  control  points  per  frame  as  well  as  variation 

l 

in  both  x  and  y  coordinates.  Despite  the  fact  that  only  a  partial  calibration  could  be  achieved, 
the  results  theoretically  obtainable  were  considered  sufficient  to  warrant  implementation  of  the 
experiment. 

3.3.2  Results 

The  experiment  was  successfully  executed  and  the  results  turned  out  to  be  consistent 
with  expectations.  They  are  reported  in  detail  in  Brown  (1967a).  Table  8  and  Figure  15 
below  are  reproduced  from  this  reference.  Our  concern  here  is  not  so  much  with  results  themselves 
as  with  the  implications  of  the  results  to  laboratory  calibrations. 

Although  a  pair  of  collimators  was  used  for  the  calibration,  the  reduction  appropriate 
to  the  experiment  was,  in  fact,  a  Stellar  SMAC  reduction-  The  pivotal  consideration  is  that  a 
pair  of  collimators  maybe  viewed  as  a  pair  of  artificial  stars  subtending  a  known  an'jle  a?*. 

One  is  therefore  at  liberty  to  assign  one  collimator  coordinates  in  right  ascension  and  declination 
of  (0,0)  and  the  second  collimator  coordinates  (Qf*,0).  Equally  admissible,  in  fact,  would  be 
assignment  of  any  arbitrary  pair  of  stellar  coordinates  (a  \ ,  8] ),  ( <*2,82 )  satisfying  the  relation 

cos  a1*  =  sin  8j  sin62  +  cosSj  cosS2  cos(ctj  -o?2). 
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fable  8.  Results  of  dual  collimator  SMAC  calibration  of  Lunar  Orbiter  Cameras 
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Now,  although  in  the  experiment  itself  the  pair  of  collimators  was  actually  rotatea*  slightly 
from  one  exposure  to  the  next,  one  could  equally  well  regard  the  collimators  as  remaining 
fixed  and  the  camera  as  having  undergone  the  rotation.  Under  this  interpretation,  the  Stellar 
SMAC  reduction  becomes  immediately  and  directly  applicable  to  calibrations  based  on  a  series 
of  photographs  of  two  or  more  collimators. 

Considering  the  skimpiness  of  the  observational  material,  the  results  produced  by  SMAC 
for  the  Lunar  Orbiter  Cameras  are  indeed  impressive.  It  should  be  appreciated  that  on  any  given 
frame  the  observations  consisted  of  merely  the  plate  coordinates  (x,y)  of  each  of  a  pair  of  artificial 
stars  (i.e.,  the  collimator  images).  These  coordinates  give  rise  to  four  observafiona!  equations 
from  the  i^  frame  involving 

a)  the  determinable  parameters  of  the  inner  cone,  namely: 

« 

(i)  c,  the  principal  distance; 

(ii)  Xp,  the  horizontal  component  of  the  principal  point; 

(Hi)  K|,'K2,  coefficient  of  radial  distortion; 

b)  the  angular  elements  of  orientation  for  the  frame:  Of;,  coj ,  k'j  . 

Consistent  with  the  SMAC  concept,  the  parameters  under  group  (a)  are  considered  to  be  common 
to  all  frames,  whereas  those  under  group  (b)  are  considered  to  vary  from  frame  to  frame  inasmuch 
as  the  orientation  of  the  camera  in  the  'stellar'  coordinate  system  is  altered  between  exposures. 
Hence,  three  of  the  four  available  equations  from  a  given  frame  are,  in  effect,  used  up  in 
establishing  the  angular  elements  of  orientation  peculiar  to  that  frame;  this  leaves  only  one 
redundant  observation  per  frame  to  be  applied  to  the  recovery  of  the  parameters  of  the  inner  cone 
that  are  common  to  all  frames.  All  told,  the  eighteen  available  frames  generate  a  total  of 
4x  18  =  72  observational  equations  involving  a  total  of  57  or  58  unknowns  consisting  of  3 x  18  =  54 
angular  elements  of  orientation  (three  per  frame)  plus  3  or  4  parameters  of  the  inner  cone  oepending 
on  whether  one  or  two  coefficients  are  needed  for  the  redial  distortion  function.  Inasmuch  as  a 
single  coefficient  turned  out  to  be  sufficient,  the  simultaneous  adjustment  of  all  observations 
generates  a  general  system  of  normal  equations  of  order  57,  involving  14  statistical  degrees  of 
freedom.  This  is  hardly  a  strong  adjustment,  but  it  is  the  best  that  can  be  done  with  the  limited 
obrr  ational  material  and  is,  as  Table  8  and  Figure  15  indicate,  sufficient  to  produce  useful 
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results.  Although  the  Lunar  Orbiter  experiment  did  not  permit  the  recovery  of  the  yp  or  of 
parameters  of  decentering  distortion,  the  converse  proposition  holds  that  these  unrecoverable 
parameters  of  the  inner  cone  can  have  no  effect  on  the  values  obtained  for  the  recoverable 
parameters  of  the  inner  cone.  Hence,  the  results  of  the  calibration  are  in  no  way  contaminated 
by  the  enforced  values  (zero)  of  the  unrecoverable  parameters. 

3.4  General  Application  of  SMAC  to  Multicoliimator  Calibrations 
3.4.1  Application  to  Conventional  Collimator  Banks 

We  have  seen  that  the  pair  of  collimators  employed  in  the  Lunar  Orbiter  calibrations 
was  successfully  treated  in  the  SMAC  reduction  as  a  pair  of  artificial  stars.  If  follows  that  a 
bank  of  collimators  employed  in  conventional  laboratory  calibrations  may  also  be  regarded  as 
an  artificial  star  field.  Accordingly,  sets  of  photographic  observations  of  mu'ticollimators  may  . 
be  directly  processed  through  a  Stellar  SMAC  reduction.  In  conventional  multicoliimator  cali¬ 
brations,  a  time  consuming  and  painstaking  process  of  autocollimation  must  be  performed  in  order 
to  orient  the  camera  so  that  the  axis  of  the  central  collimator  is  precisely  perpendicular  to  the 
focal  plane.  This  process  can  be  totally  bypassed  when  a  Stellar  SMAC  reduction  is  performed; 
here,  it  is  sufficient  merely  to  position  the  camera  so  that  its  entrance  pupil  intercepts  the  con¬ 
verging  beams  of  the  collimators.  This,  in  turn,  makes  it  practical  to  avoid  the  use  of  plates 
(if  desired)  by  exercising  the  film  magazine  (or  magazines)  of  the  camera  in  the  photography  of  the 
collimator  bank.  Because  precise  alignment  is  not  required,  a  large  number  of  frames  can  be 
quickly  exposed  with  the  camera  being  rotated  about  its  axis  between  frames.  This  procedure  would 
generate  a  very  large  number  of  collimator  images  having  an  overall  distribution  far  superior  to 
that  obtainable  from  ony  one  frame.  Carried  fn  a  SMAC  reduction,  the  measurements  of  all  such 
images  would  contribute  to  a  common  calibration  of  the  interior  projective  parameters.  The 
large  sample  of  measuring  residuals  that  could  thus  be  generated  would  provide  the  material 
necessary  for  sound  empirical  modeling  of  residual  systematic  error  and  for  the  derivation  of 
weighting  functions  appropriate  to  the  camera. 


From  the  foregoing,  it  is  clear  that  by  applying  Stellar  SMAC  to  multicoliimator 
calibrations,  one  can  not  only  greatly  simplify  set  up  procedures  but,  by  admitting  the  use  of 
film  and  the  exercise  of  several  frames  at  different  swing  angles,  one  can  also  expect  to  obtain 
a  more  meaningful,  more  accurate  and  more  thorough  calibration.  Such  a  calibration  could 
be  extended  to  take  into  account  errors  in  the  directions  of  the  collimators.  In  principle,  it 
would  suffice  if  the  angle  subtended  by  but  a  single  pair  of  collimators  v/ere  accurately  known; 
the  directions  of  all  other  collimators  could  then  be  established  within  the  reduction  itself. 

3.4.2  Implications  of  SMAC  to  the  Design  of  New  Multicoilimators 

As  we  have  just  seen,  SMAC  is  directly  applicable  to  the  reduction  of  photographs 
of  conventional  multicoilimators.  The  first  of  the  modern  collimator  banks,  designed 
and  constructed  by  the  National  Bureau  of  Standards  in  1949,  ii.corporates  a  total  of  25  colli¬ 
mators  arrayed  to  span  both  diagonals  of  the  photographic:  format.  A  multicollimator  later 
developed  for  the  U.S.  Geological  Survey  is  of  similar  general  design  but  incorporates  almost 
twice  as  many  (49)  collimators.  A  multicoliimator  recently  developed  by  the  National  Research 
Council  in  Canada  employs  a  total  of  43  collimators  evenly  spaced  at  increments  of  2p8l25 
and  arrayed  across  a  single  diagonal.  A  multicoliimator  soon  to  be  developed  for  the  U.  $.  Air 
Force  will  incorporate  about  150  collimators. 

It  is  clear  from  the  foregoing  that  the  trend  is  towards  larger  and  targer  collimator  banks 
containing  more  and  more  collimators.  It  is  thus  also  clear  that  little,  if  any,  consideration  has 
been  given  to  the  possibility  that  a  superior  calibration  could  conceivably  result  from  a  combina¬ 
tion  of  fewer  collimators  and  a  more  advanced  approach  to  data  reduction.  We  demonstrated  in 
the  Lunar  Orbiter  calibrations  that  if  decentering  distortion  is  ignored  (as  indeed  it  is  in  conven¬ 
tional  multicoliimator  calibrations)  a  single  pair  of  collimators  is  sufficient  for  the  accomplish¬ 
ment  of  a  SMAC  calibration.  The  recovery  of  decentering  distortion  requires  that  a  third  colli¬ 
mator  be  added.  In  principle,  then,  a  minimum  of  three  colli  mater  can  provide  a  full  calibra¬ 
tion  of  a  camera,  provided  that  a  moderately  large  number  of  frames  is  c-rposed  and  a  good  distri¬ 
bution  of  images  is  realized.  While  more  than  three  can  provide  a  superior  result,  it  does  not 
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follow  that  a  single  frame  containing  a  large  number  of  collimator  images  is  necessarily  to  be 
preferred  over  several  frames  each  recording  a  relatively  small  number  of  collimator  images. 

Indeed,  o,uite  the  opposite  could  be  true.  We  shall  consider  here  how  the  computational 
tradeoffs  made  possible  by  ihe  SMAC  reduction  might  be  exploited  in  the  design  of  a  simpii*- 
fied  multicollimator  capable  of  producing  superior  results.  A  SMAC  multicoliimator  suitable 
for  the  calibration  of  six  inch  mapping  cameras  is  shown  schematically  in  Figure  16.  The  nine 
collimators  in  the  array  are  considered  to  be  of  a  diffraction  limited,  catadiopfric  type  of 
moderate  (about  10  cm)  aperture.  As  shown  in  Figure  16a,  the  axes  of  the  collimators  are  all 
horizontal  and  all  are  nominally  in  the  same  plane.  The  camera  is  mounted  vertically  (i.e., 
in  its  normal  operating  position), and  the  collimalor  rays  are  reflected  to  the  entrance  pupil 
by  means  of  an  adjustable,  optically  fiat  mirror  (Figure  16b).  The  mirror  tilts  about  a  horizontal 
axis  that  is  nominally  perpendicular  to  the  axis  of  the  cenlrcl  collimator.  With  each  tilt  of  the 
mirror,  the  point  of  intersection  of  the  bundle  of  collimator  rays  assumes  a  different  position  * 
along  the  arc  PQ  (Figure  b).  A  simple  linkage  between  the  mirror  support  and  the  camera  plat¬ 
form  can  be  provided  to  reposition  the  camera  automatically  so  that  the  entrcnce  pupil  is  always 
at  the  proper  location  and  the  camera  axis  is  always  maintained  In  a  nominally  vertical  orienta¬ 
tion.  In  a  typical  calibration  ,  eleven  exposures  of  the  multicoliimator  would  be  mcde  with  the 
mirror  tilted  as  follows:  -18°,  >16°,  -12°,  -8°,  -4°,  0,  4°,  8°,  12°,  16°,  18c,  in  which  0° 
corresponds  to  the  case  in  which  the  mirror  is  inclined  45°  to  the  horizontal.  These  angles  need 
not  be  established  with  great  accuracy.  The  composite  frame  for  K  -  0°  would  have  the  appearance 
indicated  in  Figure  17a.  For  higher  accuracies,  the  process  could  be  repeated  at  different  swing 
angles  (e.g.,K  =  90°,  180°,  270°).  The  overlay  of  ell  exposures  at  K  =  0°  and  K  =  90°(Fig.  17b), 
would  generate  the  pattern  of  images  shown  in  Figure  17c.  In  slill  another  mode  of  operation,  one 
could  produce  a  supplemental  radially  symmetric  distribution  of  images  by  holding  the  camera 
fixed  in  its  midposition  of  0°  end  altering  the  swing  angle  of  the  camera  between  exposures.  Figure 
17d  shows  the  composite  pattern  of  108  images  that  would  be  generated  by  this  procedure  if  incre¬ 
ments  of  15°  in  swing  angle  were  exercised  over  the  range  K  =  0  to  K  -  165°.  From  these 
examples  it  is  clear  that  the  suggested  SMAC  Multicoliimator  can  provide  a  most  satisfactory 
overall  distribution  of  images  from  a  moderate  numbet  of  exposures.  The  addition  of  two  pairs  of 
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Figure  16.  Illustrating  concept  of  SMAC  Multicoliimator  for  calibration  of  mapping 
cameras. 


Fig.  17a.  &-(),  x  ~  0/  a:-  36° ,32° ,24’ ,  16%  Fig.  17b.  Same  as  Fig.  16a,  except-  x=90° 
8° ,  0° , -8° , - 16° ,  -24° , -32° , -36° 


Fig.  17c.  Overlay  of  Figs.  17a  and  17b. 


Fig.  I7d.  a=0°,  a>=0\x=Q,  15°, 30°,... 


Figure  17.  Illustrating  some  of  the  various  patterns  of  control  that  can  be  generated  by 
superimposing  exposures  of  SMAC  muiticollimator  illustrated  in  Fig.  16. 
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collimators  to  the  array  shown  in  Figure  16a  would  make  the  instrument  well  suited  to  the  cali¬ 
bration  of  super  wide  angle  cameras. 

Relative  to  conventional  multicollimators  subjected  to  conventional  data  reduction, 
a  specifically  designed  SMAC  Multicollimator  has  a  number  of  advantages:  . 


(a)  Because  of  the  small  number  of  collimators  involved  and  because  of  their  simple 
arrangement,  the  cost  of  the  instrument  would  be  relatively  low. 

(b)  Accordingly,  it  would  be  possible  to  invest  appreciably  more  in  collimators  of 
wider  aperture  and  superior  optical  performance. 

(c)  The  relative  simplicity  of  a  SMAC  Multicollimafor  would  facilitate  maintenance 
and  adjustment  of  the  instrument. 

(d)  Despite  the  fact  that  several  exposures  must  be  produced,  appreciably  less  time 
would  be  required  for  data  gathering  operation  becausfe  the  tedious  process  of 
autocollimation  is  dispensed  with. 

(e)  By  virtue  of  (c)  the  results  of  the  calibration  are  in  no  way  compromised  by  lack 
of  precision  in  set  up  procedures  (thus  the  'point  of  symmetry',  an  artifice  used 

to  compensate  in  part  for  the  deficiencies  of  conventional  procedures  ,  has  no  role 
in  a  SMAC  Multicollimator  calibration). 

(f)  Also  by  virtue  of  (c)  it  becomes  practical,  when  desired,  to  use  film  instead  of 
plates,  thereby  allowing  the  camera  magazine  and  platen  to  influence  the  results. 

(g)  As  in  ony  SMAC  calibration,  there  is  no  limit  to  the  number  of  frames  that  can  be 
processed  in  a  simultaneous  reduction;  thus  specific  schedules  of  exposures  can  be 
tailored  to  produce  specific  levels  of  accuracy. 

(h)  By  exercising  a  straightforward  extension  of  the  SMAC  reduction,  one  could  on 
occasion  accurately  determine  (and  thus  check)  the  directions  of  all  collimators 
relative  to  the  adopted  angle  between  an  arbitrary  pair  of  collimators;  thus  periodic 
physical  monitoring  of  collimator  angles  could  be  limited  to  measurements  of  the 
angle  between  a  selected  (preferably  widespaced)  pair  of  collimators. 

(i)  If  the  two  selected  collimators  were  designed  to  function  also  as  autocollimators, 
the  adjustment  of  their  subtended  angle  to  a  desired  value  0  could  be  rapidly 
and  accurately  accomplished  by  autocollimation  on  the  faces  of  a  prism  subtending 
an  angle  of  180°  ~Q  .  Thus  a  single  calibrated  prism  could  provide  the  ultimate 
standard  of  angular  measure  for  a  SMAC  Multicollimator. 


The  concept  of  the  SMAC  Multicollimator  provides  still  another  example  of  how  hard¬ 
ware  can  be  simplified  end  improved  when  sophisticated  techniques  of  data  reduction  are  allowed 
a  major  role  in  establishing  design.  Previous  examples  of  this  approach  include  the  Multilaterafive 
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Plato  Measuring  Comparator  (Brown.  1967c,  1968c)  and  the  Satellite  Surveying  Utility 
(Brown,  1968b).  Further  examples  are  discussed  in  Brown  1967b. 

3.5  Application  of  SMAC  to  Analysis  of  Ballistic  Camera  Stability  ; 

3.5.1  Introduction 

We  have  already  noted  that  by  exercising  SMAC  in  the  stellar  calibration  of  cameras, 
one  eliminates  all  requirements  for  precise  timing  of  exposures  and  for  physical  stability  of  the 
camera  over  the  series  of  exposures.  On  the  other  hand,  when  precise  timing  of  exposures  is 
available,  the  application  of  SMAC  enables  one  to  determine  with  great  accuracy  the  physical 
stability  of  the  camera.  Such  an  application  of  SMAC  has  particular  merit  in  ballistic  camera 
operations,  for  here  stability  is  all  important.  Accordingly,  the  evaluation  of  stability  that 
automatically  emerge:  from  the  SMAC  calibration  of  a  ballistic  camera  provides  information  . 
of  considerable  value. 

3.5.2  Example 

To  illustrate  the  concept  of  SMAC  analysis  of  stability,  we  shall  consider  the  results 
of  an  experiment  performed  on  observations  made  by  a  PC-1000  camera  in  December  1966  at 
Goddard  Space  Flight  Center.  Zenithal  exposures  were  recorded  with  the  camera  placed  on  a 
concrete  pad  .  The  exposures  were  precisely  timed  and  were  made  at  five  minute  intervals 
over  a  period  of  more  than  half  an  hour.  In  the  SMAC  reduction,  a  common  set  of  interior 
projective  parameters  (xp,  yp,  c,  Kj ,  K 2,  Pj ,  P2)  was  recovered  for  eight  successive  exposures, 
and  a  separate  set  of  angular  elements  was  recovered  for  each  of  the  eight  exposures.  The  adopted 
angular  elements  consisted  of  the  hour  angle  and  declination  (H,  6)  of  the  camera  axis  and  the 
swing  angle  i(  referred  to  the  line  of  intersection  of  the  photographic  plate  and  the  equatorial 
plane.  Approximately  25  stellar  images  from  each  exposure  were  carried  in  the  SMAC  reduction. 

Our  interest  here  is  not  with  the  interior  projective  elements  resulting  from  the  reduction, 
but  rather  with  the  variation  displayed  by  the  angular  elements.  The  variation  in  each  of  the 
angular  elements  about  its  mean  is  indicated  in  Figure  18  .  The  plotted  results  for  hour  angle  have 
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been  normalized  by  the  customary  process  of  multiplication  by  the  cosine  of  declination. 

Each  point  is  accompanied  by  a  vertical  bar  defining  its  plus  and  minus  one  sigma  confidence 
intervals.  The  one  sigma  values  for  normalized  hour  angle  and  declination  are  slightly  less 
than  0.2  seconds  of  arc;  for  sv'ing  angle,  they  are  generally  about  ten  times  greater,  averaging 
close  to  2  seconds  of  arc.  This  disparity  in  sigmas  is  attributable  to  the  fact  that  the  focal  length 
(1000  mm)  of  the  PC-1000  is  about  10  times  greater  than  the  semi-diagonal  of  the  plate  format. 

The  projective  effect  of  an  error  of  0.2  seconds  of  arc  in  the  direction  of  the  camera  axis  is 
equivalent  to  that  of  an  error  of  about  one  micron  on  the  plate.  By  the  same  token,  the  projective 
effect  of  an  error  of  2  seconds  of  arc  inswing  angle  is  equivalent  to  that  of  an  error  of  about  one 
micron  near  the  edge  of  the  plate.  Thus,  there  is  no  actual  projective  disparity  in  the  relative 
sigmas  of  the  angular  elements. 

From  Figure  18  we  see  that  several  of  the  variations  in  orientation  differ 
significantly  from  zero  at  the  two  sigmo  level.  Specifically,  these  consist  of 

(a)  hour  angles  at  t  -  5,  10,  20,  30,  and  35  minutes; 

(b)  declinations  at  t  =  0,  20,  25,  30,  and  35  minutes; 

(c)  swing  at  t  =  20  minutes. 

A  few  extreme  excursions  in  the  direction  of  the  camera  axis  are  seen  to  amount  to  as  much  as 
two  seconds  of  arc  in  a  five  minute  interval.  For  the  most  part,  departures  about  the  fitted 
straight  lines  in  the  figures  are  seen  to  be  appreciably  less  than  the  departures  about  the  means. 

This  suggests  the  existence  of  slow  secular  trends  upon  which  are  superimposed  small,  short  term 
changes  that  occasionally  assume  significance.  With  only  a  few  exceptions,  moving  arc  linear 
functions  can  accommodate  variations  over  ten  minute  intervals  to  accuracies  consistent  with  those 
of  the  plotted  points. 

Many  experiments  along  the  above  lines  would  have  to  be  performed  in  order  to  establish 
solid  conclusions.  The  present  results  indicate  that  significant,  nonabrupt  changes  in  orientation 
can  occur  over  fairly  shoit  periods  of  time.  For  this  reason,  in  routine  tracking  operations  the  total 
sequence  of  exposures  should  be  limited  to  as  shoit  an  interval  as  possible  consistent  with  acquisition 
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of  o  sufficient  number  of  control  points.  A  suitable  observational  schedule  for  satellite  observa¬ 
tions  is  recommended  in  Brown  (1964);  it  spans  a  period  of  |ust  over  five  minutes  and  allows  each 
star  to  generate  four  sets  of  punctiform  images,  it  is  true  that  the  resuits  of  our  experiment  indicate 
that  camera  orientation  can  vary  significantly  over  a  five  minute  period.  However,  it  is  aiso  true 
that  the  variation  over  such  a  short  interval  is  generally  likely  to  be  sufficiently  linear  to  admit 
the  representation  of  the  angular  elements  as 

•  a  =  a0  +  (t  -  t0)  , 

0)  =  to0  +  to0  (t  -  tj  , 

K  =  K0  +  Ke  (t  -  t0)  , 

in  which  t  denotes  the  time  of  the  exposure,  t0  is  the  mean  time  of  all  exposures,  and  d'c ,  to0 , 

K0  are  the  unknown  rates  of  change  of  the  angular  elements  about  their  respective  means  <yQ,  to6, 
K0.  The  advanced  plate  reduction  developed  in  Brown  (1964)  can  readily  be  extended  to  incorpor- 
are  6t0  ,  to0  ,  k0  as  additional  parameters  to  be  recovered.  Our  experience  over  the  past  few 
years  indicates  that  such  an  extension  would  be  especially  worthwhile  if  accompanied  by  a  statistical 
testing  procedure  designed  to  suppress  6fot  toc  ,  or  k0  automatically  to  zero  if  it  were  found  to 
differ  insignificantly  from  zero.  In  such  a  reduction,  instantaneous  orientation  matrices  generated 
by  the  a,  to,  and  k's  corresponding  to  times  of  satellite  observations  would  be  used  in  computing 
satellite  directions. 

We  believe  that  instability  is,  in  many  instances,  the  predominant  source  of  error  in 
PC-1000  observations.  When  virtually  perfect  stability  is  maintained  throughout  the  observational 
period  and  average  atmospheric  shimmer  is  experienced,  one  can  expect  the  PC-1000  to  produce 
directional  accuracies  of  about  0.6  to  0.7  seconds  of  arc  for  individual  flashes  from  a  satellite 
(Brown,  1966).  Indications  are,  however,  that  this  level  of  accuracy  is  actually  achieved  in 
general  operations  only  with  about  one  plate  in  thiee,  the  others  being  compromised  to  a  small, 
but  significant  degree  by  instability.  The  extended  reduction  discussed  above  would,  in  most 
cases,  "ndo  the  damage  induced  by  instability  and  would  thus  upgrade  PC-1000  accuracies. 

The  implementation  of  such  a  reduction  is  therefore  to  be  recommended.  In  addition,  a  statistical 
study  of  instability  by  means  of  SMAC  reductions  of  a  large  number  of  plates  exposed  over  extended 
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periods  should  be  undertaken  to  establish  a  clearer  undemanding  of  the  problem  and  to  generate 
realistic  a  priori  constraints  to  bs  exercised  for  the  angular  rates  6t0  ,  c0o  ,  . 


CONCLUSIONS  AND  RECOMMENDATIONS 


Our  earlier  statement  concerning  the  flexibility  and  universality  of  the  SMAC  reduction 
is  well  supported  by  the  various  examples  given  in  Section  3-  We  hare  seen  that  SMAC  is  equally 
applicable  to  aerial  calibrations,  to  stellar  calibrations,  and  to  laboratory  calibrations*  In  a 
special  application  to  laboratory  calibrations,  we  have  seen  that  SMAC  can  provide  the  basis  for 
the  design  of  much  simplified,  yet  possibly  more  effective,  multicollimcitors* 

Because  SMAC  con  be  exploited  as  a  common  method  of  reduction  for  all  of  the  basic 
approaches  to  calibration,  it  provides  the  logical  basis  for  the  intercomparison  of  different  pro¬ 
cesses,  We  recommend  therefore  that  further  studies  be  undertaken  to  determine  the  precise  degree 

l 

of  validity  of  various  approaches  to  laboratory  and  stel!  ur  calibration  by  employing  comprehensive 
operational  Aerial  SMAC  calibrations  as  standards.  The  limited  comparisons  made  so  far  suggest  . 
that  conventional  laboratory  calibrations  of  radial  distortion  CGn  be  effectively  in  error  by  as  much 
as  ten  microns  near  the  corners  of  the  format.  Some  of  this  error  may  well  be  attributable  to  insuffi¬ 
ciently  flat  plates,  some  to  the  small  number  of  available  control  points  (i.e.,  small  relative  to 
the  expanded  number  made  possible  through  a  SMAC  reduction),  some  to  errors  in  alignment  of  the 
camera,  and  some  to  the  method  of  data  reduction.  We  suggest  that  in  future  tests  laboratory 
calibrations  be  performed  not  only  in  the  usual  manner,  but  also  in  c  manner  appropriate  to  a 
SMAC  reduction.  Inasmuch  as  further  flight  testing  of  the  USQ-28  System  is  now  being  planned, 
we  recommend  that  one  phase  of  the  testing  be  devoted  to  a  thorough  comparative  analysis  of  as 
many  different  methods  of  calibration  os  is  practicable. 

Specifically,  wc  recommend  that  a  pair  of  KC-6A  Cameras  be  subjected  to  the  following 
testing  program: 

A.  Laboratory  Calibrations: 


(1)  Sfanchic!  calibrations  on  at  hast  three  different  nvul'icoliimaiors  (e.g., 
those  at  National  Bureau  of  Standards,  U.  $,  Geological  Survey  and 

r-  •  -i*i  j\ 
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(2)  SMAC  calibrations  with  some  multicoliimators  as  (1)  using  film  and  a  series 
of  exposures  on  different  frames  at  15°  increments  in  swing  angle  from 

K  =  0°  to  K  =  345°  (a  total  of  24  frames  per  camera). 

B.  Stellar  SMAC  Calibrations 

(1)  Stellar  calibration  with  camera  facing  upward,  using  film  and  the  following 
operational  sequence;  on  the  first  frame,  expose  star  trails  for  60  seconds, 
close  the  shutter  for  30  seconds,  expose  stars  for  one  second;  advance  to 
next  frame,  rofare  canieta  15°  in  swing  angle  and  repeat  same  exposures 

os  first  frame;  continue  in  this  manner  at  15°  increments  in  swing  angle 
until  a  total  of  24  frames  has  been  exposed. 

(2)  Same  procedures  as  above  but  with  camera  facing  downward  into  large, 
optically  flat  mirror  (at  least  24"  in  diameter)  or  a  large,  well  damped 
mercury  poo! . 

C.  Aerial  SMAC  Calibrations:  Clover  leaf  flights  over  McLure  or  Casa  Grande 

Ranges  in  accord  v/ith  following  specifications:  * 

(la)  Night  flight,  illuminated  targets,  ballistic  camera  tracking,  Shiran  tracking, 
cycling  rare  of  one  frame  per  three  seconds,  image  Motion  Compensation 
exercised,  positive  x  axis  of  film  always  nominally  aligned  with  direction 
of  flight. 

Same  as  (la) ,  but  performed  on  a  different  night. 

Same  as  (la),  but  with  IMC  not  exercised. 

Same  as  (la),  but  with  camera  rotated  180°  in  swing  angle  between  direct 
and  reverse  runs  on  each  path . 

Daytime  flight,  no  ballistic  camera  tracking  or  illuminated  targets,  but 
otherwise  some  as  (la). 

On  each  of  these  five  fests/24  frames  (6  from  each  flight  leg)  are  to  be  used 
in  the  Aerial  SMAC  Calibration. 


.  (Ib) 
(2) 

(3) 

(4) 


The  wealth  of  material  gathered  in  such  a  testing  program  would  permit  a  most  thorough  evalua¬ 
tion  of  the  various  processes  of  camera  calibration.  Each  of  the  SMAC  calibrations  would  yield 
a  sufficiently  largo  sJmple  of  residuals  to  permit  both  the  extraction  of  empirical  models  for  residual 
systematic  error  and  the  determination  of  empirical  weighting  functions.  The  degree  of  fesr-to-fes! 
consistency  obtained  for  such  empirical  models  and  weighting  functions  v/ould  help  to  establish  the 
degree  of  confidence  to  be  placed  in  these  concepts. 
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In  all  of  the  above  tests  involving  film,  the  first  ten  frames  of  the  roll  should  be 
expended  immediately  prior  to  the  data  gathering  exposures  in  accordance  with  recommenda¬ 
tions  made  by  Carman  (1968).  Carman  found  that  the  first  few  frames  on  a  roll  may  be  subject 
to  excessive  film  deformation  attributable  to  effects  of  atmospheric  humidity.  Normal  film 
deformation  should  present  no  problems  in  the  recommended  testing  program,  for  if  can  be 
thoroughly  compensated  for  by  the  platen  reseau  of  the  KC-6A. 

If  if  could  be  Gr ranged,  we  would  recommend  that  still  another  operational  approach 
to  camera  calibration  be  evaluated  in  conjunction  v/ifh  the  above  resting  program.  This 
approach  involves  the  recovery  of  radial  and  decenfering  distortion  parameters  within  the 
process  of  aerotriangulatiori.  An  adequate  test  could  be  based  on  a  49  photo  block  consisting 
of  7  strips  each  containing  7  photos  with  60%  forward  overlap  and  60%  side  overlap.  A  pattern 
of  25  pass  points  v/ould  be  measured  on  each  photo.  All  absolute  control  would  be  withheld, 
except  for  the  minimum  needed  for  determinacy.  As  is  discussed  in  Brown  (1968a),  the  normal  ' 
equations  for  the  simultaneous  adjustment  of  conventional  photogrammetric  blocks  can  be  made 
to  assume  a  patterned,  bana'ed-bordered  form  as  indicated  in  Figure  19  .  Although  the  system 
of  normal  equations  for  the  proposed  experiment  wcula'  involve  the  simultaneous  recovery  of 
6  x  49  =  294  unknown  elements  of  exterior  orientation  plus  perhaps  as  many  as  6  distortion 
parameters,  its  solution  can  be  effected  with  extraordinary  efficiency  by  means  of  a  special 
algorithm  cal'ed  Recurrent  Partitioning.  If  it  should  turn  out,  as  we  suspect  might  well  be  the 
case,  that  a  satisfactory  camera  calibration  can  be  produced  as  an  incidental  byproduct  of 
analytical  aerotriangulafion  of  a  suitably  des’gned  block,  such  a  process  might  ultimately  emerge 
as  the  method  generally  to  be  preferred  in  aerial  phofogremroetry.  In  view  of  this,  we  urge  that 
calibration  by  aerotriangulafion  be  tried  and  evaluated  as  pari  of  the  recommended  testing  pro¬ 
gram. 

Our  expectation  of  achieving  success  with  the  method  just  described  is  based  on  the 
fact  that  over  the  past  six  years  we  have  successfully  employed  a  variation  of  if  with  plate 
cameras  used  for  photogrammetric  trianguiatlon  of  points  cn  large  structures  such  as  radio  tele¬ 
scopes.  This  application  differs  from  the  proposed  application  to  aerial  cameras  in  that  (a)  highly 
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convergent  geometry  is  exercised,  (b)  only  from  three  to  five  plates  having  100  percent  overlap 
are  carried  in  the  reduction,  and  (c)  a  minimum  of  about  a  hundred  pass  points  per  plate  are 
employed  in  the  analytical  aerotiiangulatior..  The  proposed  test  would  establish  whether  or  not 
equally  satisfactory  results  could  be  obtained  from  a  block  of  vertical  aerial  photographs  having 
fev/er  points  per  photo  (e.g.,  the  25  suggested),  less  overlap,  but  many  more  photos  (e.g., 
the  49  suggested). 

Our  final  recommendation  is  that  the  concept  of  the  SMAC  Multicollimator  be  explored 
more  deeply  in  future  investigations.  In  particular,  a  computer  simulation  should  be  performed  to 
determine  the  minimum  number  of  frames  and  distribution  of  images  needed  for  the  self-calibration 
of  the  directions  of  the  collimators  relative  to  the  adopted  angle  between  a  selected  pair  of  colli¬ 
mators.  Other  simulations  should  be  performed  to  ascertain  accuracies  to  be  expected  from  various 
distributions  of  images.  Should  the  outcome  of  such  computer  studies  be  favorable,  we  would 

I 

recommend  that  the  approach  be  tested  further  by  means  of  a  temporary  laboratory  set  up  of  mini¬ 
mal  cost.  Such  an  investigation  is  especially  warranted,  we  feel,  by  the  fact  that  a  large  scale 
multicollimator  currently  being  considered  by  the  Air  Force  could  well  cost  close  to  a  million 
dollars.  In  our  view,  a  SMAC  Multicollimator  costing  a  small  fraction  of  this  amount  could 
conceivably  produce  equal,  if  not  superior,  results. 
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appendix  a 


LINEARIZATION  OF  THE  PROJECTIVE  EQUATIONS 


We  assume  that  possibly  nonzero  approximations  are  available  for  the  distortion 
parameters  and  set 

x  =  x  H  x(K?°<*  +  K|°r*  +K°°rfi) 


(A.  1) 

A 

y 


+  [P°°(r3+2^)  +  2P|°xy][l  +P|°rs] 


y +y(K°°  r3  +  K°°  r4  +  K^i*) 


+  [2P°  0 x y  +  P°  0  (rs  +  2y3 )  ][  1  +  P°  0  ra  3 


I 


in  which 

x  —  x°  -  x°® 

p 

(A. 2)  y  =  y°  -  y°° 

r  =  (if3  +”y2}^. 


The  quantities  x,y  then  represent  values  of  the  measured  plate  coordinates  referred  to  the 
approximate  principal  point  and  corrected  approximately  for  radial  and  decentering  distortion. 

The  orientation  matrix  to  be  used  in  the  projective  equations  can  assume  various 
forms  depending  on  the  choice  of  angular  elements  of  orientation.  We  shall  limit  our 
consideration  here  to  two  specific  forms,  namely  to  the  expressions  given  by  equations 
(2.92)  and  2.98)  in  Volume  1  of  ;he  Third  Edition  of  the  Manual  of  Photcgrammetry.  These 
expressions  are 
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The  orientation  matrix  is  applied  to  rotate  the  direction  cosines  v  defining  a  ray  in 

object  space  into  the  corresponding  direction  cosines  m,  n,  q  in  image  space: 
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In  the  process  of  linearization  of  the  observational  equations,  v/e  shall  require  the  matrix 


(A.o)  SfeSdL 


3m 

3m 

3m 

5a 

3<o 

3x 

Bn 

3 n 

3n 

3a 

3w 

3x 

2i 

sq 

3q 

3a 

3a> 

3x 

When  the  orientation  matrix  is  defined  by  (A. 3)  this  matrix  Is  of  the  form 


i 
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/a  y\  3(m,n,q) 

(  }  db,a>,x) 


rqcosx 


~Cjji  +  By 


n 


/ 


qsin  x  **C,j2+  B'y  -m 

p  -Ffj  +  Ey  0 


in  which 

(A. 8)  p  =  xcosa  +^sino;sina)  -  ysinacos^ . 
Alternatively/  when  (A. 4)  is  used/  it  becomes 


(A  d(m/n/q) 

{  •  }  hicxrurx) 


BX  “ 

am 

-qsinx 
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B'X  “ 

A'm 

-qcosx 
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/ 

EX  - 

dm 

P 

0 

in  which 

(A.  10)  p  =  -xsinasir.co  -  ycosoisinej  +  ycosto* 

At  this  point/  we  shall  assume  that  approximations  a00 ,  ojDO  ,x00  have  been 
employed  in  the  evaluation  of  the  orientation  matrix  and  that  the  resulting  orientation 
matrix  is  used  in  (A. 5)  to  generate  approximate  direction  cosines  m00,  n00,  q00 ,  We 
then  define 

x00  =  c00  (m00/q00} 

(A.  11) 

y°o  =  c00(n00/q00) 

in  which  cc0  denotes  an  appi  amotion  to  the  principal  distance.  The  quantities  x00/y00 
denote  the  computed  values  of  the  plate  coordinates  referred  to  the  principal  point.  The 
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differences  between  (A,  1)  and  (A,  11)  generate  the  observational  equations  for  the 
.adjustment 


f°  «  x  -x00 

(A.  12) 

«S  =  y-r  '  • 


The  elements  of  the  discrepancy  vector  of  the  iinearized  observational  equations  are  given 
simply  by 


(A.  13)  € 


The  elements  of  the  A  matrix  are  provided  by 


(A.  14) 


affiVg) 

?<x°,y°) 


,+£n 

0+J  1+f 

.  s28j 


in  which  the  quantities  £  are  first  order  functions  of  the  approximate  distortion  parameters 
and  are  obtained  from  the  differentiation  of  (A.  1)  with  respect  to  x°/y°.  Inasmuch  as  these 
first  order  quantities  generate  second  o;  Jer  quantities  when  multiplied  by  the  first  order 
variations  in  the  plate  coordinates  (i,e.,  the  residuals  vy,vy),  they  may  ordinarily  be 
neglected. 

To  express  the  elements  of  the  B  matrix  compactly,  we  find  it  convenient  to 
introduce  the  auxiliaries: 


M  =  m00/q00 
(A.  15)  N  n00/q00 
Q  =  c00/q00 
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U  =  Pj00^  +  2x«)  +  2  P”0  x  y 
(A .  16)  V  =  2  P ? 0  x  y  +  P°  0  (i*  +  2y  3) 
W  =  1  -t  P°°  r8 


The  expressions  for  the  elements  of  the  B  matrix  then  become 


(A.  17)  BT  = 
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in  which  the  £  are  the  same  negligible  quantities  as  those  appearing  in  the  A  matrix. 

*  $ 

•  • 

Those  elements  of  the  B  matrix  pertaining  to  the  angular  elements  of  orientation 
can  be  computed  from 
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in  which  the  superscript  'so'  denotes  that  the  elements  of  the  matrix  are  computed  from 
approximations. 


-:13- 


In  cases  where  the  coordinates  of  the  exposure  station  are  to  be  recovered,  the 
«• 

remainder  of  the  B  matrix  can  be  computed  from 


in  which  R00  denotes  the  approximate  distance  between  the  control  point  and  the  exposure 
station,  namely 

R00  -  {CX-(Xc)00]2-»  CY-(Yc)0032  +  CZ-(Zc)00]2  f  . 


APPENDIX  B 

RESULTS  OF  AERIAL  SMAC  CALIBRATIONS  TOR 
KC-6A  CAMERAS  005  AND  002 


In  the  mein  body  of  tho  report  the  results  of  two  Aerial  SMAC  calibrations  of 
Camera  006  were  discussed  and  analyzed  in  considerable  detail.  Though  of  secondary 
interest  in  the  flight  testing  program,  Cameras  005  and  008  were  also  subjected  to  Aerial 
SMAC  calibrations.  Key  results  of  those  calibrations  are  summarized  below  in  Table 
and  Figure 

For  the  most  port,  the  results  for  005  and  008  are  comparable  with  those  obtained 
for  Camera  006.  Thus,  rather  large  displacements  of  the  principal  point,  particularly  in 
xp/  are  experienced,  and  one  sigma  accuracies  of  one  micron  or  better  are  obtained 
for  the  distortion  functions  over  most  of  the  format.  However,  unlike  the  decentering 
distortion  for  Camera  005  which  had  been  found  to  be  unusually  low,  that  for  Camera  005 
(Fig.  20)  turned  out  to  be  decidedly  on  the  high  side,  growing  to  as  much  as  15  microns 
within  the  farmet.  Decentering  distortion  for  Camera  003  (Fig.  20)  was  found  to  be  more 
nearly  typical  of  the  average  modern  mapping  lens. 

Comparisons  between  results  from  Aerial  SMAC  and  results  from  conventional 
multicollimator  calibrations  are  given  in  the  main  body  of  the  report.  There,  too,  are 
given  the  results  of  the  Stellar  SMAC  calibration  of  Camera  008. 
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TABLE  9.  Summary  of  Results  of  Aerial  SMAC  Calibrations  of  Cameras  005  and  008. 

r - - - - - - 

I  CAMERA  005  (Test  No.  2)  CAMERA  008  (Test  No.  1) 


Parameter 
xp  (mm) 

Xp  (mm) 

c  (mm) 

Kj  (mm/mm3) 
K2  (mm/mmB) 

Jj  (mm/mm2) 

(deg) 


Value 

.1674.004 
-.045*. 004 
152.5394.002 


Parameter 

xp  (mm) 

Yv  (mm) 
c  (mm) 


(~.  I8ftfc. 018)x  10  Kj  (mm/mm3) 

(  ,503±.079)x  10"  “  Kg  (mm/mm6) 

(6 . 80^4 . 434)  x  10  Jx  (mm/mm2) 

142.243.7  $0  (deg) 


.1754.003 

-.0524.003 

154.6264.002 

(-.385±.015)x  W~l 
(1.34/4. 063)  x  10"  J 

(2. 5504. 359)  x  10~? 
19.747.7 


POINTS  ON  RADIAL  AND  DECENTERING  DISTORTION  CURVES 


CAMERA  005 

~6^> 

(Microns) 

0 

■2.40.01 
4.040.11 
4.2-40.20 
0,67 
-0.  iil.03 
■!|.j±1.31 
~i34l. 41  1 


Pr 

(Microns) 

0 

0.2740.02 
1- 0940.07 
2.4540. 16 
4. 3640.28 
6.840.43 
9*804-0.63 
13*34-40.98 


(u)  Mjfei»p#i  to  vdlua  of  e  of  ISi.262nun 
(b)  refetttd  |q  ya|ue  af  c  of  151,255mm 


CAMERA  008 


(Mic.ons) 

0 

4.040.01 
6.40.09 
5.840.27 
2.  40.56 
-3.340.88 
-6.94l.12 
-2*741.20 


(Microns) 


0.  RJ40.C1 
O.4I40. 06 
0.92-0. 13 
1.63±0.i3 
2.5540.36 
3. 6/4  0.52 
5,0040.70 
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